| (861-862) 


CZASOPISMO 
POLSKIEGO TOWARZYSTWA CHEMICZNEGO 


Praca dotowana ze środków Uniwersytetu Wrocławskiego 


RADA REDAKCYJNA 
RYSZARD ADAMIAK, IRENA BARANOWSKA, ANDRZEJ BARAŃSKI, 
BOGUSŁAW BUSZEWSKI (PRZEWODNICZĄCY), TADEUSZ GORECKI, 
MIETEK JARONIEC, ANATOL KOJŁO, TADEUSZ M. KRYGOWSKI, 

JERZY LESZCZYNSKI, KRZYSZTOF MATYJASZEWSKI, PIOTR PANETH, 
JANUSZ PAWLISZYN, K. MICHAŁ PIETRUSEWICZ, DARIUSZ POGOCKI, 
MAREK POTRZEBOWSKI, SŁAWOMIR RUBINSZTAJN, GRZEGORZ SCHROEDER, 
ANDRZEJ W. SOKALSKI, ARTUR P. TERZYK 


KOMITET REDAKCYJNY 


MARCIN DRĄG, ADAM JEZIERSKI, LESZEK KĘPIŃSKI, 
LUDWIK KOMOROWSKI, WITOLD RYBA-ROMANOWSKI, SŁAWOMIR SZAFERT, 
ANDRZEJ TROCHIMCZUK, KAZIMIERA WILK 


REDAKTOR NACZELNY 
ZDZISŁAW LATAJKA 


SEKRETARZ REDAKCJI 
BEATA ŚWIĄTEK-TRAN 
e-mail: beata.swiatek-tran@chem.uni.wroc.pl 
DAGMARA JACEWICZ 


e-mail: czasopisma@ptchem.pl 


Korespondencję należy kierować pod adresem: 
Redakcja „Wiadomości Chemicznych” 
ul. F. Joliot-Curie 14, 50-383 Wrocław 
e-mail: czasopisma@ptchem.pl 


© Copyright by Redakcja „Wiadomości Chemicznych”, Wrocław 2019 
pISSN 0043-5104 
eISSN 2300-0295 


Przygotowanie do druku i druk: 
Firma Wydawnicza K2, e-mail: k2Qdruk-k2.pl 


WIADOMOŚCI 2019.73, 3-4 
chemiczne PL ISSN 0043-5104 


W OSIEMDZIESIĄTĄ ROCZNICĘ URODZIN 
PROFESORA ALEKSANDRA KOLLA 


Aleksander Koll urodził się w Wilnie 3 czerwca 
1939 roku. Po wojnie rodzina przeniosła się na „zie- 
mie odzyskane” i zamieszkała w Głogowie. Tutaj 
ukończył szkołę podstawową i średnią, a w roku 1957 
rozpoczął studia na Uniwersytecie Wrocławskim 
na Wydziale Matematyki, Fizyki i Chemii. Studia 
ukończył w 1962 roku przedkładając pracę magister- 
ską wykonaną pod opieką prof. Lucjana Sobczyka. 
W tym samym roku rozpoczął studia doktoranckie. 
Stopień doktora uzyskał w 1968 roku na podstawie 
dysertacji zatytułowanej „Struktura i właściwości 
związków (pirydyl),Y”. Promotorem pracy był profesor Lucjan Sobczyk. Kolejne 
szczeble kariery naukowej to habilitacja w 1986 roku, stanowisko docenta w 1988 
roku, profesora nadzwyczajnego w 1992 roku i tytuł profesora w roku 1996. 

W swojej działalności naukowej Profesor Koll interesuje się różnorodnymi 
aspektami oddziaływań międzycząsteczkowych, a zwłaszcza wiązania wodorowego. 
Zainteresowanie tą tematyką rozpoczęło się już u początku Jego drogi naukowej. 
Praca magisterska dotyczyła polarności wiązań wodorowych między fenolami 
i aminami, praca doktorska poświęcona była korelacjom między strukturą i wła- 
ściwościami bipirydyli, zaś głównym wątkiem habilitacji była polarność i struktura 
wewnątrzcząsteczkowych wiązań wodorowych. W swoich badaniach Profesor Koll 
wykorzystuje zarówno zaawansowane techniki doświadczalne (pomiary momentów 
dipolowych, badania spektroskopii IR, UV-Vis, NMR) jak i obliczenia teoretyczne. 
Interesuje się zarówno wiązaniami wodorowymi międzycząsteczkowymi, jak 
i wewnątrzcząsteczkowymi. Silne międzycząsteczkowe wiązania wodorowe badał na 
przykład w kompleksach halogeno-wodorów i halogeno-deuterów z akceptorami 
wodoru/deuteru zawierającymi tlen. Bardzo ciekawe rezultaty przyniosły badania 
wewnątrzcząsteczkowych wiązań wodorowych w zasadach Schiffa i Mannicha pro- 
wadzone wspólnie z prof. Sobczykiem, prof. Rospenk, prof. Filarowskim i innymi. 
Zarówno rezultaty doświadczalne, jak i obliczenia teoretyczne pozwoliły powiązać 
właściwości i parametry termodynamiczne wewnątrzcząsteczkowych wiązań wodo- 
rowych ze strukturą i konformacyjną cząsteczek. Do ważnych osiągnięć należy udo- 
wodnienie, że charakterystyczna zależność polarności wiązania wodorowego od 
DpKa dotyczy nie tylko międzycząsteczkowych wiązań wodorowych, ale i wiązań 
wewnątrzcząsteczkowych. W swoich badaniach wiele uwagi poświęca efektom ste- 
rycznym i rezonansowym, na przykład w układach podstawionych bifenyli. 
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Profesor Koll prowadzi swoje badania we współpracy z licznym gronem 
naukowców z Polski i z zagranicy. W czasie stażu odbytego w Kanadzie współpra- 
cował z prof. Fyfem zajmując się badaniami NMR i kinetycznymi aspektami pod- 
stawienia nukleofilowego. W czasie pobytu w na Uniwersytecie w Leuven (Belgia) 
nawiązał kontakty naukowe z prof. Huyskensem, prof. Zeeger-Huyslens i dr Helle- 
mansem. Szczególnie żywe kontakty naukowe wiążą profesora Kolla z prof. Wol- 
schannem i prof. Karpfenem z Uniwersytetu Wiedeńskiego, z którymi prowadzi 
szeroko zakrojone badania teoretyczne i doświadczalne wewnątrzcząsteczkowych 
wiązań wodorowych. Ożywione i owocne kontakty naukowe wiążą Profesora Kolla 
z Uniwersytetami w St. Petersburgu i Tjumeniu. Współpraca z profesorami, Deni- 
sowem, Shchepkinem, Bułaninem, Burejką, Rutkowskim i Melikową zaowocowała 
cyklem fundamentalnych prac na temat wiązania wodorowego i oddziaływań mię- 
dzycząsteczkowych, a wśród nich unikalne badania widm oscylacyjnych prowa- 
dzone w skroplonych gazach szlachetnych. Dorobek publikacyjny profesora Kolla 
to ponad 200 oryginalnych prac opublikowanych w czasopismach o zasięgu mię- 
dzynarodowym, które cytowane były około 3000 razy, indeks Hirscha wynosi 31. 

Bardzo bogata jest działalność dydaktyczna i organizacyjna Profesora Kolla. Do 
przejścia na emeryturę, w roku 2009 wypromował on 11 doktorów i kilkudziesięciu 
magistrów. Prowadził wiele wykładów kursowych i monograficznych, jest współ- 
autorem podręcznika Eksperymentalna Chemia Fizyczna. W latach 1996-2002 był 
prodziekanem na Wydziale Chemii Uniwersytetu Wrocławskiego, od 1997 do 2009 
kierował Zespołem Badawczym Struktury i Oddziaływań Międzycząsteczkowych, 
przez wiele lat sprawował funkcje kierownika Zespołu Dydaktycznego Chemii 
Fizycznej. W latach 1998-2003 był przewodniczącym Sekcji Fizykochemii Związ- 
ków Organicznych Polskiego Towarzystwa Chemicznego. Wielokrotnie organizował 
i współorganizował doroczne Szkoły Fizykochemii Organicznej zapoczątkowane 
przez Profesora Sobczyka, a następnie kontynuowane w ramach Sekcji Fizykoche- 
mii Związków Organicznych. 

Profesor Koll był wielokrotnie nagradzany i odznaczany, w roku 1966 otrzy- 
mał Nagrodę Polskiego Towarzystwa Chemicznego, Nagrody Ministra w 1969, 1983 
i 1990, wielokrotnie otrzymywał nagrody Rektora Uniwersytetu Wrocławskiego, 
w roku 2007 otrzymał Medal Zawidzkiego za wybitne osiągnięcia z zakresu chemii 
fizycznej i nieorganicznej. 

Od 2009 roku sprawuje funkcję Rektora Niepublicznej Wyższej Szkoły Medycz- 
nej we Wrocławiu. Jest redaktorem czasopisma Kosmetologia Estetyczna. 
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ABSTRACT 


The structures and vibrational frequencies of a series of complexes of benzal- 
dehyde, N-benzylideneamine and N-benzylidenemethylamine with simple hydro- 
gen bond donors and with dihalogens have been investigated at the M06-2X/6- 
311++G(2d,2p) level. It is shown that the well-known blue shift of formyl C-H 
stretching frequencies encountered upon complex formation is a direct consequ- 
ence of the internal coupling of the formyl C=O bond with the vicinal C-H degree of 
freedom of the free aldehyde and of the imine C=N bond with the vicinal C-H bond 
in the two imine derivatives. These blue shifts occur without directly involving the 
CH bonds in hydrogen or halogen bonding. 


Badano technikami obliczeniowymi na poziomie M06-2X/6-311++G(2d,2p) 
strukturę i częstotliwości drgań w seriach kompleksów aldehydu benzoesowego, 
N-benzylidenoaminy i N-benzylidenometyloaminy z prostymi donorami protonu 
i ze związkami dihalogenowymi. Pokazano, że dobrze znane krótkofalowe prze- 
sunięcie drgania rozciągającego C-H grupy formylowej, obserwowane w wyniku 
utworzenia kompleksu, jest bezpośrednią konsekwencją wewnętrznego oddziały- 
wania wiązania C=O grupy formylowej z wicynalnym wiązaniem C-H w dwóch 
pochodnych iminowych. Wspomniane przesunięcie występuje mimo tego, że grupa 
C-H nie jest bezpośrednio zaangażowana w wiązanie wodorowe lub halogenowe. 


Keywords: blue-shifted hydrogen bonding, halogen bond, density functionals, ben- 
zaldehyde derivatives 

Słowa kluczowe: wiązanie wodorowe typu „blue-shifted, wiązanie halogenowe, 
funkcjonały gęstości, pochodne aldehydu benzoesowego 
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INTRODUCTION 


Weak hydrogen bonds A-H---X, involving less polar A-H groups, like e.g. C-H, 
are the object of systematic experimental and intense theoretical studies. Some of 
these A-H---X hydrogen bonds, in the majority C-H---X contacts, exhibit peculiar 
properties such as a contraction of the C-H bond, a blue shift of the C-H stretching 
mode, and a decrease of the infrared intensity [1-8]. These properties contrast with 
the well-established features of conventional hydrogen bonding. Therefore, C-H---X 
hydrogen bonds displaying this behavior are often termed as blue-shifting hydrogen 
bonds. 

Blue shifts of A-H or C-H stretching frequencies may, however, also occur when 
the donor groups are not directly involved in hydrogen bonding but are in a vicinal 
position to a A-Y or C-Y hydrogen bond acceptor, e.g. C=O or C-Hal or related 
electron-withdrawing groups [9-13]. For these cases the name blue-shifted has been 
coined. Undoubtedly, red-shifting as well as blue-shifting or blue-shifted hydrogen 
bonding can be rationalized by the conventional theory of intermolecular interac- 
tion. Since the latter are mostly encountered in cases of weak or modestly strong 
intermolecular complexes, the most promising and at the same time also predictive 
explanations for the understanding of blue vs red shifts are based on an analysis 
of intramolecular properties of the molecule that displays the frequency shifts in 
the complex rather than on a posteriori wave function analysis of the already for- 
med complex. The origin of the blue shift must have intramolecular roots. Attempts 
in this direction have been to investigate the structural response to external fields 
[12, 15-17], the sign of the dipole moment derivative, du/dr(X-H) [2, 16], and more 
directly connected to the potential energy surface (PES), the coupling between vici- 
nal A-H and Y-A bonds as calculated either with the aid of selected scans of A-H 
and Y-A bonds [9, 14], or equivalently by conventional normal coordinate analysis 
[17-21]. 

In previous calculations the phenomenon of blue-shifted hydrogen bonding 
has been investigated in complexes with formaldehyde and with substituted alde- 
hydes [11, 12, 19]. In this work, we chose simple complexes of benzaldehyde, 1, 
and of two of its imine derivatives, the Schiff bases N-benzylideneamine, 2, and 
N-benzylidenemethylamine, 3, (Fig. 1) with the hydrogen bond donors HF and H,O 
and with the halogen bond donors Cl, and CIF. Intra- and intermolecular hydrogen 
of phenol to 2 and 3 have been studied experimentally [22]. To the best knowledge 
of the authors there are neither experimental nor theoretical studies available that 
discuss the blue-shifted properties of the formyl C-H or nitrogen-analogous N-H 
bond stretchings in these complexes. 
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Figure 1. M06-2X/6-311++G(2d,2p) optimized structures of the molecules benzaldehyde, N-benzylidene- 
amine and N-benzylidenemethylamine 

Rysunek 1. | Zoptymalizowane na poziomie M06-2X/6-311++G(2d,2p) struktury molekularne aldehydu ben- 
zoesowego, N-benzylidenoaminy i N-benzylidenometyloaminy 


1. COMPUTATIONAL METHODOLOGY 


All calculations on the isolated molecules and complexes were performed 
within the M06-2X [23] density functional method using the 6-311++G(2d,2p) as 
implemented in the GAUSSIAN 09 suite of programs [24]. This computational level 
is certainly sufficient for our purposes and should yield reliable results for structu- 
res, energies and harmonic vibrational frequencies of the molecules and comple- 
xes. Geometry optimizations and frequency calculations were performed without 
explicit counterpoise corrections (CP) to the basis set superposition error (BSSE). 
Harmonic force constant analyses in internal coordinates were carried out with the 
program GAR2PED [25]. 
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2. RESULTS AND DISCUSSION 


2.1. STRUCTURES AND STABILIZATION ENERGIES 


The optimized structures of the isolated molecules benzaldehyde, N-ben- 
zylideneamine and N-benzylidenemethylamine, as obtained at the M06-2X/6- 
311++G(2d,2p) level are shown in Figure 1. From the structural as well as the vibra- 
tional spectroscopic point of view the most characteristic property of aldehydes is 
the long formyl C-H bond. Concomitantly, the formyl C-H stretching frequency is 
found at substantially lower wavenumbers than e.g. the aromatic C-H stretching fre- 
quencies. A second characteristic feature of the formyl group is the strong coupling 
between C=O and C-H stretching degrees of freedom. It has been shown in a large 
number of theoretical investigations on aldehydes [11, 12, 18] and on a variety of 
other molecules [13, 17, 19] that one of the key properties that allows to predict the 
appearance of blue shifts of C-H, or more generally of A-H stretching frequencies in 
intermolecular complexes, is the positive sign of the intramolecular coupling con- 
stant between C-H and C=O bonds of the formyl group. In Table 1, the off-diagonal 
harmonic force constants f(C-H, C=X) and f(C-C,C=X), (X=O,N) for 1, 2, and 3 
are listed. All are positive. Interestingly, f(C-C, C=X) is by about a factor of two 
larger. The largest couplings occur for 1. An equivalent way to show the mutual 
interdependence of C-H or C-C distances with the C=O or C=N degrees of freedom 
can be achieved by scanning C=O or C=N and plotting the optimized C-H and 
C-C bond lengths. This has been done in Figure 2. The negative slope indicates that 
r(C-H) shrinks, when r(C=O) or r(C=N) widens. The same applies for r(C-C) with 
a distinctly steeper slope. 


Table 1. Harmonic off-diagonal force constants f(C-H, C=X) and f(C-C,C=X)' 
Tabela 1. Pozadiagonalne harmoniczne stałe siłowe f(C-H, C=X) i f(C-C,C=X)' 
Compound f(C-H,C=X)* f(C-C,C=X) 

Benzaldehyde 0,59 0,99 
N-benzylideneamine 0,30 0,74 
N-benzylidenemethylamine 0,33 0,71 

1 X=0,N. 

* All values in mdyn Å”. 


When the C=O bond of benzaldehyde acts as acceptor, two different minima 
have been obtained for the complexes with each of the four donor molecules. For the 
two derivatives with C=N bonds only one minimum per donor has been localized. 
The resulting complex structures are displayed in Figure 3. 

The shortest intermolecular distances and selected intramolecular bond 
distance changes are reported in Tables 2 and 3. In all sixteen complexes, we observe 
an increase of C=O or C=N distances, as expected larger for HF than for H,O, and 


BLUE-SHIFTED FORMYL C-H STRETCHING FREQUENCIES 171 


larger for CIF than for Cl,, a consequence of the stronger hydrogen bond donor 
properties of HF and the better halogen bond donor properties of CIE. 


Table 2. M06-2X/6-311++G(2d,2p) calculated intermolecular bond distances (A), and intramolecular bond 
distance changes (mA) in hydrogen-bonded complexes of benzaldehyde, N-benzylideneamine 
and N-benzylidenemethylamine with hydrogen bond donors H,O and HF 

Tabela 2. Obliczone na poziomie M06-2X/6-311++G(2d,2p) długości wiązań miedzyczasteczkowych (A) 
i zmiany długości wiązań wewnątrzcząsteczkowych (mA) w kompleksach z wiązaniem wodoro- 
wym złożonych z aldehydu benzoesowego, N-benzylidenoaminy lub N-benzylidenometyloaminy, 
i donora wiązania wodorowego: H,O lub HF 

Benzaldehyde r(O---H) Ar(H-X)1 Ar(C-C) Ar(C-H) Ar(C=O) 
+H,O, la 1,93 7,8 -6,7 -1,8 4,9 
+H,O, 1b 1,95 8,1 -5,4 -3,0 6,7 
+HE, Ic 1,67 21,3 -11,0 -4,0 8,5 
+HE, 1d 1,66 22,6 -12,3 -4,0 10,6 
N-benzylideneamine r(N---H) Ar(H-X) Ar(C-C) Ar(C-H) Ar(C=N) 
+H,O, ża 1,94 13,5 -4,1 -1,2 2,7 
+HE, 2b 1,62 42,7 -7,5 -2,6 3,5 
N-benzylidenemethylamine r(N---H) Ar(H-X) Ar(C-C) Ar(C-H) Ar(C=N) 
+H,0O, 3a 1,92 14,5 -2,6 -1,5 2,3 
+HĘ 3b 1,60 47,8 -6,4 -3,4 3,7 
! X=O,F. 
Table 3. M06-2X/6-311++G(2d,2p) calculated intermolecular bond distances (A), and intramolecular 


bond distance changes (mA) in halogen-bonded complexes of benzaldehyde, N-benzylideneamine 
and N-benzylidenemethylamine with halogen bond donors Cl, and CIF 


Tabela 3. Obliczone na poziomie M06-2X/6-311++G(2d,2p) długości wiązań międzycząsteczkowych (A) 
i zmiany długości wiązań wewnątrzcząsteczkowych (mA) w kompleksach z wiązaniem halogeno- 
wym złożonych z aldehydu benzoesowego, N-benzylidenoaminy lub N-benzylidenometyloaminy 
i donora wiązania halogenowego: Cl, lub CIF 

Benzaldehyde r(O---CI) Ar(Y-Cl)' Ar(C-C) Ar(C-H) Ar(C=O) 
"CL. le 2,73 2,02 -4,2 -1,3 3,0 
+CL, 1f 2,66 2,02 -5,5 -1,2 4,1 
+CIĘ 1g 2,53 1,65 -7,2 -2,2 5,0 
+CIĘ lh 2,47 1,66 -9,0 -2,6 6,9 
N-benzylideneamine r(N---CI) Ar(Y-Cl) Ar(C-C) Ar(C-H) Ar(C=N) 
+CL, 2c 2,66 23,4 -3,7 -0,9 1,1 
+CIĘ 2d 2,26 64,5 -9,8 -2,3 2,9 
N-benzylidenemethylamine r(N--Cl) Ar(Y-Cl) Ar(C-C) Ar(C-H) Ar(C=N) 
+CL, 3c 2,63 27,9 -2,9 -1,3 1,4 
+CIĘ 3d 2,20 80,3 -7,7 -3,9 5,2 


'Y=CLE 
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We also calculated that Ar(C-C) and Ar(C-H) are negative, in line with the sign of 
the coupling force constants in the isolated molecules (Fig. 1) and with the intra- 
molecular scans (Fig. 2). Moreover, the reduction of r(C-C) is distinctly larger than 
that of r(C-H) in agreement with the relative sizes of the numerical derivatives of 
Ar(C-C)/Ar(C=X) and Ar(C-H)/Ar(C=X). 


1.125 


1.115 
T 
© 1.10577 
D 

1.095 

1.085 

1.16 1.21 1.26 1.31 1.36 
R(C=X) 


1.16 1.21 1.26 1.31 1.36 
R(C=X) 


Figure 2. M06-2X/6-311++G(2d,2p) optimized A-H distances and C-C distances as obtained upon scan- 
ning the C=O distance in benzaldehyde (©) and the C=N distance in N-benzylideneamine (CH 
and N-benzylidenemethylamine (A). All distances in A 

Rysunek 2. Zoptymalizowane na poziomie M06-2X/6-311++G(2d,2p) odległości A-H i C-C uzyskane 
przy skanowaniu długości wiązania C=O w aldehydzie benzoesowym (©) oraz wiązania C=N 
w N-benzylidenoaminie (L) i N-benzylidenometyloaminie (A). Wszystkie wartości w A 
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Figure 3. M06-2X/6-311++G(2d,2p) optimized A-H structures of complexes of benzaldehyde, 
N-benzylideneamine and N-benzylidenemethylamine with the hydrogen bond donors HF and 
H,O and the halogen bond donors Cl, and CIF 

Rysunek 3. Zoptymalizowane na poziomie M06-2X/6-311++G(2d,2p) struktury molekularne kompleksów 
złożonych z aldehydu benzoesowego, N-benzylidenoaminy lub N-benzylidenometyloaminy 
i donora wiązania wodorowego: H,O lub HE, albo donora wiązania halogenowego: Cl, lub CIF 


"p 


The calculated stabilization energies of the complexes are compiled in Tables 4 
and 5. In agreement with comparable hydrogen-bonded dimers, the complexes with 
HF are substantially more stable than those with H,O as interaction partner. Simi- 
larly, the halogen-bonded complexes with CIF are much stronger than those with 
Cl 
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Table 4. M06-2X/6-311++G(2d,2p) calculated bare AE, and zero-point energy corrected binding ener- 
gies, AE(ZPE), in hydrogen-bonded complexes of benzaldehyde, N-benzylideneamine and 
N-benzylidenemethylamine with hydrogen bond donors H,O and HE All values in kcal mof" 

Tabela 4. Obliczone na poziomie M06-2X/6-311++G(2d,2p) energie elektronowe, AE, oraz energie skory- 
gowane o energię drgań zerowych, AE(ZPE), w kompleksach z wiązaniem wodorowym złożonych 
z aldehydu benzoesowego, N-benzylidenoaminy lub N-benzylidenometyloaminy i donora wiąza- 
nia wodorowego: H,O lub HF 

Benzaldehyde AE AE(ZPE) 
+H.O, la -6,6 -5,0 
+H,O, Ib -6,5 -4,8 
+HE, Ic -10,5 -8,2 
+HE, 1d -10,3 -8,3 
N-benzylideneamine AE AE(ZPE) 
+H,O, 2a -7,5 -5,4 
+HE, 2b -14,0 -11,5 
N-benzylidenemethylamine AE AE(ZPE) 
+H,0, 3a -7,4 -5,4 
+HE, 3b -14,2 -11,8 


Wszystkie wartości w kcal mol. 


Table 5. M06-2X/6-311++G(2d,2p) calculated bare AE, and zero-point energy corrected binding ener- 
gies, AE(ZPE), in halogen-bonded complexes of benzaldehyde, N-benzylideneamine and 
N-benzylidenemethylamine with halogen bond donors CL, and CIF. All values in kcal mol” 

Tabela 5. Obliczone na poziomie M06-2X/6-311++G(2d,2p) energie elektronowe, AE, oraz energie skory- 
gowane o energię drgań zerowych, AE(ZPE), w kompleksach z wiązaniem halogenowym z ald- 
ehydu benzoesowego, N-benzylidenoaminy i N-benzylidenometyloaminy lub donora wiązania 
halogenowego: Cl, lub CIF 

Benzaldehyde AE AE(ZPE) 
+CI, le -4,2 -3,6 
+Cl, If -4,8 -4,2 
+CIĘ 1g -6,5 -5,7 
+CIĘ lh -7,8 -6,4 
N-benzylideneamine AE AE(ZPE) 
+Cl, 2c -6,0 -5,2 
+CIĘ 2d 10,8 -9,5 
N-benzylidenemethylamine AE AE(ZPE) 
+Cl, 3c -5,6 -6,0 
+CIĘ 3d -12,4 -10,8 


Wszystkie wartości w kcal mol. 
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2.2. SHIFTS OF STRETCHING FREQUENCIES 


The calculated shifts of the intramolecular stretching modes of hydrogen bond 
and halogen bond donors and of the C-H, C-C and C=O (C=N) stretching modes of 
the acceptors are reported in Tables 6 and 7. Turning first to the complexes 1a to 1d, 
we encounter the expected behavior for hydrogen bonding to the carbonyl group 
of aldehydes: sizable red shifts of the hydrogen bond donor stretches, red shifts of 
the carbonyl stretching frequencies, and blue shifts of the formyl C-H stretches, the 
latter are not involved in hydrogen bonding. The C-C stretching frequencies are 
more difficult to interpret because the C-C degree of freedom cannot be described 
by a single normal mode but takes part in several energetically far separated normal 
modes. We have selected the three normal modes with the largest contributions of 
C-C stretch. Remarkably, all three show a small blue shift. 


Table 6. M06-2X/6-311++G(2d,2p) calculated harmonic stretching frequency shifts in hydrogen-bonded 
complexes of benzaldehyde, N-benzylideneamine and N-benzylidenemethylamine with hydrogen 
bond donors H,O and HF. All values in cm" 

Tabela 6. Obliczone na poziomie M06-2X/6-311++G(2d,2p) przesunięcia położeń harmonicznych drgań 
rozciągających w kompleksach z wiązaniem wodorowym złożonych z aldehydu benzoesowego, 
N-benzylidenoaminy lub N-benzylidenometyloaminy i donora wiązania wodorowego: H,O lub 
HE Wszystkie wartości w cm” 


Benzaldehyde Av(H-X)' Av(C-C) Av(C-H) Av(C=O) 
+H,O, la -103 2,4,5 17 -17 
+H,O, 1b -114 6, 6, 2 40 -25 
+HĘ 1c -510 2,6,5 43 -31 
+HF, 1d -530 12,2 50 -40 
N-benzylideneamine Av(H-X) Av(C-C) Av(C-H) Av(C=N) 
+H,O, 2a -345 6, 6, 6 16 -3 
+HE, 2b -992 1, 14,5 34 -1 
N-benzylidenemethylamine Av(H-X) Av(C-C) Av(C-H) Av(C=N) 
+H,O, 3a -278 LS 21 -1 
+HE, 3b -1078 2525.7 36 1 


'X=0, F 
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Table 7. M06-2X/6-311++G(2d,2p) calculated harmonic stretching frequency shifts in halogen-bonded 
complexes of benzaldehyde, N-benzylideneamine and N-benzylidenemethylamine with halogen 
bond donors Cl, and CIE All values in cm™ 

Tabela 7. Obliczone na poziomie M06-2X/6-311++G(2d,2p) przesunięcia położeń harmonicznych drgań 
rozciągających w kompleksach z wiązaniem halogenowym złożonych z aldehydu benzoesowego, 
N-benzylidenoaminy lub N-benzylidenometyloaminy i donora wiązania halogenowego: CI, lub 
CIE Wszystkie wartości w cm" 


Benzaldehyde Av(Y-Cl)' Av(C-C) Av(C-H) Av(C=O) 
+CL, le -9 11,3 12 -12 
+Cl,, 1f -23 5,6,2 10 -17 
+CIĘ 1g -39 3,5,4 23 -20 
+CIĘ lh -53 9:9;.2 26 -29 
N-benzylideneamine Av(Y-Cl) Av(C-C) Av(C-H) Av(C=N) 
+CL, 2c -32 3,253 13 -2 
+CIĘ 2d -151 11, 11,5 32 -6 
N-benzylidenemethylamine Av(Y-Cl) Av(C-C) Av(C-H) Av(C=N) 
+CL, 3c -56 6, 12, 4 14 -2 
+CIĘ 3d -197 5, 1,3 48 -9 
'Y=Cl, E 


In the hydrogen bonded complexes 2a, 2b, 3a, and 3b, the red shifts of the H-X 
stretching frequencies are much larger than in benzaldehyde complexes, in line with 
the larger stabilization energies. The C-H and C-C blue shifts are comparable to the 
benzaldehyde case, the C=N shifts are, however, close to zero. 

For the halogen-bonded complexes, le to 1h, and 2c, 2d, 3c, and 3d, we cal- 
culate the same patterns for the shifts of the stretching frequencies: red shifts of the 
dihalogens, blue shifts of the C-H and C-C stretching frequencies. In addition, the 
C=N stretches are now also red-shifted. 


SUMMARY AND CONCLUSIONS 


In the present work, we have investigated the structures and harmonic vibra- 
tional frequencies of complexes of benzaldehyde, N-benzylideneamine and N-ben- 
zylidenemethylamine with the hydrogen bond donors H,O and HĘ and with the 
halogen bond donors Cl, and CIF. Hydrogen and halogen bonding to the formyl 
C=O or to the imine C=N lead to widening of these acceptor groups and to contrac- 
tions of vicinal C-H and C-C bonds. We have shown that knowledge of the intramo- 
lecular coupling between these bonds, as gained by appropriate scans of the relevant 
sections of the intramolecular potential surfaces and/or by conventional normal 
coordinate analysis of the acceptor molecules, is a useful means to predict and to 
rationalize the sign and the size of structural distortions in the complexes as well as 
the frequency shifts encountered in the complexes. In addition to the red shifts of 
hydrogen and halogen bond donors and of the acceptor groups, blue shifts are cal- 
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culated not only for the C-H, but also for the C-C stretching frequencies, although 
these vicinal groups are not directly involved in hydrogen or halogen bonding. 
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miki molekularnej. 
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ABSTRACT 


The bridged proton dynamics has been an object of theoretical investigations 
in selected aromatic Mannich base-type compounds. Here, we would like to sum- 
marize results obtained on the basis of Car-Parrinello molecular dynamics (CPMD) 
and path integral molecular dynamics (PIMD) methods. The simulations were 
performed in the gas phase (for isolated molecule), in the crystalline phase (based 
on experimentally measured crystal structures) and in solution (where CCl, was 
used as a solvent). Therefore, the environmental influence on the bridged proton 
dynamics was taken into account during our studies. We mainly focused on time- 
-evolution of metric parameters (interatomic distances) of intramolecular hydrogen 
bonds analyzing the CPMD and PIMD trajectories. Another aspect of our studies 
was devoted to the reproduction of spectroscopic parameters based on the CPMD 
trajectories. Our attention was paid to computations of infrared spectra (IR) using 
atomic velocity and dipole moments obtained during the CPMD simulations. In 
addition, the “snapshot-envelope” method was applied to enhance the spectroscopic 
studies. The review article will close with the discussion of 1D and 2D potential 
of mean force (pmf) computed for selected Mannich bases. The computationally 
obtained results were verified using experimental data available. 


Keywords: CPMD, PIMD, potential of mean force, intramolecular hydrogen bond, 
Mannich base derivatives 

Słowa kluczowe: CPMD, PIMD, potencjał siły średniej, wewnątrzcząsteczkowe wią- 
zanie wodorowe, pochodne zasad Mannicha 
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WWW 
RAHB 
CPMD 
PIMD 
pmf 


1D 
2D 


WYKAZ STOSOWANYCH SKRÓTÓW 


wewnątrzcząsteczkowe wiązanie wodorowe (ang. 
intramolecular hydrogen bond) 

wiązanie wodorowe wspomagane rezonansem (ang. 
Resonance Assisted Hydrogen Bond) 

dynamika molekularna Cara-Parrinello (ang. Car- 
-Parrinello molecular dynamics) 

dynamika molekularna w formalizmie całek po tra- 
jektoriach (ang. Path Integral Molecular Dynamics) 
potencjał siły średniej (ang. potential of mean force) 
jednowymiarowy (ang. one-dimensional) 
dwuwymiarowy (ang. two-dimensional) 
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WPROWADZENIE 


Dynamika protonu w mostku wodorowym w wybranych, aromatycznych zasa- 
dach Mannicha stała się obiektem naszych badań teoretycznych na przestrzeni kilku 
lat. Zasady Mannicha, podobnie jak zasady Schiffa, są uznawane za modelowe układy 
w badaniach nad wewnątrzcząsteczkowym wiązaniem wodorowym (WWW) [1-4]. 
Obecność takiego wiązania daje możliwość wewnętrznej reorganizacji w molekule 
(np. obecność zjawiska przeniesienia protonu), a co z tym jest związane wpływa 
na jej właściwości fizyko-chemiczne. Dodatkowo, tworzy się tzw. pierścień quasi- 
-aromatyczny, którego obecność powoduje dodatkową stabilizację konformacyjną 
molekuły i reorganizację struktury elektronowej. Właściwości fizyko-chemiczne 
mostków wodorowych są modulowane poprzez efekty steryczne i/lub indukcyjne 
[5] wywołane wprowadzeniem podstawników w „sąsiedztwo” grupy iminowej, bądź 
też do pierścienia fenylowego. Nie bez znaczenia jest też sprzężenie 7r-elektronowe, 
pojawiające się wraz z obecnością pierścienia aromatycznego w molekule [6, 7]. 
Efekt ten może w sprzyjających warunkach prowadzić do znacznej modyfikacji cha- 
rakteru WWW, co określane jest mianem wiązań wodorowych wspomaganych rezo- 
nansem (RAHB). W obecnym artykule chcielibyśmy podsumować rezultaty badań 
teoretycznych, wykonanych dla wybranych zasad Mannicha, w oparciu o metodę 
dynamiki molekularnej Cara-Parrinello (CPMD) [8], a także dynamikę moleku- 
larną w formalizmie całek po trajektoriach (PIMD) [9-11]. Symulacje były prowa- 
dzone w fazie gazowej (dla pojedynczej molekuły), w fazie krystalicznej (w oparciu 
o zmierzone eksperymentalnie struktury krystaliczne) i w rozpuszczalniku (w nie- 
polarnym CC1,). A zatem wpływ otoczenia na dynamikę mostka wodorowego został 
w badaniach kwantowo-chemicznych uwzględniony w szerokim zakresie. Zasto- 
sowanie w badaniach metod dynamiki molekularnej ab initio (CPMD i PIMD), 
pozwoliło na prześledzenie zmian parametrów metrycznych WWW w zależności 
od czasu symulacji. Innym istotnym aspektem badań była reprodukcja parametrów 
spektroskopowych (zwłaszcza widma podczerwieni (IR)) w oparciu o analizy tra- 
jektorii otrzymanych z symulacji metodą CPMD. W tym celu zastosowane zostały 
trzy podejścia metodologiczne: 

I. Otrzymanie widm mocy z prędkości atomowych 

II. Otrzymanie widm IR z trajektorii momentów dipolowych 

III. Zastosowanie metody “snapshot-envelope” [12], która umożliwiła skwanto- 
wanie drgania O-H. 

Artykuł zakończymy dyskusją rezultatów otrzymanych 1D- i 2D-potencjałów 
siły średniej (pmf) wyznaczonych dla wybranych zasad Mannicha. Rezultaty symu- 
lacji komputerowych zostały porównane z danymi eksperymentalnymi. 
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1. ZASTOSOWANIE METODY DYNAMIKI CARA-PARRINELLO 
(CPMD) I DYNAMIKI MOLEKULARNEJ W FORMALIZMIE 
CAŁEK PO TRAJEKTORIACH (PIMD) W BADANIACH 
WEWNĄTRZCZĄSTECZKOWYCH WIĄZAŃ WODOROWYCH 
W WYBRANYCH POCHODNYCH ZASAD MANNICHA 


Omawiane pochodne aromatycznych zasad Mannicha posiadają WWW, które 
zostały sklasyfikowane jako silne, bądź też średniej mocy [13-17] i wykazujące dużą 
anharmoniczność. Warto zauważyć, że — w odróżnieniu od zasad Schiffa — pojedyn- 
cze wiązanie C-N ogranicza wpływ pierścienia aromatycznego na atom akceptora. 
Na Rysunku 1 zostały przedstawione struktury molekularne pochodnych aroma- 
tycznych zasad Mannicha, dla których wykonano badania teoretyczne z zastosowa- 
niem metod dynamiki ab initio. 
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Rysunek 1. Struktur molekularne omawianych pochodnych aromatycznych zasad Mannicha. Linią przery- 
waną zaznaczono obecność wewnątrzcząsteczkowego wiązania wodorowego 

Figure 1. Molecular structures of discussed aromatic derivatives of Mannich bases. Dotted line indicates the 
presence of an intramolecular hydrogen bond 


Pochodna zasady Mannicha, 3,5,6-trimetylo-2(N,N-dimetylo-aminometylo) 
fenol (1) posiada WWW średniej mocy, które jest jednak silniejsze porównując 
z innymi zasadami Mannicha. Efekt ten przypisuje się efektom sterycznym wpro- 
wadzonym przez obecność grup CH, w bliskim sąsiedztwie atomu tlenu (donora) 
[13]. Na Rysunku 2 została przedstawiona w sposób graficzny ewolucja parametrów 
metrycznych mostka wodorowego w funkcji czasu [6]. Symulacje CPMD były pro- 
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wadzone w fazie gazowej (pojedyncza molekuła), jak również w fazie krystalicznej 
(uwzględnione więc zostały efekty towarzyszące obecności otoczenia, a zwłaszcza 
sąsiadujących molekuł w komórce krystalicznej). 
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Rysunek 2. Ewolucja w czasie parametrów metrycznych mostka wodorowego w 3,5,6-trimetylo-2-(N, N-di- 
metyloaminometylo)-fenolu (1). Lewy panel — rezultaty symulacji w fazie gazowej; prawy panel 
- rezultaty symulacji w fazie krystalicznej 

Figure 2. Time-evolution of metric parameters of atoms involved in the hydrogen bond formation in 
3,5,6-trimethyl-2-(N,N-dimethylaminomethyl)-phenol (1). Left panel — results of the gas phase 
simulations; right panel — results of the crystalline phase simulations 


W czasie trwania symulacji nie obserwujemy zjawiska przeniesienia protonu. 
Rezultaty symulacji wykonanych w fazie gazowej wskazują na małe zmiany dłu- 
gości wiązania O-H w czasie trwania symulacji. Natomiast odległości miedzyato- 
mowe O--N i H--N charakteryzowały się znacznymi fluktuacjami. Średnia wartość 
odległości O---N wyniosła 2,684 A. Analiza trajektorii otrzymanych dla symulacji 
w fazie krystalicznej wykazała większą mobilność protonu w mostku wodorowym. 
Odległości O-H i H---N niejednokrotnie miały zbliżone wartości (atom wodoru był 
współdzielony przez atomy donora i akceptora, czyli był zlokalizowany pośrodku 
mostka wodorowego). Zaobserwowano (prawy panel na Rysunku 2) również, że 
odległość pomiędzy atomami O---N nie zmienia się znacząco, natomiast odległość 
HN zmienia się znacznie silniej, niż w czasie symulacji w fazie gazowej. Obliczona 
średnia odległość O---N wynosi 2,600 A i nie różni się znacząco od danych 
uzyskanych eksperymentalnie tj. 2,628 A [13]. Przyjrzyjmy się teraz spektralnej ana- 
lizie eksperymentalnej [13] i uzyskanej z symulacji CPMD. Wartości drgania O-H 
zasady Mannicha (1) eksperymentalne i obliczone zostały zestawione w Tabeli 1. 
Wartości eksperymentalne obejmują wyniki otrzymane w CCL, jak również w fazie 
krystalicznej. Rezultaty symulacji otrzymanych dla fazy gazowej zostały porównane 
z częstościami drgania O-H otrzymanymi w CCL,, gdyż jest to rozpuszczalnik o małej 
polarności (cząsteczka CCl, z uwagi na symetrię nie ma momentu dipolowego, lecz 
wyższe momenty elektryczne) i możemy go potraktować, jako realistyczne przybli- 
żenie rezultatów pomiarów widm w fazie gazowej. 
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Tabela 1. Zestawienie eksperymentalnych i obliczonych (CPMD) wartości drgania O-H (w cm”) 
w 3,5,6-trimetylo-2-(N,N-dimetyloaminometylo)-fenolu (1) 
Table 1. Comparison of experimental and computed (CPMD) values of the O-H stretching (in cm”) 


in 3,5,6-trimethyl-2-(N,N-dimethylaminomethyl)-phenol 


Faza gazowa Faza krystaliczna 
Drganie O-H wartości eksperymentalne 
2200 + 2400 cm" ca 1700 + 3100 cm" 
CPMD (wartości z widm mocy z prędkości atomowych) | 2350 + 3250 cm 1900 + 3000 cm™ 
CPMD (wartości z widm momentów dipolowych) 2200 + 3250 cm™ 1800 ~ 2900 cm" 


* Dane eksperymentalne wzięte z Ref. [13] 


Parametry techniczne symulacji CPMD zostały dobrane w taki sposób, aby 
możliwe było otrzymanie trajektorii z momentów dipolowych, a także z prędkości 
atomowych. Na Rysunku 3 zostały przedstawione widma otrzymane jako rezultat 
analiz trajektorii momentów dipolowych. Obliczone widma z trajektorii momen- 
tów dipolowych pozwalają na poprawny opis intensywności drgań i liczb falowych, 
natomiast niemożliwa jest dekompozycja widma na poszczególne składowe. 
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Rysunek 3. Widma IR otrzymane z trajektorii momentów dipolowych dla 3,5,6-trimetylo-2-(N,N-dimetylo- 
aminometylo)-fenolu (1). Lewy panel - wyniki symulacji w fazie gazowej; prawy panel — wyniki 
symulacji w fazie krystalicznej. Doświadczalne widmo IR - wstawka w prawym panelu 

Figure 3. Predicted IR spectra obtained from the dipole moments trajectories of 3,5,6-trimethyl-2-(N,N- 
dimethylaminomethyl)-phenol (1). Left panel - gas phase simulations result; right panel — crystal- 
line phase simulations result. The experimental IR spectrum is presented in the inset 


Jak widać na Rysunku 3 (lewy panel), widmo otrzymane dla (1) jako rezul- 
tat symulacji w fazie gazowej, charakteryzuje się silnymi intensywnościami drga- 
nia O-H. Zakres drgań O-H obejmuje drgania od 2200 cm™ do 3250 cm”. Eks- 
perymentalny zakres drgania O-H obejmuje przedział od 2200 cm™ do 3400 em" 
[13]. A zatem otrzymany model obliczeniowy wykazuje dobrą ilościową zgodność 
z danymi doświadczalnymi. W prawym panelu Rysunku 3 zostały zamieszczone 
rezultaty symulacji i pomiarów IR w fazie krystalicznej [6]. Jak widzimy intensyw- 
ności drgania O-H również są wyraźnie widoczne na obu widmach. W przypadku 
widma obliczonego, widoczny jest zakres drgań od 1800-2900 cm™ z maksimum 
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przy częstości 2450 cm '. Drugim rodzajem analizy spektralnej w oparciu o symu- 
lacje CPMD, jest otrzymanie widm IR z prędkości atomowych (widm mocy). Na 
Rysunku 4 zostały przedstawione widma mocy otrzymane dla (1), jako rezultat 
symulacji w fazie gazowej i krystalicznej. Widma te odtwarzają położenie pasm 
widma oscylacyjnego, ale nie intensywności drgań. Zaletą tego typu analiz jest moż- 
liwość dekompozycji otrzymanych obliczonych widm. Jak widzimy na Rysunku 4 
znajduje się widmo otrzymane dla cząsteczki 3,5,6-trimetylo-2-(N,N-dimetyloami- 
nometylo)-fenolu (1), a także widmo otrzymane tylko dla ruchów protonu uczest- 
niczącego w WWW to ostatnie jest zdominowane przez wkład rozciągający O-H, 
ale wykazuje cechy modulowania przez ruchy atomów ciężkich. 
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Rysunek 4. Widma mocy z prędkości atomowych 3,5,6-trimetylo-2-(N,N-dimetyloaminometylo)-fenolu 
(1). Lewy panel — rezultat symulacji w fazie gazowej; prawy panel - rezultat symulacji w fazie 
krystalicznej. Kolor ciemnoszary — wszystkie atomy; kolor jasnoszary - wkład protonu mostka 
wodorowego 

Figure 4. Atomic velocity power spectra of 3,5,6-trimethyl-2-(N,N-dimethylaminomethyl)-phenol (1). Left 
panel — gas phase simulations result; right panel — crystalline phase simulations result. Dark grey 
— all atoms, Light grey — bridged proton contribution 


Na widmie mocy otrzymanym jako rezultat symulacji w fazie gazowej (Rys. 4, 
lewy panel), wyróżnić możemy trzy obszary drgań „protonu mostkowego”: 
850-1000 cm”, 1200-1600 cm” i 2350-3250 cm™. Ostatni zakres koresponduje 
z zakresem drgań grupy O-H i maksimum intensywności drgań jest w zakresie 
2650-2750 cm”. Na widmie mocy otrzymanym jako rezultat symulacji w fazie 
krystalicznej (Rys. 4, prawy panel) możemy wyróżnić dwa obszary absorpcji drga- 
nia „protonu mostkowego”: 800-1600 em `, 1900-3000 cm”. Rezultaty symulacji 
CPMD (w dwóch stanach skupienia) w postaci widm mocy z prędkości atomowych 
i momentów dipolowych pokazały, że zakres drgań O-H jest zależny od otoczenia 
chemicznego. Warto przypomnieć, że otrzymany opis spektralny jest opisem kla- 
sycznym — dynamika CPMD traktuje bowiem jądra atomowe jako cząstki klasyczne. 

Pochodna zasady Mannicha, 4,5-dimetylo-2-(N,N-dimetyloaminometylo)- 
-fenol (2), od związku (1) różni brakiem podstawienia pierścienia aromatycznego 
w pozycji orto względem grupy hydroksylowej. Symulacje CPMD dla tego związku 
zostały wykonane w trzech stanach skupienia: gazowym, ciekłym i stałym [18, 19]. 
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WWW obecne w związku (2) zostało sklasyfikowane jako średniej mocy [13]. Nie 
zaobserwowano zjawiska przeniesienia protonu w czasie trwania symulacji. Na 
Rysunku 5 zostały przedstawione rezultaty analizy zmian parametrów metrycznych 
mostka wodorowego w funkcji czasu. W fazie gazowej i roztworze proton w czasie 
trwania symulacji znajduje się po stronie atomu donora. Długość wiązania O-H nie 
ulega znacznym zmianom w czasie trwania symulacji. Widzimy jednak, że odległo- 
ści międzyatomowe pomiędzy donorem (atom tlenu) a akceptorem (atom azotu), 
a także w wiązaniu wodorowym HN. zmieniają się znacznie mocniej. 


3.5 r T T T T 3.5 


h hum? Ww] AMI WY ik My Ay byl 


25 | - asf} 


IW 
ro 

H 
D 


TT E TN d + Fav egal LAng taa oaa anoa a: 


05 1 1 1 1 1 1 05 


4 2 3 
Czas / Time [ps] Czas / Time [ps] 
35 


3 


n A i 


Czas / Time [ps] 


Rysunek 5. Ewolucja w czasie parametrów metrycznych mostka wodorowego w 4,5-dimetylo-2-(N,N-dime- 
tyloaminometylo)-fenolu (2). Lewy góry panel — rezultaty symulacji w fazie gazowej; prawy górny 
panel — rezultaty symulacji w roztworze CCl; środkowy dolny panel — rezultaty symulacji w fazie 
krystalicznej 

Figure 5. Time-evolution of metric parameters of atoms involved in the hydrogen bond formation in 4,5-di- 
methyl-2-(N,N-dimethylaminomethyl)-phenol (2). Left upper panel — results of the gas phase si- 
mulations; right upper panel — results of the liquid-phase simulations; middle lower panel - results 
of the crystalline phase simulations 


W fazie krystalicznej zaobserwowano, że proton w czasie trwania symulacji 
znajdował się w środku mostka wodorowego (tzw. „proton sharing event”). Inna 
dynamika mostka wodorowego związana jest z obecnością molekuł w komórce 
krystalicznej (komórka krystaliczna zawiera 4 molekuły), a także efektami elektro- 
statycznymi [18]. Duża mobilność protonu w mostku wodorowym, związana jest ze 
skróceniem odległości O---N. Widzimy również, że potencjał O-H silnie zależy od 
odległości O---N. Otrzymane średnie odległości międzyatomowe wynoszą 2,626 A 
dla O--N, 1,048 A dla O-H i 1,662 A dla H--N. Obliczenia w fazie krystalicznej, 


MOBILNOŚĆ PROTONU W MOSTKACH WODOROWYCH 189 


zwłaszcza odtworzenie odległości O---N, która różni się od wartości zmierzonej rent- 
genograficznie tylko o 0,02 A (wartość eksperymentalna to: 2,645 A [13]), pokazuje, 
że symulacje CPMD są w stanie odtworzyć poprawnie dynamikę mostka wodoro- 
wego, w układach posiadających wewnątrzcząsteczkowe wiązania wodorowe. Na 
Rysunku 6 został przedstawiony model zastosowany do obliczeń w roztworze. Jest 
to tzw. dyskretny model rozpuszczalnikowy. W sąsiedztwie związku (2) widoczne 
są cząsteczki rozpuszczalnika - CCl,. Szczegóły techniczne symulacji znalazły się 
w Ref. [18, 19]. Rezultaty CPMD otrzymane dla związku (2) zostały również czę- 
ściowo przedstawione w Ref. [20]. 


Rysunek 6. Model do obliczeń CPMD w roztworze - pojedyncza molekuła 4,5-dimetylo-2-(N,N-dimetylo- 
aminometylo)-fenolu (2) w niepolarnym rozpuszczalniku (CCI,) 

Figure 6. CPMD simulations model in solution - single molecule of 4,5-dimethyl-2-(N,N-dimethylamino- 
methyl)-phenol (2) in nonpolar solvent (CCL) 


Analiza spektralna (symulacje widm IR) została wykonana w oparciu o tzw. 
widma mocy z prędkości atomowych. Rezultaty wykonanych analiz zostały przed- 
stawione na Rysunku 7. 
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Rysunek 7. Widma mocy z prędkości atomowych 4,5-dimetylo-2-(N,N-dimetyloaminometylo)-fenolu (2). 
Lewy góry panel — wyniki symulacji w fazie gazowej; prawy górny panel - wyniki symulacji w roz- 
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Figure 7. Atomic velocity power spectra of 4,5-dimethyl-2-(N,N-dimethylaminomethyl)-phenol (2). Left 


upper panel — results of the gas phase simulations; right upper panel — results of the simulations 
in solution; middle lower panel - results of the crystalline phase simulations. The experimental IR 


spectrum is presented in the inset 


Otrzymane widma wskazuja, ze zakres drgan grupy O-H jest zalezny od oto- 
czenia chemicznego. Rezultaty symulacji w fazie gazowej wskazaly istnienie szero- 
kiego pasma drgań O-H w zakresie 2500-3400 cm '. W roztworze obserwujemy 
poszerzenie pasma drgań O-H. Zakres absorpcji wynosi 2300-3400 cm”. Najin- 
tensywniejszą absorpcję obserwujemy w zakresie drgań od 2600 cm™ do 2800 cm”. 
W fazie krystalicznej zakres drgań grupy O-H obejmuje obszar absorpcji od 
1700 cm™ do 2900 cm™. Zmierzony eksperymentalnie obszar drgań grupy O-H 
obejmuje zakres 2200-3400 om [13]. W celu dokładniejszego opisu zakresu drgań 
grupy O-H dokonano kwantyzacji w oparciu o metodę ,,snapshot-envelope” [12]. 
Na Rysunku 8 przedstawiono rezultaty wykonanych symulacji w różnych stanach 
skupienia, a także dla protonu i deuteru. 
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Rysunek 8.  Obwiednia pasma drgania rozciągające O-H mostka wodorowego według metody „snapshot- 
-envelope” w 4,5-dimetylo-2-(N,N-dimetyloaminometylo)-fenolu (2). Lewy góry panel - rezul- 
taty symulacji w fazie gazowej; prawy górny panel — rezultaty symulacji w roztworze; środkowy 
dolny panel - rezultaty symulacji w fazie krystalicznej 

Figure 8. Simulated hydrogen bridge O-H stretching envelope of 4,5-dimethyl-2-(N,N-dimethylamino- 
methyl)-phenol (2). Left upper panel — results of the gas phase simulations; right upper panel 
— results of the liquid-phase simulations; middle lower panel - results of the crystalline phase 
simulations 


Co zatem uwzględnienie efektów kwantowych zmieniło w opisie spektralnym 
omawianej zasady Mannicha ? Rezultaty otrzymane dla fazy gazowej wskazały obszar 
absorpcji O-H od 1900-3200 cm™. Pasmo O-H przesunięte jest nie tylko w stronę 
wyższych częstości, ale jest ono węższe [18]. W roztworze CCl, widoczny jest wpływ 
rozpuszczalnika na szerokość pasma. Jest ono szersze i obejmuje (porównując z fazą 
gazową) zakres drgań od 1600 cm™ do 3400 cm” [19]. W fazie krystalicznej zostało 
otrzymane szerokie pasmo drgań O-H w zakresie 800-2900 cm `. Wykonane symu- 
lacje CPMD wykazały, że najsilniej na właściwości fizyko-chemiczne omawianej 
zasady Mannicha wpływa obecność środowiska polarnego, natomiast wpływ sił 
dyspersyjnych wydaje się być znikomy. 

Pochodna zasady Mannicha, dla której chcemy omówić rezultaty symulacji 
CPMD i PIMD - o-dimetyloaminometylofenol (3), posiada najprostszą budowę 
chemiczną w porównaniu do pozostałych molekuł omawianych w niniejszej pracy. 
Pierścień fenylowy nie posiada podstawników, a zatem efekty steryczne mające 
wpływ na siłę wiązania wodorowego pochodzą wyłącznie od grupy aminowej. 
Symulacje CPMD i PIMD były prowadzone w fazie gazowej. Na Rysunku 9 w spo- 
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sob graficzny zilustrowano rezultaty analiz trajektorii CPMD [21] Rezultaty analizy 
parametrów metrycznych odległości miedzyatomowych w mostku wodorowym 
(lewy panel, Rys. 9) wykazały, że w czasie trwania symulacji proton zlokalizowany 
jest po stronie atomu donora. Nie zaobserwowano znaczących zmian w długości 
wiązania O-H (średnia długość tego wiązania jest równa 1,018 A). Większą zmien- 
ność odległości międzyatomowych w czasie trwania 12 ps symulacji odnotowano 
dla odległości donor-akceptor (średnia wartość tej odległości wynosi 2,769 A) 
i H--N (średnia wartość jest równa 1,857 A). Analiza spektralna została wykonana 
w oparciu o widma mocy (prawy panel, Rys. 9). 
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Rysunek 9. Lewy panel - ewolucja w czasie parametrów metrycznych mostka wodorowego w o-dimetylo- 
aminometylofenolu (3). Prawy panel - widma mocy z prędkości atomowych. Rezultaty symulacji 


CPMD w fazie gazowej 

Figure 9. Left panel - time-evolution of metric parameters of hydrogen bridge in o-dimethylaminome- 
thylphenol (3). Right panel - atomic velocity power spectra. CPMD results for the gas phase 
simulations 


Obliczony zakres absorpcji grupy O-H znajduje się w przedziale od 2550 om" 
do 3300 cm i koresponduje z danymi eksperymentalnymi. Z pomiarów IR wynika, 
że drgania grupy O-H wskazują na obszar od 2600 cm™ do 3450 cm™ [14]. Obli- 
czone widma CPMD są przesunięte w stronę niższych częstości, co ma związek ze 
specyfiką metody CPMD [8]. W celu dokładniejszego opisu WWW w omawianej 
zasadzie Mannicha, zostały wykonane symulacje w oparciu o dynamikę kwantową 
— metoda PIMD. 
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Rysunek 10. Porównanie położenia protonu w mostku wodorowym w fazie gazowej dla związku (3). Lewy 
panel — rezultaty CPMD, prawy panel -rezultaty PIMD 

Figure 10. Comparison of the proton position in the hydrogen bridge of (3) in the gas phase. Left panel — 
CPMD results, right panel - PIMD results 


Otrzymane rezultaty zostały przedstawione w postaci 2D histogramów 
(Rys. 10). Jak widzimy, skwantowanie ruchów jąder atomowych w czasie symu- 
lacji wpłynęło na opis parametrów metrycznych mostka wodorowego. Na histo- 
gramach zostały przedstawione zależności zmiany odległości międzyatomowych 
pomiędzy atomami ciężkimi tworzącymi mostek wodorowy O*''N i długością 
wiązania O-H. Skróceniu uległa odległość międzyatomowa O---N, która w czasie 
symulacji CPMD fluktuowała pomiędzy 2,50-3,05 A. W przypadku rezultatów 
z metody PIMD, odległość ta jest w zakresie 2,50-3,00 A. Jak wynika z histo- 
gramu, odwrotną zależność zaobserwowano w przypadku wiązania O-H. Analiza 
trajektorii symulacji PIMD wskazała na zakres od 0,90 A do 1,20 A, a rezultaty 
symulacji CPMD 0,98-1,10 A. Podsumowując, wprowadzenie efektów kwanto- 
wych wskazało na większą delokalizację protonu w mostku wodorowym. 

Ostatnim z omawianych, ale bardzo interesującym związkiem, jest N-tlenek 
zasady Mannicha 2-(N-dietyloamino-N-oksymetylo)-4,6-dichlorofenol (4). Symu- 
lacje CPMD dla tego związku zostały wykonane w fazie gazowej i w fazie krysta- 
licznej [22]. Eksperymentalne badania strukturalne i spektroskopowe wykazały, że 
związek ten posiada bardzo krótkie i silne WWW [15, 17]. Odległości O---O, wyzna- 
czone neutronograficznie, wynoszą odpowiednio 2,400 A i 2,423 A [16]. Zastoso- 
wanie modelu dynamicznego CPMD pozwoliło na dokładną analizę dynamicznej 
natury protonu w mostku wodorowym. Silna delokalizacja mostkowego protonu 
posłużyła za inspirację do dalszych badań strukturalnych i spektroskopowych [22]. 


194 A. JEZIERSKA, J.J. PANEK 


kr 
* 


d 


—> 


= 


Mi 


012 


W 


Rysunek 11. Ewolucja w czasie parametrów metrycznych mostka wodorowego w 2-(N-dietyloamino-N- 
oksymetylo)-4,6-dichlorofenolu (4). Lewy panel - rezultaty symulacji w fazie gazowej; prawy 
panel — rezultaty symulacji w fazie krystalicznej 

Figure 11. | Time-evolution of metric parameters of atoms involved in the hydrogen bond formation in 
2-(N-diethylamino-N-oxymethyl)-4,6-dichlorophenol. Left panel — results of the gas phase simu- 
lations; right panel — results of the crystalline phase simulations 
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Na Rysunku 11 zostały przedstawione rezultaty ewolucji parametrów metrycz- 
nych w funkcji czasu mostka wodorowego. Średnia odległość O--O, otrzymana 
jako rezultat obliczeń CPMD w fazie gazowej wynosi 2,510 A, natomiast średnie 
odległości O-H i H---O wynoszą odpowiednio 1,095 A i 1,431 A. Obliczona średnia 
wartość odległości O---O jest znacznie większa, niż zmierzona neutronograficznie 
[16]. Analiza trajektorii CPMD wskazała, że proton zlokalizowany jest po stronie 
donora w czasie trwania symulacji, jednak w trakcie zaobserwowano nieliczne 
zdarzenia, w których proton znajdował się pośrodku mostka wodorowego, a nawet 
przenosił się na stronę akceptora (lewy panel, Rysunek 11). Rezultaty symulacji 
CPMD w fazie krystalicznej wykazały, że proton w mostku wodorowym jest silnie 
zdelokalizowany i stale porusza się pomiędzy atomami donora (atom O) i akceptora 
(atom O). Z danych eksperymentalnych [16] wynika, że proton w krótszym mostku 
wodorowym jest zlokalizowany po stronie donora. W oparciu o wyniki symula- 
cji CPMD (model dynamiczny) widzimy, że proton przesunięty w stronę atomu 
akceptora. Otrzymane z obliczeń CPMD w fazie krystalicznej, średnie odległości 
miedzyatomowe O---O wynoszą odpowiednio 2,446 A i 2,451 A [22]. Na Rysunku 12 
zostały przedstawione rezultaty analizy spektralnej -obliczone widma mocy z pręd- 
kości atomowych. 
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Rysunek 12. Widma mocy z prędkości atomowych 2-(N-dietyloamino-N-oksymetylo)-4,6-dichlorofenolu 
(4). Lewy panel — rezultaty symulacji w fazie gazowej; prawy panel -rezultaty symulacji w fazie 
krystalicznej 

Figure 12. Atomic velocity power spectra of 2-(N-diethylamino-N-oxymethyl)-4,6-dichlorophenol (4). Left 
panel — results of the gas phase simulations; right panel -results of the crystalline phase simulations 


Na widmie mocy, otrzymanym z symulacji w fazie gazowej (lewy panel, Rys. 12), 
widoczne jest szerokie pasmo drgań O-H---O w zakresie od 1500 do 2700 cm". Sze- 
roki zakres absorpcji drgania O-H wskazuje, że jest ono silnie sprzężone z pozo- 
stałymi drganiami wewnątrz molekuły. Rezultaty analiz spektralnych w fazie kry- 
stalicznej (prawy panel, Rys. 12), również pokazały szeroki zakres absorpcji drgań 
pasma O-H w zakresie 500-1800 cm", co jest zgodne z danymi eksperymentalnymi 
[17]. Wyniki obliczeń CPMD wskazują na silną delokalizację protonu w mostku 
wodorowym, zwłaszcza w fazie krystalicznej, co jest w zgodzie ze spektroskopowymi 
danymi eksperymentalnymi, które wskazują na umiejscowienie protonu w środku 
mostka wodorowego (tzw. „proton shared hydrogen bond”). Rozważania o dynamice 
mostka wodorowego w omawianym N-tlenku zasady Mannicha, zostały również 
częściowo przedstawione w Ref. [20]. 


2. POTENCJAŁ ŚREDNIEJ SIŁY (PMF) 


Przedstawione w poprzednim rozdziale przebiegi odległości związanych 
z WWW podkreślają dynamiczny charakter tego oddziaływania. Parametry ener- 
getyczne — np. chwilowy potencjał odczuwany przez proton mostkowy — są silnie 
modulowane przez otoczenie, zwłaszcza przez odległość donor-akceptor. Zmiany 
w obrębie układu aromatycznego mają mniejszy chwilowy wpływ na potencjał 
ruchu protonu, ponieważ pojedyncze wiązanie C-N szkieletu zasady Mannicha 
ogranicza efekt rezonansowy (RAHB). Jednak podstawniki modyfikujące inne poza 
WWW fragmenty cząsteczki, w tym układ aromatyczny, także zmieniają rowno- 
wagową strukturę molekularną, a zatem w sposób pośredni modyfikują całościowe 
zachowanie się WWW. Nie jest możliwe oddanie tych efektów za pomocą pojedyn- 
czej krzywej energii potencjalnej dla ruchu protonu przy równowagowej struktu- 
rze cząsteczki. Jednak zastosowanie schematów dynamiki molekularnej, w naszym 
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przypadku CPMD, umożliwia odtworzenie przebiegu energii swobodnej wzdłuż 
wybranej współrzędnej. Taka zależność nosi nazwę potencjału średniej siły (pmf), 
co ma podkreślać statystyczną naturę tej wielkości, pochodzącej z próbkowania 
przestrzeni fazowej. Dzięki temu uwzględniony jest wpływ różnych konfiguracji 
otoczenia, zwłaszcza — dla badań w roztworze — zmiennej, dynamicznej otoczki sol- 
watacyjnej. 

W pracy [18] wykorzystaliśmy metodę ,,snapshot-envelope” nie tylko do 
odtworzenia obwiedni pasm rozciągających WWW O-H (O-D) dla 4,5-dimetylo- 
-2-(N,N-dimetyloaminometylo)-fenolu (2), ale też do wyznaczenia rozkładu gęsto- 
ści prawdopodobieństwa położenia protonu / deuteronu. Uśrednienie serii takich 
rozkładów dało możliwość wyznaczenia przebiegu 1D pmf dla (2) i są one przed- 
stawione na Rysunku 13. Jak widać, WWW w związku (2) ma wyraźnie zaznaczoną 
anharmoniczność, szczególnie w fazie stałej, dla której widoczne jest też nieznaczne 
przesunięcie położenia minimum pmf w kierunku wyższych długości O-H (O-D). 
Zwiększona masa deuteronu, a zatem bardziej skupiona postać oscylacyjnych funk- 
cji falowych, prowadzi do widocznej różnicy między przebiegiem pmf dla protonu 
i deuteronu. Jednak, zgodnie z wcześniejszym opisem parametrów strukturalnych, 
WWW w (2) nie wykazuje skłonności do trwałego przeniesienia protonu na stronę 
akceptora. 
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Rysunek 13. Przebieg potencjału średniej siły (pmf) dla ruchu protonu / deuteronu w mostku wodorowym 
związku (2). 
Figure 13. ` Potential of mean force for proton / deuteron motion in the hydrogen bridge of compound (2) 


Ciekawym rozwinięciem metodologicznym było, dokonane przez nas w pracy 
[21], rozszerzenie tego jednowymiarowego podejścia na przypadek dwuwymia- 
rowy. Jako zmienne zbiorowe do opisu mostka wodorowego prostej zasady Man- 
nicha, o-dimetyloaminometylofenolu (3), użyto długości O-H oraz kąta O-H---N. 
Dzięki temu, jak pokazuje Rysunek 14, uwidoczniono znaczną kątową labilność 
protonu. Nawet jeśli zakres swobodnych ruchów protonu obejmuje głównie odle- 
głości donor-proton do 1,2-1,3 A, to za to w obszarze kątów jest to od 120° do 170°. 
Szkielet kowalencyjny cząsteczki uniemożliwia preferencję dla położenia liniowego, 
ale widoczne jest lekkie rozszerzanie się obszaru swobodnego ruchu protonu ku 
górze rysunku (ku większym kątom O-H---N). 
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Rysunek 14. Mapa 2D pmf dla ruchu protonu w WWW o-dimetyloaminometylofenolu (3) 
Figure 14. | A 2D pmf map for the proton motion in the hydrogen bridge of o-dimethylaminemethylphenol (3) 


Analogiczne mapy 2D pmf uzyskaliśmy dla N-tlenka zasady Mannicha, 
2-(N-dietyloamino-N-oksymetylo)-4,6-dichlorofenolu (4). W tym przypadku, jak 
pokazano w rozdziale 2, jest to wiązanie bardzo labilne, krótkie, z częstymi przesko- 
kami protonu na stronę atomu akceptora. Jest to dobrze widoczne na mapie pmf, 
pokazanej na Rysunku 15. Przy długości mostka O---O wynoszącej ok. 2,5 A, odleg- 
łości donor-proton sięgające 1,3 A to już przeniesienie protonu na stronę atomu 
akceptorowego. Także pozycja minimum (ok. 1,15 A) jest znacznie przesunięta ku 
środkowi mostka. Te dane oznaczają, że WWW ma charakter studni potencjału 
z pojedynczym minimum zlokalizowanym blisko środka mostka. Nie występuje 
przeniesienie protonu od minimum przy atomie donora do drugiego minimum 
przy atomie akceptora. WWW związku (4) jest też znacznie bardziej liniowe niż 
dla opisanej powyżej zasady Mannicha (3), a proton często znajduje się wewnątrz 
quasi-pierścienia (domkniętego przez linię atom donora — atom akceptora); takie 


przypadki odpowiadają na Rysunku 15 umownym wartościom kąta O-H---O powy- 
żej 180°. 
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Rysunek 15. Mapa 2D pmfdla ruchu protonu w WWW 2-(N-dietyloamino-N-oksymetylo)-4,6-dichlorofenolu 
(4) - wyniki symulacji CPMD w fazie gazowej 

Figure 15. A 2D pmfmap for the proton motion in the hydrogen bridge of2-(N-diethylamino-N-oxymethyl)- 
4,6-dichlorophenol (4) — results of CPMD run in the gas phase 
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PODSUMOWANIE 


Przedstawione wyniki obliczeń metodami dynamiki molekularnej ab initio 
pokazują, że dynamika WWW pochodnych zasad Mannicha jest bardzo zróżnico- 
wana. W o-hydroksypochodnych (1), (2) i (3) przeważnie obserwuje się pozostawa- 
nie protonu przy atomie donora, lecz z silnie anharmonicznym potencjałem umożli- 
wiającym przenikanie w obszar atomu akceptora. Wyraźnie widoczny jest też wpływ 
środowiska, zwłaszcza krystalicznego. Natomiast dla N-tlenka zasady Mannicha (4) 
WWW jest opisane potencjałem silnie przesuniętym ku centrum mostka i z bardzo 
szerokim minimum, zwłaszcza w fazie stałej. Dzięki temu w fazie gazowej częste są 
przeskoki protonu ku atomowi akceptora, zaś w krysztale równowagowe położenie 
protonu znajduje się pośrodku mostka (uwspólnianie protonu). Te wyniki, zgodne 
z danymi doświadczalnymi, pokazują potencjał dynamiki CPMD, schematu PIMD 
oraz analiz trajektorii dynamiki, przy badaniu mikroskopowej natury właściwości 
fizyko-chemicznych pochodnych zasad Mannicha. 
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ABSTRACT 


In the work, after preliminary discussion of the complexity of the phenomenon 
of melting chemical compounds, two basic research methods are presented: calo- 
rimetric methods and thermooptical method. The physical basis of the five main 
calorimetry techniques is now detailed (adiabatic calorimetry, differential thermal 
analysis - DTA, differential scanning calorimetry with heat compensation - DSC 
heat compensated, differential scanning calorimetry with heat flow - DSC heat flux 
and differential scanning calorimetry with temperature modulation - MDSC) and 
thermo-optical techniques used in phase transitions investigations. The advantages 
and disadvantages of these methods are shown in numerous examples and the accu- 
racy attainable by the individual measuring techniques is compared. 


Keywords: adiabatic calorimetry, DSC calorimetry with heat compensation, DSC 
calorimetry with heat flow, MDSC modulation calorimetry, differential thermal 
analysis DTA, liquid crystals, mesogens, rotational phases 

Słowa kluczowe: kalorymetria adiabatyczna, kalorymetria DSC z kompensacją 
ciepła, kalorymetria DSC z przepływem ciepła, modulacyjna kalorymetria MDSC, 
różnicowa analiza termiczna DTA, ciekłe kryształy, mezogeny, fazy rotacyjne 
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WPROWADZENIE - TOPNIENIE JAKO NAJBARDZIEJ 
ZŁOŻONY RODZAJ PRZEMIANY FAZOWEJ 


Wszystkie opisane w literaturze związki chemiczne występują w trzech stanach 
skupienia: w fazie stałej, fazie ciekłej i fazie gazowej [1]. Jeśli obserwujemy brak 
którejś z tych faz, to jest to rzadki przypadek związany z brakiem trwałości tego 
związku w zakresie temperatur, w którym spodziewamy się występowania tej fazy. 

Przemiany fazowe z termodynamicznego punktu widzenia najczęściej opisu- 
jemy korzystając z teorii Ehrenfesta [2]. Zgodnie z tą teorią w pierwszym przybli- 
żeniu możemy wszystkie przemiany fazowe podzielić na dwie grupy: przemiany 
pierwszego (I) rodzaju oraz drugiego (II) rodzaju. Przewidywanych przez tę teorię 
przemian fazowych wyższych rzędów nie zaobserwowano, a nawet pokazano w lite- 
raturze, że nie mogą występować [3]. Natomiast w najprostszym energetycznym 
opisie problem przemian fazowych pierwszego rodzaju jest bardziej złożony. Przez 
przemiany te rozumiemy takie przejścia fazowe, którym towarzyszy niezerowe 
ciepło przemiany fazowej oraz niezerowa zmiana objętości. Parametrem, który ilo- 
ściowo różnicuje te przejścia fazowe najczęściej jest wielkość towarzyszącej ental- 
pii, nazywanej od czasów Blacka ciepłem utajnionym (ang. laten heat) [4]. Spośród 
obserwowanych procesów największą wartością ciepła przemiany fazowej charakte- 
ryzuje się proces parowania. Niemniej ze względu na dużą jednorodność tego pro- 
cesu oraz fakt, że ilościowo proces ten możemy opisać teorią Troutona [1], badania 
procesów parowania pozostały już tylko domeną inżynierów. Natomiast najbardziej 
interesującym problemem badawczym i wciąż nie w pełni zrozumianym okazuje się 
proces topnienia. Jest to najbardziej energetyczna przemiana fazowa, nie wliczając 
w to oczywiście procesu parowania. 

Topnienie jest procesem niezwykle złożonym. W jego wyniku wszystkie mole- 
kuły fazy stałej zyskują 3 stopnie swobody translacji oraz 3 stopnie swobody rotacji 
(bardziej poprawnie nazwanych reorientacjami [5]). Oczywiście tak jest w przy- 
padku klasycznym. Problemem otwartym pozostaje fakt jak rozumiemy fazę stałą 
oraz czy wszystkie stopione ciecze są jednakowe pod względem dynamiki molekuł 
je tworzących. Oba zagadnienia są złożone. 

Przez ciało stałe będziemy rozumieć fazę krystaliczną o ścisłe określonej struk- 
turze najczęściej opisywanej za pomocą symetrii grupy przestrzennej oraz wymia- 
rów komórki elementarnej. W ramach fazy stałej możliwy jest polimorfizm, czyli 
możliwość tworzenia faz krystalicznych o różnych strukturach co objawia się naj- 
częściej zmianą właściwości fizycznych i najczęściej zapisujemy ten fakt jako zmianę 
grupy przestrzennej. Takim przemianom fazowym często towarzyszy ciepło prze- 
miany fazowej. Te przemiany fazowe nazywamy również przemianami typu porzą- 
dek-nieporządek. Różnią się więc one stopniem uporządkowania molekuł w sieci 
krystalicznej. Właściwsze wydaje się nazywanie tych pośrednich faz krystalicznych 
mezofazami zgodnie z nomenklaturą Friedela [6] wprowadzoną prawie sto lat temu 
w 1922 roku. Przemiany takie charakteryzują się temperaturą przemiany fazowej, 
ciepłem przemiany fazowej oraz często zmianą właściwości optycznych (np. krysz- 
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tału optycznie dwuosiowego na jednoosiowy czy zmianą kierunku osi optycznych). 
Trudnym problemem jest kinetyka tego procesu. Bardzo często takie przemiany 
pomiędzy dwoma formami krystalicznymi, zwykle nazywanymi polimorfami, nie 
zachodzą do końca, co znacząco obniża dokładność wyznaczania ciepeł przemian 
fazowych. Poza tym możliwa jest histereza temperatury, która jest źródłem niedo- 
kładności w opisie temperatury przemian fazowych podczas cyklu schładzania. 
Oczywiście właściwe, termodynamiczne, parametry przemian fazowych otrzymu- 
jemy jedynie w procesie ogrzewania podczas topnienia, lecz nie wszystkie mezofazy 
są tam widoczne. Dzięki procesowi histerezy temperatury topnienia, jej wartość 
można obniżyć nawet o 20 K. W ten sposób odsłania się obszar, w którym mogą 
wystąpić przejścia fazowe do dodatkowych mezofaz, widocznych tylko podczas pro- 
cesu chłodzenia. Fazy takie nazywamy monotropowymi. 

Dodatkowym problemem podczas obserwacji procesów topnienia i krzepnię- 
cia są procesy zeszklenia. Te procesy są najbardziej charakterystyczne dla cieczy 
utworzonych przez makromolekuły. Proces ten jest również obserwowany wśród 
molekuł o małej masie molowej, w szczególności o rozbudowanej strukturze. Pro- 
cesy zeszklenia mogą dotyczyć nie tylko normalnej izotropowej cieczy, ale również 
tych częściowo nieuporządkowanych faz krystalicznych. Te rodzaje amorficznej 
fazy stałej są obecnie niezwykle intensywnie badane [7-9]. 

Pojęcie fazy ciekłej również wymaga uściślenia. Faza ciekła powstaje po uru- 
chomieniu ruchów translacyjnych oraz rotacyjnych (reorientacyjnych) w molekułach 
tworzących fazę krystaliczną. Natomiast natura oddziaływań molekularnych, która 
była źródłem sił odpowiedzialnych za stabilność sieci krystalicznej, pozostawała 
w fazie ciekłej. Dlatego trudno porównywać właściwości soli stopionych (czyli tzw. 
cieczy jonowych), ciekłych metali (np. rtęci) z ciekłymi węglowodorami. Wg Ubel- 
lode [10] bardziej właściwsze niż nazwa ciecz byłoby używanie określenia stopiony 
stan materii (the molten state of matter) co oddaje podobieństwo w sensie natury 
oddziaływań do stanu krystalicznego. Mówiąc o cieczach spotykamy się również 
w literaturze z pojęciem cieczy anizotropowych. Są nimi fazy utworzone przez mole- 
kuły, które posiadają pełną (lub częściowo ograniczoną) swobodę translacji, nato- 
miast występuje pewne ograniczenie w ich reorientacji. Preferowany jest pewien 
kierunek ułożenia molekuł, zwany directorem. Faza taka odkryta przez Reinitzera 
w 1888 roku [11] została nazwana przez Lehmanna ciekłym kryształem [12] ze 
względu na płynność i jednocześnie anizotropię właściwości optycznych. Nazwę 
tę możemy utożsamiać, z pewnymi wyjątkami, z pojęciem mezofazy. Natomiast 
nomenklatura IUPAC [13] wyraźnie pokazuje występujące pomiędzy nimi różnice. 
Bardziej ogólna jest nazwa mezofaza, obejmuje ona również kryształy ze swobodą 
rotacji (reorientacji), czyli fazy rotacyjne. Mezofaza to faza pośrednia, w której mole- 
kuły posiadają przynajmniej jeden stopień uporządkowania i przynajmniej jeden 
stopień nieuporządkowania (translacyjnego lub rotacyjnego). Natomiast pojęcie 
fazy ciekłokrystalicznej wyklucza przypadek pełnego uprządkowania translacyj- 
nego, z zachowaniem nieuporządkowania rotacyjnego. W pozostałych przypadkach 
jest synonimem wyrażenia mezofaza. 
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1. METODY KALORYMETRYCZNE W BADANIACH 
PRZEMIAN FAZOWYCH 


Badania procesu topnienia jak również innych przejść fazowych przeprowadza 
się przede wszystkim metodami kalorymetrycznymi. Nie tylko przemiany fazowe I 
rodzaju są tam wyraźnie widoczne. Również metody te często umożliwiają badanie 
przejść fazowych II rodzaju dla tych przypadków, gdy w temperaturze przemiany 
fazowej dochodzi do mierzalnej zmiany pojemności cieplnej badanego materiału. 
W praktyce możemy stosować różne metody i techniki badawcze. Świętosławski 
w swoim podręczniku z roku 1928 „Chemja Fizyczna” [14, 15] wyróżnia 3 rodzaje 
kalorymetrów: izotermiczne, adiabatyczne oraz zwykłe (obecnie za Kubaschewskim 
[16] nazywane izoperibolicznymi). Bardzo szczegółowo ten problem dyskutował 
Zielenkiewicz w swoich licznych publikacjach [17, 18]. 

Z dzisiejszego punktu widzenia w badaniach kalorymetrycznych przejść fazo- 
wych za najczęściej stosowane (i być może najważniejsze) należy wyróżnić dwie 
podstawowe techniki. Metodę kalorymetrii adiabatycznej oraz metodę kalorymetrii 
skaningowej z modyfikacjami takimi jak: różnicowa analiza termiczna (DTA), róż- 
nicowa kalorymetria skaningowa (DSC) z kompensacją ciepła (heat compencated 
DSC), różnicowa kalorymetria skaningowa (DSC) z przepływem ciepła (heat flux 
DSC) oraz modulacyjna różnicowa kalorymetria skaningowa (MDSC). 


1.1. METODA KALORYMETRII ADIABATYCZNEJ 


Metoda ta obecnie uchodzi za najdokładniejszą, mierzącą w sposób bezpo- 
średni pojemności cieplne badanych materiałów, entalpie przemian fazowych oraz 
ich temperatury. Jej wadami są: bardzo długi czas pomiaru (obecnie mierzony 
w dniach i tygodniach) oraz duża ilość substancji (około 2 g) niezbędna do precyzyj- 
nego wyznaczenia parametrów termicznych. Metodę tę wynalazł Richards w roku 
1905 [19, 20], pomimo iż najczęściej metodę tę kojarzymy z niskotemperaturowymi 
badaniami pojemności cieplnej metali w laboratorium Nernsta [21]. Należy tu pod- 
kreślić, że kalorymetr adiabatyczny został zbudowany już w pierwszym okresie dzia- 
łalności profesora Świętosławskiego w Warszawie w latach 20. ubiegłego wieku [15]. 
Większość stosowanych kalorymetrów adiabatycznych są konstrukcjami własnymi, 
zaprojektowanymi na potrzeby własnych specyficznych badań. W Polsce jest wiele 
takich konstrukcji. Do bardziej znanych należą kalorymetry w Zakładzie Badania 
Materii Miękkiej IFJ PAN w Krakowie [22, 23] oraz w Oddziale Niskich Temperatur 
i Nadprzewodnictwa INTiBS PAN we Wrocławiu [24, 25]. 

Najważniejszą cechą tego urządzenia jest minimalizacja straty ciepła związana 
z jej ucieczką poprzez ścianki naczynia pomiarowego. Praktycznie realizuje się to 
poprzez wprowadzenie szeregu warstw pośrednich (ekranów), odizolowanych ciepl- 
nie, których temperatura nadąża za temperaturą komory pomiarowej. Dodatkowo 
obniża się możliwość ucieczki ciepła poprzez wprowadzenie próżni pomiędzy ekra- 
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nami. Bardzo wiele takich konstrukcji zostało opisanych w literaturze. Rysunek la 
pokazuje przykładową budowę prostego kalorymetru adiabatycznego, który jednak 
zawiera wszystkie niezbędne elementy [26]. Są nimi: pompa wysokiej próżni (A), 
stalowa rama naczynia Dewara (B), zbiornik ciekłego azotu (C), naczynie Dewara 
(D), naczynie próżniowe (E), zewnętrzny płaszcz adiabatyczny (F), wewnętrzny 
płaszcz adiabatyczny (G), naczynie pomiarowe z próbką (H), termopary (I). 
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Rysunek 1. a) Schemat budowy typowego kalorymetru adiabatycznego (opis w tekście) [26]. b) Molowa 
pojemność cieplna 4-pentyloksybenzylideno-4 -chloroaniliny uzyskana metodą kalorymetrii 
adiabatycznej [27] 

Figure 1. a) Construction diagram of a typical adiabatic calorimeter (description in the text) [26]. b) Molar 
heat capacity of 4-pentyloxybenzylidene-4 -chloroaniline obtained by adiabatic calorimetry [27] 


Oczywiście urządzenia te są bardziej skomplikowane jeśli mają pracować 
w temperaturach „helowych”. Istotą współczesnych konstrukcji jest możliwość 
automatycznego śledzenia warunków równowagi termicznej za pomocą licznych 
czujników temperatury. Należy tu zaznaczyć, że obecnie pomiarem sterują kompu- 
tery, co dodatkowo pokazuje niezwykłą pracowitość osób mierzących pojemności 
cieplne tą metodą w okresie przed automatyzacją procesu zbierania danych. Typowy 
wynik uzyskany tą techniką pokazuje Rysunek 1b. Przedstawia on temperaturową 
zależność molowej pojemności cieplnej 4-pentyloksybenzylideno-4 -chloroaniliny 
[27]. Wykres ten przedstawia przede wszystkim typowy kształt zależności tempe- 
raturowej pojemności cieplnej fazy stałej wyjaśniony ilościowo przez Debyea [28]. 
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Natomiast z dzisiejszego punktu widzenia najbardziej istotne są te fragmenty, które 
odbiegają od zachowania klasycznego, czyli nieciągłości zmian pojemności ciepl- 
nej w pobliżu temperatury przemiany fazowej. Zaprezentowane pomiary pozwalają 
również bardzo precyzyjnie wyznaczyć entalpie przemian fazowych oraz ich tempe- 
ratury. W tym eksperymencie zaobserwowano dwie mezofazy (SmB,,, SmA,), trzy 
fazy krystaliczne (KI, KII i KIII) oraz oczywiście fazę izotropową. Identyfikację tych 
faz przeprowadzono na drodze rentgenografii (XRD) oraz optycznej mikroskopii 
polaryzacyjnej (POM) co zostało precyzyjnie udokumentowane w cytowanej pracy [27]. 


1.2. RÓŻNICOWA ANALIZA TERMICZNA (DTA) 


Metoda ta jest modyfikacją analizy termicznej opisanej po raz pierwszy 
w literaturze przez Le Chateliera [29] i polegającej na pomiarze zmian tempera- 
tury w czasie dla próbki podgrzanej wstępnie do wysokiej temperatury, a następ- 
nie schładzanej. Ta technika była powszechnie stosowana do badania diagramów 
fazowych. Jej rozwinięciem i udokładnieniem było wprowadzenie drugiego pustego 
naczyńka porównawczego i pomiar różnicy temperatur pomiędzy nimi, w warun- 
kach umieszczenia ich w sposób symetryczny w piecu o regulowanej temperaturze. 
Technikę tę nazywamy dziś różnicową analizą termiczną (DTA). Wprowadził ją do 
literatury Robert Austen w roku 1899 [30]. Typowy schemat budowy głowicy DTA 
przedstawia Rysunek 2a. 

Dużym postępem było rozseparowanie obu bliźniaczych naczyń pomiarowych 
oraz umieszczenie termopary pomiarowej pod stolikami na których umieszczano 
badane próbki [30]. Przedłużyło to czas „życia” termopar pomiarowych oraz umoż- 
liwiło właściwą ich kalibrację. Schemat budowy takiej głowicy przedstawia Rysu- 
nek 2b. Ta modyfikacja konstrukcji dała początek technice DSC z przepływem cie- 
pła (heat flux). 

Technika DTA nadaje się do pomiarów wstępnych, jak również do pomiarów 
w temperaturach bardzo wysokich. Jej głównym problemem jest szerokość powsta- 
jących „pików” w punktach anomalii termicznych, które nie zanikają zbyt szybko ze 
względu na użycie materiałów o małym przewodnictwie cieplnym, którymi najczęś- 
ciej są ceramiki. 
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Rysunek 2. a) Schemat klasycznego przyrządu DTA. Termopary umieszczone są w próbce pomiarowej oraz 
referencyjnej. b) Schemat urządzenia DTA typu Boersma [31] z termoparami umieszczonymi pod 
dnem rozseparowanych naczyń pomiarowych (wg [32]) 

Figure 2. a) Diagram of a classic DTA instrument. Thermocouples are placed in the measuring and refe- 
rence sample. b) Diagram of a Boersma type DTA [31] with thermocouples located under the 
bottom of separated measuring vessels (according to [32]) 
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Rysunek 3. Termogram DTA (górny wykres) oraz DSC z przepływem ciepła (heat flux) substancji mezogen- 
nej bis-(4 -heptyloksybenzal)-1,4-fenylodiaminy 

Figure 3. DTA thermogram (upper graph) and DSC with heat flux of mesogenic Bis (4 -heptyloxybenzal) 
-1,4-phenyldiamine 


Jakość otrzymywanych wyników klasyczną metodą DTA ilustruje Rysunek 3. 
Przedstawiono na nim dwa termogramy tej samej substancji, bis-(4 -heptyloksyben- 
zal)-1,4-fenylodiaminy. Wykres górny został otrzymany w pomiarze na klasycznym 
DTA, który w sprzężeniu z termowagą stanowi zestaw SETARAM SetSys. Drugi, 
dolny wykres, otrzymano z pomiaru na przyrządzie Mettler Toledo DSC3 metodą 
DSC z przepływem ciepła (heat-flux). Z tych wykresów wynika, że obie metody 
identyfikują wszystkie 7 przejść fazowych I rodzaju, jak również iż temperatury 
tych przejść fazowych są podobne. Na wykresie widać również podstawowe wady 
metody DTA: dużą szerokość „pików” oraz nakładanie się blisko siebie znajdujących 
się przemian fazowych co uniemożliwia oszacowanie ich entalpii. 
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1.3. ROZNICOWA KALORYMETRIA SKANINGOWA (DSC) 
Z KOMPENSACJĄ CIEPŁA (ANG. HEAT COMPENSATED) 


Metoda ta jest podobna do metody DTA, ale w przeciwienstwie do niej celem 
głównym nie jest pomiar pojawiającej się różnicy temperatur między próbką pomia- 
rową a próbką referencyjną, natomiast metoda ta wykorzystuje tę różnicę tempera- 
tur do sterowania mocą podawaną na próbkę, aby „wyzerować” mierzalny gradient 
temperatur. Czyli jest mierzony rzeczywisty efekt cieplny dostarczony do próbki, aby 
obie próbki (pomiarowa i referencyjna) miały przez cały czas tę samą temperaturę. 
W ten sposób mierzone dodatkowe ciepło jest miarą pojemności cieplnej oraz cie- 
pła przemian fazowych I rodzaju. Schemat budowy głowicy kalorymetru typu DSC 
z kompensacją ciepła przedstawia Rysunek 4a. Istotną cechą budowy jest obecność, 
w obu platynowych podstawkach, dwóch izolowanych taśm platynowych, z których 
jedna jest czujnikiem temperatury, a druga jest grzałką. Dzięki takiej konstrukcji 
gradient temperatury możemy kompensować podgrzewając każdą z podstawek, co 
umożliwia pomiar mocy cieplnej podczas grzania i chłodzenia oraz dla procesów 
egzo i endotermicznych. 

Metoda ta powstała w laboratoriach firmy Perkin Elmer i za jej twórcę uzna- 
jemy O'Neilla [33-35]. Technika ta jest prawdziwą metodą kalorymetryczną, gdyż 
mierzy rzeczywiste efekty cieplne. 


a) 


Próbka b) LT Probka 


f g Próbka 
` referencyjna pomiarowa referencyjna 


+1—| Termometr stolika 


bk 
+ — i 
O Grzałka stolika i | | Sygnał DSC 

i Heat - flux 


me ANA ne AW 


Probka 
pomiarowa 


Rysunek 4. a) Schemat budowy głowicy pomiarowej kalorymetru skaningowego DSC z kompensacją ciepła 
(heat compensated). b) Schemat budowy głowicy pomiarowej kalorymetru skaningowego DSC 
z przepływem ciepła (ang. heat-flux) 

Figure 4. a) Construction scheme of DSC scanning calorimeter with heat compensated. b) Schematic dia- 
gram of the DSC scanning calorimeter measuring head with heat flow (heat-flux) 


1.4. ROZNICOWA KALORYMETRIA SKANINGOWA (DSC) 
Z PRZEPŁYWEM CIEPŁA (ANG. HEAT FLUX) 


Metoda ta rozwinęła się w oparciu o technikę DTA, którą po pewnych udo- 
skonaleniach konstrukcyjnych oraz dzięki rozwojowi elektroniki oraz informatyki 
udało się doprowadzić do poziomu zapewniającego konkurencyjność w stosunku do 
metody DSC z kompensacją ciepła. Schemat głowicy tego urządzenia przedstawia 
Rysunek 4b. Najważniejszą zmianą, jest to iż obie platynowe podstawki umieszczone 
są na wspólnym metalowym podłożu dającym możliwość w miarę szybkiego 
przepływu ciepła, co zapewnia znaczące zmniejszenie szerokości połówkowej 
„pików” termicznych. Grubość tego podłoża nie może jednak być zbyt duża, aby 
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nie zlikwidować sygnału gradientu temperatury. Z głowicy wychodzi tylko jeden 
sygnał, który jest poddany obróbce cyfrowej zgodnie z parametrami kalibracji. 
Wykorzystuje się tu efekt zależności tego sygnału kalorymetrycznego (AT) od 
masy próbki, ciepła przemiany fazowej oraz aktualnej temperatury. Wszystkie 
nieliniowości korygowane są numerycznie. Dzięki temu wyniki tą metodą w pełni 
przypominają wyniki otrzymane techniką DSC z kompensacją ciepła. Wniosek 
ten ilustruje Rysunek 5, na którym przedstawiono wyniki pomiarów nonanolanu 
4-oktyloksy-4 -hydroksyazobenzenu (MP-8-9) obu metodami DSC. Wykres górny 
otrzymano za pomocą techniki DSC z kompensacją ciepła na przyrządzie PERKIN 
ELMER DSC7, a dolny zmierzony został techniką DSC z przepływem ciepła na 
przyrządzie METTLER TOLEDO DSC3. Praktycznie nie widać różnicy w jakości 
otrzymanych wyników. Obie metody identyfikują 6 przemian fazowych o tych 
samych parametrach (temperaturze oraz entalpii). 
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Rysunek 5. _ Termogramy nonanolanu 4-oktyloksy-4 -hydroksyazobenzenu (MP-8-9) otrzymane metodą ka- 
lorymetrii DSC z kompensacją ciepła (krzywa górna) oraz metodą kalorymetrii DSC z przepły- 
wem ciepła (krzywa dolna) 

Figure 5. Thermograms of 4-octyloxy-4 -hydroxyazobenzene nonanoade (MP-8-9) obtained by DSC calo- 
rimetry with heat compensation (upper curve) and DSC calorimetry with heat flow (bottom 
curve) 


1.5. RÓŻNICOWA MODULACYJNA KALORYMETRIA SKANINGOWA (MDSC) 


Metoda ta jest najnowszą modyfikacją aparaturową ale również software-ową 
jaka przeniknęła w obecnym wieku do klasycznej aparatury DSC. Jej źródłem jest 
tzw. kalorymetria AC, którą wprowadził do współczesnej metodologii pomiarowej 
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Kraftmakher [36]. W swoich publikacjach przegladowych [37, 38] podkresla, ze 
pomysł tych pomiarów powstał jeszcze w latach 20. w laboratorium Corbino [39], 
natomiast ze względu na ówczesny stan elektroniki rozwijanie tej metodologii nie 
było możliwe. Dopiero w latach 70. ubiegłego wieku w MIT, w oparciu o wcześniej- 
sze wyniki Kraftmakhera, skonstruowano aparaturę pozwalającą tą techniką badać 
zagadnienia związane z krytycznością przemian fazowych [40-42]. W technice 
tej można było bardzo precyzyjnie mierzyć temperaturę oraz pojemności cieplne 
w pobliżu przemian fazowych co jest warunkiem niezbędnym do opisu przejść fazo- 
wych za pomocą indeksów krytycznych. 

Istota metody jest następująca. Do bardzo cienkiej próbki (np. w postaci 
pastylki) przymocowuje się z jednej strony małą grzałkę, a z drugiej strony czuły 
termometr. Następnie do grzałki podłącza się sygnał zmienny o ściśle określonej 
częstotliwości. Eksperyment polega na pomiarze amplitudy zmian temperatury oraz 
jej fazy względem podawanego sygnału napięciowego. Z tych danych możemy pre- 
cyzyjnie wyznaczyć pojemność cieplną tej próbki. Niestety dokładność zmniejsza 
się wraz z przybliżeniem się ku temperaturze przemiany fazowej, gdyż amplituda 
zmian temperatury jest odwrotnie proporcjonalna do pojemności cieplnej. Dlatego 
niekiedy w najbliższym otoczeniu przemiany fazowej zmienia się technikę pomiaru 
na podawanie ściśle określonej porcji energii i pomiarze przyrostu temperatury, 
podobnie jak w pomiarach kalorymetrii adiabatycznej [43]. W metodzie kalory- 
metrii AC częstość sygnału na grzałce dobiera się doświadczalnie, aby straty ciepła 
były minimalne. Aparatura pomiarowa w kalorymetrii AC była wytwarzana w wielu 
laboratoriach i nie jest komercyjną. 

Technikę tę firmy komercyjne przeniosły do aparatury DSC tworząc tzw. 
metodę MDSC. Praktycznie wszystkie aparaty DSC mogłyby stosować tę technikę. 
Polega to na nałożeniu na liniowy wzrost temperatury głowicy sygnału zmiennego 
o małej amplitudzie. Wyniki uzyskiwane w tych pomiarach są bardzo mocno uza- 
leżnione od dodatkowego oprogramowania, co powoduje, że ta technika jest wciąż 
bardzo mało popularna [44]. 

Postęp w elektronice jest tak olbrzymi, iż obecnie produkowane są komercyjnie 
chipy, np. XEN-LCM2506, które są głowicą pomiarową do metody kalorymetrii AC 
znanej pod nazwą nanokalorymetrii. Dzięki temu można samemu skonstruować 
aparat do badań kalorymetrycznych w skali nano [45]. 


2. METODA OPTYCZNEJ MIKROSKOPII POLARYZACYJNEJ 
W BADANIACH PRZEMIAN FAZOWYCH 


Jak już było wspomniane powyżej przemianom fazowym towarzyszą nie tylko 
anomalie termiczne, ale również zmiana struktury krystalicznej, która szczególnie 
dobrze widoczna jest w świetle spolaryzowanym. Obrazy takich tekstur można 
spotkać już w literaturze XIX-wiecznej. Przykładem może być zestaw dla ogrom- 
nej grupy kryształów związków organicznych opublikowany przez Lehmanna 
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w 2-tomowym dziele Molekular Physik wydanym w latach 1988-1989, czyli przed 
erą badań rentgenowskich [46]. Wynika z nich, iż śledzenie zmian obrazów tek- 
stur jest również dobrą metodą badania przemian fazowych. Wciąż najlepszym 
przyrządem optycznym jest oko ludzkie i stąd najczęściej stosowane obserwacje 
procesów topnienia za pomocą powszechnie stosowanego stolika Boetiusa. Jest to 
bardzo żmudna praca wymagająca cierpliwości i nie zawsze pozwalająca precyzyj- 
nie wyznaczyć temperaturę przejścia fazowego. Dlatego gdy tylko rozwój elektroniki 
osiągnął pewien poziom rozpoczęto poszukiwanie metod zastąpienia obserwacji 
pomiarami aparaturowymi. 


2.1. ZASADY POMIARÓW METODĄ ANALIZY TERMOOPTYCZNEJ (TOA) 


Automatyzację obserwacji procesów topnienia wprowadził Barrall II w latach 
60-tych, w okresie powstawania metody DSC - z kompensacją ciepła [47-49]. Sto- 
sowana metoda była bardzo prosta i została przedstawiona na Rysunku 6a [47, 48]. 
Stolik grzewczy z regulacją temperatury oraz czujnikiem temperatury umieszczano 
w świetle mikroskopu polaryzacyjnego. Dodatkowo w wiązce światła wychodzącej 
z obiektywu mikroskopowego umieszczano fotodiodę. Oba sygnały wprowadzano 
na rejestrator x-y uzyskując termogram, czyli zależność temperaturową intensyw- 
ności światła wychodzącego z próbki, opisujący procesy przejść fazowych. Taki ter- 
mogram ilustruje Rysunek 6b [49]. 
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Rysunek 6. (a) Schemat aparatury Baralla II do badań termooptycznych (TOA). Mikroskop (A), fotokomórka 
(B), polaryzator (C), analizator(D), stolik grzewczy (E), regulator temperatury (F), rejestrator x-y 
(G), wzmacniacz fotokomórki (H), próbka (1), filtr (J), źródło światła (K) [49]. (b) Termogram 
TOA anizaldazyny wraz termogramem DTA wg Baralla II [49] 

Figure 6. (a) Diagram of Barall IP apparatus for thermo-optic research (TOA). Microscope (A), photocell 
(B), polarizer (C), analyzer (D), heating table (E), temperature controller (F), recorder xy (G), 
photocell amplifier (H), sample (I), filter (J), the light source (K) [49]. (b) Thermogram of TOA 
anizaldazine with a DTA thermogram according to Barall II [49] 
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Postęp elektroniki pozwolił zastąpić rejestrator x-y komputerem, co w spo- 
sób radykalny ułatwiło rejestrację cyfrową wyników oraz umożliwiło stosowanie 
metod statystycznych do interpretacji obserwowanych procesów. Urządzenie do 
takich obserwacji było nawet produkowane komercyjnie jako dodatkowy osprzęt 
do stolika grzewczego typu Mettler F80. Metoda ta została również uznana za jedną 
z metod analizy termicznej i opisana w klasycznym dziele „Handbook of Calori- 
metry Vol. I” [50]. 

Obecnie technika ta jest stosowana w bardzo wielu laboratoriach ciekłokry- 
stalicznych [52-55]. W szczególności w badaniach efektów elektrooptycznych 
faz ciekłokrystalicznych. Wiele takich pomiarów wykonano np. w laboratorium 
prof. Colesa na Uniwersytecie w Manchester [51]. Obecnie w Polsce powszechnie 
tę technikę stosuje, jako metodę podstawową i komplementarną do DSC, laborato- 
rium prof. Chruściela na Uniwersytecie Podlaskim [55] oraz Laboratorium Ciekło- 
krystaliczne na Uniwersytecie Wrocławskim [53]. 
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Rysunek 7. Schemat budowy aparatury TOA stosowane w Laboratorium Cieklych Kryształów na 
Uniwersytecie Wrocławskim [53, 62] 

Figure 7. Construction diagram of TOA equipment used in the Laboratory of Liquid Crystals at the 
University of Wroclaw [53, 62] 
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Rysunek 8. Termogram TOA (E)-4-((4-dodecyloksyfenylo)diazenylo)benzoesanu tetradecylu 
Figure 8. TOA thermogram (E) tetradecyl 4-((4 -dodecylphenyl)diazenyl)benzoate 


W dzisiejszych czasach zestaw do pomiarów TOA można zbudować w nie- 
zwykle prosty sposób (Rys. 7). Wymagany jest optyczny mikroskop polaryzacyjny 
(POM), stolik grzewczy z regulatorem temperatury oraz multimetr. Stolik grzewczy 
może być własnej konstrukcji lub komercyjny (np. LINKAM, METTLER, INSTEC). 
W światło mikroskopu możemy wprowadzić fotodiodę (np. OSRAM BPW 21), 
którą podłączamy do multimetru cyfrowego. Multimetr ten jak również regulator 
temperatury stolika grzewczego łączymy z komputerem i za pomocą odpowied- 
niego oprogramowania rejestrujemy zależność temperaturową intensywności 
światła przechodzącego przez mikroskop. Schemat takiego zestawu stosowanego 
na Uniwersytecie Wrocławskim przedstawia Rysunek 7 [53], a typowy termogram 
Rysunek 8. Przedstawia on właściwości termiczne (E)-4-((4-dodecyloksyfenylo) 
diazenylo)benzoesanu tetradecylu. 

Rysunek 8 wyraźnie opisuje właściwości termiczne związku ciekłokrystalicz- 
nego o stosunkowo prostej budowie. Przede wszystkim wyraźnie widać, że wystę- 
pujące dwie mezofazy są monotropowe. Punkty czerwone opisują proces topnie- 
nia. Tylko 3 punkty przedstawiają moment współistnienia dwóch faz: krystalicznej 
i izotropowej. Ta zależność jest jedną z najważniejszych cech tej metody badawczej 
(TOA) i pozwala potwierdzić czystość badanego związku, w tym również brak śla- 
dów rozpuszczalnika. W metodzie tej możemy przyjąć, że związek jest rzeczywiście 
czysty, gdy różnica temperatur pomiędzy punktami A i A zawiera się w przedziale 
0,3 K. Najważniejszą cechą tej metody jest to, że te pomiary możemy przeprowadzać 
bez utraty dokładności mierzonego sygnału nawet przy niezwykle małych szyb- 
kościach „skanu temperatury” np. 0,2 K/min eliminując w ten sposób problemy 


ANALIZA TERMO-OPTYCZNA (TOA) JAKO NARZEDZIE BADAN ZJAWISK TOPNIENIA 217 


przewodnictwa ciepła. W technice tej jest możliwe również automatyczne wyzna- 
czenie temperatury topnienia w postaci tzw. „onset”-u, powszechnie stosowanego 
w metodach DSC. Podczas chłodzenia (niebieskie punkty) wyraźnie widać, że około 
1 Kponiżej temperatury topnienia pojawia się pierwsza mezofaza, która została jed- 
noznacznie zidentyfikowana metodą tekstur ciekłokrystalicznych (POM) jako faza 
smektyczna A. W tym przypadku niska wartość intensywności światła jest związana 
z powstaniem tzw. tekstury homeotropowej, czyli prostopadłego ułożenia molekuł 
względem powierzchni stolika, a równoległego względem padającej wiązki światła. 
Około 82°C pojawia się faza smektyczna C. Jest to bardzo często przejście II rodzaju, 
niewidoczne metodą DSC. Natomiast w tej geometrii pomiaru przemiana fazowa 
jest widoczna niezwykle wyraźnie. Liniowy wzrost intensywności jest związany 
z powstaniem tekstury smugowej, w miejsce homeotropowej, z pojawieniem się 
coraz większego kontrastu. Ten proces jest prawie liniowy i dobrze opisany linią DF 
na rysunku. Linia E przedstawia proces krystalizacji. Bardzo często ten proces nie 
jest tak ostry, ale zachodzi w pewnym zakresie temperatur. Widać wtedy wyraźnie 
tzw. front fazowy. Natomiast w przypadku faz przechłodzonych, będących w stanie 
metatrwałym, przejście do fazy krystalicznej, stabilnej termodynamicznie, odbywa 
się bardzo szybko co widzimy jako ostra przemiana fazowa. 

W pomiarach tych nie jest istotna wartość obserwowanej intensywności światła, 
ale jej zmiana. Dlatego być może lepszym opisem byłoby rejestrowanie temperatu- 
rowej zależności pochodnej intensywności po temperaturze (podobnie jak w meto- 
dzie EPR). Nie wnosi to jednak istotnych danych w badaniu przemian fazowych. 
Nawet rejestracja zwykłej intensywności światła przechodzącego przez mikroskop 
polaryzacyjny jest wystarczająca do właściwej interpretacji sytuacji fazowej. 

Należy jednak zauważyć, że dzięki powszechnemu zastosowaniu kamer cyfro- 
wych w badaniach ciekłokrystalicznych przemian fazowych w literaturze obser- 
wujemy próby stosowania technik statystycznej obróbki obrazu [56-59]. Bardzo 
ciekawe jest sprzężenie pomiaru termooptycznego (TOA) z jednoczesnym użyciem 
orientującego zmiennego pola elektrycznego. Taki eksperyment może być wykorzy- 
stany do badania polimorfizmu mezofaz [60, 61]. 


2.2. ZASTOSOWANIE METODY TOA W BADANIACH PROCESÓW TOPNIENIA 


Metoda TOA charakteryzuje się następującymi właściwościami: 

« Pozwala dokładnie wyznaczyć temperaturę przemiany fazowej zarówno 
I rodzaju jak i II rodzaju pod warunkiem występowania wyraźnych zmian 
tekstury w świetle spolaryzowanym. Daje również możliwość automatyzacji 
tej procedury. 

u Umożliwia badanie przemian fazowych przy bardzo małych szybkościach 
grzania i chłodzenia (w metodzie DSC ta procedura znacząco obniża war- 
tość sygnału kalorymetrycznego). 
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= Umożliwia dokładnie wyznaczać zakres występowania faz wąskozakreso- 
wych (np. AT=1°C) poprzez obniżenie szybkości skanu do 1°C/min oraz 
niżej. Nie jest to możliwe w metodzie DSC. 

u Daje możliwość badania (przy odpowiedniej geometrii tekstury) mezogen- 
nych przejść fazowych, również tych II rodzaju, niewidocznych metodami 
kalorymetrii DSC (np. SmA-SmC przy teksturze homeotropowej). 

= Pozwala dokładnie wyznaczyć zakres obszaru dwufazowego, co umożliwia 
automatyzację określania rzeczywistej temperatury przemiany fazowej jak 
również określenia czystości badanego związki. 

e Podstawową wadą tej metody jest brak powtarzalności wartości sygnału 
(intensywność przechodzącego światła) dla tej samej substancji w innej 
celce pomiarowej. 


PODSUMOWANIE 


Metoda termooptyczna (TOA) jest niezwykle przydatną metodą szybkiej reje- 
stracji (i często identyfikacji) przemian fazowych w fazie skondensowanej. Przy 
codziennym, rutynowym jej stosowaniu staje się niezwykle cenna, w szczególności 
w badaniach faz ciekłokrystalicznych. Najważniejsze, że jest to metoda komplemen- 
tarna do stosowanych powszechnie w tym celu metod kalorymetrycznych (głównie 
metod DSC). Ponieważ w badaniach przemian fazowych najczęściej dany proces 
chcielibyśmy zrozumieć pod kątem mechanizmu pojawiającego się zjawiska, dla- 
tego tę metodę (TOA) możemy traktować jako komplementarny do metody kalory- 
metrycznej DSC sposób potwierdzenia obecności przemiany fazowej. Na podstawie 
tych dwóch metod możemy projektować dalsze metody badań, często dużo bardziej 
czasochłonne lub wymagające znacznie większej ilości badanego materiału. 
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ABSTRACT 


Palladium immobilized in metal-organic frameworks (MOF) exhibit promis- 
ing catalytic properties in hydrogenation of different unsaturated substrates. Due to 
the specific porous and crystalline structure MOFs can contribute in bonding and 
activation of organic substrates, increasing catalytic efficiency of PA@MOF compos- 
ites. The superior tunability of MOFs structures enables to design highly selective 
catalysts for hydrogenation of different substrates, such as olefins, esters, ketones, 
alcohols or alkynes. Due to the synergistic effects of palladium and MOF not only 
high activity but also high selectivity can be achieved. 

The article presents representative examples of MOF-based palladium catalysts 
for hydrogenation to illustrate perspectives, also technological, of their application. 


Keywords: metal-organic frameworks, catalysis, palladium, hydrogenation, reaction 
selectivity 
Stowa kluczowe: sieci metaliczno-organiczne, kataliza, pallad, uwodornienie, selek- 


tywność 
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WYKAZ STOSOWANYCH SKROTOW 


MOF - ang. Metal Organic Framework 

XAS - ang. X-ray absorption spectroscopy 
XRD - ang. X-Ray Diffraction 

MIL - ang. Material of Institut Lavoisier 

ZIF - ang. Zeolitic Imidazolate Framework 
UiO - ang. University of Oslo 

PVP - poliwinylopirolidon 


TTACP — tetratio-azacyklopentadekan 
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WPROWADZENIE 


Sieci metaliczno-organiczne typu MOF (ang. Metal Organic Frameworks) sa 
to krystaliczne i porowate trójwymiarowe polimery koordynacyjne, zbudowane 
z jonów metali, lub częściej ich klastrów połączonych ze sobą wiązaniem koordy- 
nacyjnym za pośrednictwem organicznego liganda [1-4]. W strukturze sieci meta- 
liczno-organicznych MOF znajdują się dwa centralne elementy budulcowe, którymi 
są węzły i łączniki. Węzłami są zwykle klastry lub jony metali przejściowych, metali 
głównych (np. Al, Sn) lub lantanowców [5]. Łącznikami mogą być związki orga- 
niczne (ligandy) posiadające co najmniej dwie grupy funkcyjne zdolne do tworzenia 
wiązań kowalencyjnych z jonami metalu „na obu końcach” [1, 2]. Do syntezy sieci 
typu MOF najczęściej używane są związki zawierające grupy karboksylowe, sulfo- 
nowe i fosfonowe, a także imidazolowe czy pirazolowe [6]. 

W porównaniu do innych materiałów porowatych, jak zeolity, krzemionka czy 
węgle aktywne, materiały typu MOF wykazują szereg charakterystycznych cech 
ważnych z punktu widzenia ich potencjalnych zastosowań. Wśród tych cech warto 
wymienić dużą powierzchnię właściwą, dużą objętość porów czy niską gęstość [7], 
które otwierają możliwości wykorzystania związków typu MOF w katalizie. Duża 
powierzchnia właściwa wynosząca często ponad 1000 mř/g, i znaczna objętość 
porów ułatwiają dyfuzję substratów do centrów katalitycznie aktywnych, co może 
zwiększać efektywność reakcji. Duży potencjał materiałów typu MOF jako nośni- 
ków katalizatorów metalicznych wiąże się z szerokimi możliwościami ich mody- 
fikacji strukturalnych poprzez zmianę rozmiaru porów i kształtowanie otoczenia 
chemicznego centrów aktywnych wewnątrz porów dzięki doborowi odpowiednich 
bloków budulcowych lub modyfikacje post-syntetyczne. Obecność porów o odpo- 
wiednich rozmiarach i charakterze chemicznym umożliwia związanie kompleksów 
metali wewnątrz struktury MOF, o ile pozwalają na to rozmiary ich cząsteczek. 
Ponadto, używając materiał o wielkości porów dobranej do rozmiarów substra- 
tów reakcji katalitycznej można sterować jej selektywnością [3]. W tym miejscu 
warto zauważyć, że struktura związków typu MOF jest znacznie bardziej elastyczna 
w porównaniu do zeolitów czy tlenków metali. Oznacza to, że materiały typu MOF 
po zaadsorbowaniu cząsteczek związku chemicznego mogą się kurczyć lub roz- 
szerzać [9, 10]. Kolejną cechą, która może decydować o potencjalnej przydatności 
związków typu MOF w katalizie jest ich homochiralność, co otwiera możliwość roz- 
wijania heterogenicznej katalizy asymetrycznej [11]. Ponadto, wpływ na aktywność 
katalityczną mogą mieć grupy funkcyjne ligandów, jako miejsca koordynacji aktyw- 
nych katalitycznie jonów metali lub nanocząstek metali. 

Reakcja katalityczna z udziałem związków typu MOF jako nośników kataliza- 
torów metalicznych może przebiegać na dwa sposoby przedstawione schematycznie 
na Rysunku 1. W pierwszym przypadku za aktywność katalityczną odpowiadają jony 
metali (lub klastry) obecne w węzłach polimeru koordynacyjnego MOF (Rys. 1a). 
Drugą, częściej spotykaną możliwością, jest immobilizowanie kompleksów metali 
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lub nanocząstek metali (Rys. 1b) wewnątrz struktury związku typu MOF lub na jego 
powierzchni [12]. 


Produkty 
a) b) 


Substraty 


Produkty 


Substraty 


ży Klaster lub jon metalu w Immobilizowane kompleksy lub nanocząstki metali 


Rysunek 1. Centra aktywne katalitycznie w polimerach koordynacyjnych MOF 
Figure 1. Catalytic centers in MOF polymers 


1. WŁAŚCIWOŚCI KATALITYCZNE [Pd(2-PYMO)], 


Przykładem związku typu MOĘ w którym reakcja katalityczna przebiega na 
węźle metalicznym jest związek [Pd(2-pymo), |, gdzie 2-pymo (pymo = jon 2-piry- 
midynolanowy) jest ligandem, w którym atomy azotu koordynują do jonu Pd” two- 
rząc trójwymiarową porowatą strukturę [13]. W 2007 r. w zespole Cormy przepro- 
wadzono reakcję uwodornienia 1-oktenu, stosując [Pd(2-pymo)], jako katalizator 
i uzyskano 99% konwersję substratu z 59% selektywnością do n-oktanu w czasie 
40 min. Drugim produktem reakcji był 2-okten, utworzony w wyniku izomeryzacji 
1-oktenu. Zastosowanie cyklododekenu jako drugiego substratu w reakcji uwodor- 
nienia pozwoliło na pokazanie zależności między wielkością substratu a jego reak- 
tywnością. Cyklodeken nie ulegał reakcji uwodornienia, co doprowadziło autorów 
do wniosku, że jest on zbyt duży aby dostać się do wnętrza struktury [Pd(2-pymo),],. 
Tym samym założono, że aktywne centra palladowe są zlokalizowane wewnątrz- 
porów. [14] Dalsze badania nad mechanizmem reakcji uwodornienia przedstawił 
w 2012 roku Klemm [15, 16]. Rozważał on dwie możliwe drogi reakcji, różniące się 
mechanizmem aktywacji cząsteczki H. odpowiednio homolitycznym lub hetero- 
litycznym. W pierwszym przypadku utleniająca addycja cząsteczki H, do palladu 
prowadzi do utworzenia diwodorkowego kompleksu Pd** (1) (Rys. 2) do którego 
koordynuje alken (2). Po tym etapie następuje migracyjna insercja i utworzenie 
kompleksu alkilowego oraz ponowna koordynacja liganda (3). W wyniku redukcyj- 
nej eliminacji następuje odszczepienie cząsteczki produktu i powrót do pierwotnej 
formy Pd” katalizatora (4). 
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Rysunek 2. Mechanizm reakcji uwodornienia alkenu katalizowanej przez [Pd(2-pymo),]„ z homolitycznym 
rozszczepieniem wodoru 
Figure 2. Mechanism of alkene hydrogenation catalyzed by [Pd(2-pymo),], with homolytic splitting of H. 
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Rysunek 3. ` Mechanizm reakcji uwodornienia alkenu katalizowanej przez [Pd(2-pymo),], z heterolitycznym 
rozszczepieniem wodoru. 
Figure 3. Mechanism of alkene hydrogenation catalyzed by [Pd(2-pymo),] with heterolytic splitting of H. 


228 A.W. AUGUSTYNIAK, A.M. TRZECIAK 


Mechanizm oparty na heterolitycznym rozczepieniu czasteczki H, (Rys. 3) 
zakłada, ze po związaniu cząsteczki H, wiązanie H-H zostaje rozerwane w ten 
sposób, ze ligand ulega protonacji i powstaje monowodorkowy kompleks Pd* (1). 
Następnie ma miejsce przyłączenie alkenu i migracyjna insercja (2-3). W końco- 
wym etapie katalizator przyłącza drugą cząsteczkę wodoru, następuje usunięcie 
alkanu ze sfery koordynacyjnej i odtwarza się wodorkowy kompleks Pd”'(4). 

W celu uzyskania bardziej szczegółowych danych dotyczących mechanizmu 
reakcji uwodornienia oraz aktywnej formy palladu przeprowadzono pomiary in situ 
za pomocą rentgenowskiej spektroskopii absorpcyjnej XAS (ang. X-ray Absorption 
Spectroscopy) wspomagane przez pomiary dyfraktometrii proszkowej XRD (ang. 
X-ray Powder Diffraction) oraz spektroskopii w podczerwieni IR (ang. infrared 
spectroscopy). Stwierdzono, że ligand pymo ulega protonowaniu już po 30 min eks- 
pozycji na H,. Później, w ciągu 3 h następuje stopniowe uwodornienie pierścienia 
aromatycznego liganda. Redukcja Pd”* do Pd’ w węźle [Pd(2-pymo) |], rozpoczyna 
się po ok. 3-5 h działania H,, co może wskazywać na to, że jest ona indukowana 
uwodornieniem i protonowaniem liganda. Równocześnie ma miejsce zniszczenie 
struktury związku MOF jednak układ jest nadal aktywny katalitycznie dzięki obec- 
ności nanoczastek Pd’. Dlatego po dłuższym czasie, 4 h, cyklodeken także ulega 
uwodornieniu. Tym samym destrukcja struktury [Pd(2-pymo),], powoduje zanik 
obserwowanej początkowo selektywności związanej z rozmiarem porów i substra- 
tów. Uwodornienie 1-oktenu następuje prawdopodobnie na centrach Pd”*, ponieważ 
pierwsze formy Pd” są identyfikowane dopiero po zakończonej reakcji uwodornie- 
nia. Nie udało się jednak zidentyfikować kompleksów wodorkowych Pd* stosując 
pomiary XAS, prawdopodobnie ze względu na ich niskie stężenia. 


2. SYNTEZA KATALIZATORÓW M©MOF 


W literaturze można znaleźć wiele przykładów, w których porowaty materiał 
typu MOF spełnia rolę nośnika katalizatora metalicznego [17, 18]. Znane są różne 
metody immobilizacji związków metali na materiałach typu MOF, które pozwalają 
otrzymać kompozyty MGMOF Można tu wyróżnić metodę ex situ polegającą na 
impregnacji związków typu MOF kompleksami metali lub nanocząstkami metali. 
Druga metoda, in situ, polega na tworzeniu materiału typu MOF z prekursorów 
w roztworze, w którym znajduje się katalizator w postaci kompleksu metalu lub 
nanocząstek. W tym miejscu warto podkreślić, że nanocząstki metalu mogą również 
powstać in situ zarówno w trakcie syntezy kompozytu jak również podczas reakcji 
katalitycznej. 
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Rysunek 4. Sposoby syntezy materiałów kompozytowych typu M@MOF 
Figure 4. Methods of synthesis of M@MOF composites 


3. KATALIZATORY DWUFUNKCYJNE 


Atrakcyjną cechą katalizatorów immobilizowanych na materiałach typu MOF 
jest ich potencjalna dwufunkcyjność, która umożliwia równoczesną aktywację 
dwóch substratów. Przykładem takiego katalizatora jest PA@MIL-101, zastosowany 
do syntezy ketonu metylowo-izobutylowego z acetonu i H, [19]. Wykonanie syntezy 
ketonu w jednym etapie było możliwe dzięki połączeniu funkcji kwasu Lewisa, którą 
pełnił MIL-101 z aktywnością nanocząstek Pd’ w reakcji uwodornienia. 

Zastosowanie dwufunkcyjnego katalizatora PA@MIL-101 zawierającego niena- 
sycone koordynacyjnie centra o charakterze kwasu Lewisa (Cr**) oraz nanocząstki 
Pd” umożliwiło opracowanie jednoetapowej transformacji cytronellalu do mentolu 
(Rys. 5) [20]. Pierwszy etap reakcji prowadzony w atmosferze N,, polegał na cykli- 
zacji cytronelalu do izospugenolu i nie wymagał obecności palladu. Materiał typu 
MOE, Cr '-MIL-101, umożliwiał wykonanie cyklizacji, chociaż reakcja była wolniej- 
sza niż w obecności kompozytu P€@MIL-101. Drugi etap wymagał atmosfery H, 
oraz udziału palladu jako katalizatora uwodornienia wiązania C=C. Warto zauwa- 
żyć, że zastosowanie H, od początku reakcji powodowało spadek selektywności do 
mentolu, zamiast którego powstawał głównie 3,7-dimentylooktanol (Rys. 5). Można 
to wyjaśnić większą szybkością uwodornienia wiązań C=C i C=O, podczas gdy izo- 
meryzacja przebiegała w tych warunkach wolniej. 
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Rysunek 5. Schemat reakcji cyklizacji oraz uwodornienia cytronelalu 
Figure 5. Scheme of the cyclization and hydrogenation reactions of citronellal 


Dwufunkcyjny katalizator PA@MIL-101 użyto w wieloetapowej syntezie tan- 
demowej, w której substratami były związki nitrowe a produktami drugorzędowe 
aryloaminy (Rys. 6) [21]. Zaletą opracowanej metody było współdziałanie w reakcji 
uwodornienia kwasowości Lewisa materiału MIL-101 i aktywności katalitycznej 
nanoczastek Pd”. Katalizator dwufunkcyjny okazał się lepszy niż katalizatory pal- 
ladowe osadzone na nośnikach o znikomej kwasowości, takich jak węgiel aktywny 
czy ALO.. 
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Rysunek 6. _ Tandemowa synteza aryloamin 
Figure 6. Tandem synthesis of arylamines 


4. NANOCZĄSTKI PALLADU IMMOBILIZOWANE 
W MATERIAŁACH TYPU MOF 


Do tej pory dużo zastosowań w katalizie znalazły układy zawierające 
nanocząstki metali związane z materiałem typu MOF (nanocząstki metalu@ 
MOF), które można podzielić na dwie grupy [22]. Pierwszą grupę stano- 
wią kompozyty, w których średnica nanocząstek jest większa niż rozmiar 
porów materiału typu MOF i są one osadzone na powierzchni polimeru. 
W drugiej grupie układów wielkość nanocząstek jest dopasowana do rozmiaru 


KATALIZATORY PALLADOWE IMMOBILIZOWANE W MATERIALACH TYPU MOF 231 


porów materiału typu MOF i znajdują sie one wewnątrz porowatej struktury poli- 
meru (Rys. 7). 


Pory MOF 


E | Nanocząstki metali Powierzchnia MOF 


Rysunek 7. Rodzaje kompozytów typu nanocząstki metalu@MOF 
Figure 7. Examples of metal nanoparticles@MOF composites 


Na aktywność i selektywność takich układów katalitycznych mogą wpływać 
różne czynniki, takie jak wielkość, kształt, dyspersja oraz położenie nanocząstek 
metalu, a także charakter i budowa samego nośnika. 

Przykładem katalizatora, w którym nanocząstki Pd są związane wewnątrz 
porów jest katalizator PAQZIF-8 (ang. Zeolitic imidazolate framework) [23]. Synteza 
tego kompozytu przeprowadzona była w roztworze prekursorów materiału typu 
MOF i wcześniej przygotowanych nanocząstek palladu stabilizowanych PVP. Co 
ważne, nanocząstki palladu w trakcie syntezy znalazły się wewnątrz porów, co bar- 
dzo często jest niemożliwe dla wstępnie przygotowanych nanocząstek ze względu 
na ich większy rozmiar w stosunku do wielkości porów związku MOE Pd@ZIF-8 
okazał się dobrym katalizatorem reakcji redukcji p-nitrofenolu do p-aminofenolu 
co może być związane z bardzo dobrą dyspersją nanocząstek palladu znajdujących 
się w porach ZIF-8 (Rys. 8). 


Rysunek 8. Obrazy TEM nanoczastek palladu immobilizowanych w ZIF-8 
Figure 8. TEM pictures of the Pd nanoparticles immobilized in ZIF-8 


232 A.W. AUGUSTYNIAK, A.M. TRZECIAK 


Pd@MOF-5 jest przykładem katalizatora, w którym nanoczastki palladu znaj- 
dują się na powierzchni związku typu MOF [24]. Aktywność katalityczną tego 
materiału sprawdzono w reakcji uwodornienia cynamonianu etylu. Interesujące 
jest to, że badania sorpcji N, wykazały utratę porowatości już po pierwszym użyciu 
katalizatora, który mimo to wykazywał większą aktywność w kolejnych 2 cyklach 
niż w pierwszym. Można wnioskować, że nanocząstki palladu prawdopodobnie nie 
znajdują się wewnątrz porów MOF-5, a są związane na jego powierzchni. 

W naszej grupie badawczej opracowano katalizatory oparte na związku MOF 
tytanu(IV), MIL-125-Ti-NH, (MIL= ang. Material of Institut Lavoisier) oraz nano- 
cząstkach palladu (Rys. 9) [25]. 


Pd(OAc), ,MeOH 
30 min, RT 


Rysunek 9. Schemat syntezy i obrazy TEM nanoczastek palladu immobilizowanych w MIL-125-Ti-NH, 
Figure 9. Scheme of the synthesis and TEM pictures for the Pd nanoparticles immobilized in 
MIL-125-Ti-NH, 


Synteze prowadzono w lagodnych warunkach, a dodatek KBr pozwolit otrzy- 
mać nanoczastki palladu o regularnych kształtach kubicznych (Rys. 9). Badania 
aktywności katalitycznej tych związków pokazały, że szybkość reakcji redukcji 
p-nitrofenolu zależy od morfologii nanocząstek palladu. Dla nanocząstek o kształ- 
tach kubicznych uzyskano wyższe stałe szybkości niż dla nanocząstek o mniej zde- 
finiowanych kształtach. Większa aktywność katalizatora zawierającego kubiczne 
nanocząstki Pd spowodowana jest prawdopodobnie większą powierzchnią aktywną 
Pd’, co zostało potwierdzone eksperymentalnie za pomocą pomiaru chemisorpcji 
CO. Ważnym rezultatem przedstawionych badań jest stwierdzenie, że katalizatory 
na bazie MOF zawierające nanocząstki Pd można otrzymać w łagodnych warun- 
kach. Ponadto, stosując odpowiednie czynniki modyfikujące można wpływać na 
kształt nanocząstek, a w konsekwencji także na ich aktywność katalityczną. 

Katalizator bimetaliczny zawierający nanocząstki Pd i Ni otrzymano wyko- 
rzystując jako nośnik związek UiO-66-NH, (UiO = University of Oslo) funkcjona- 
lizowany tetratio-.azacyklopentadekanem (TTACP) (Rys.10) [26]. Funkcjonaliza- 
cja była możliwa dzięki obecności grup N-H na powierzchni nośnika, a cykliczny 
ligand, znajdujący się na zewnątrz struktury związku MOF posłużył do immobili- 
zacji nanocząstek palladu i niklu o strukturze typu rdzeń-otoczka (ang. Coreshell). 
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Ui0-66-NH,/TTACP Ui0-66-NH,/TTACP/NiQPd 


Rysunek 10. Schemat syntezy kompozytu UiO-66-NH,/TTACP/Ni@Pd 
Figure 10. Scheme of UiO-66-NH,/TTACP/Ni@Pd composite synthesis 


Aktywność katalityczną UiO-66-NH, /TTACP/Ni@Pd badano w reakcji reduk- 
cji 4-nitrofenolu i 2-nitroaniliny do odpowiednich związków aminowych, stosując 
NaBH, jako donor wodoru. Parametry kinetyczne pokazały, że stała szybkości reak- 
cji dla tego katalizatora (k = 0,0142 s™) jest większa niż dla większości badanych 
wcześniej katalizatorów. Dodatkową zaletą są właściwości magnetyczne kataliza- 
tora związane z obecnością niklu, które ułatwiają jego wydzielenie z mieszaniny 
reakcyjnej przy pomocy zwykłego magnesu. Reakcję z tą samą próbką katalizatora 
powtórzono dziesięciokrotnie, bez znaczącego spadku konwersji. Zgodnie z zapro- 
ponowanym schematem aktywacji substratów na powierzchni katalizatora, kluczo- 
wym etapem jest adsorpcja BH, a następnie rozerwanie wiązania B-H i utworzenie 
wodorku Pd-H. 

Li opisał jednoetapową syntezę katalizatora PA@UiO-67, zawierającego nano- 
cząstki Pd” enkapsulowane w strukturze MOF [27]. Do syntezy prowadzonej in situ 
użyto ZrCl, i kwas 2,2 '-bipirydyno-5,5 -dikarboksylowy oraz PdCl,(CH,CN),. 
Zastosowanie jako łącznika kwasu z fragmentem bipirydynowym zamiast kwasu 
4,4,-difenylodikarboksylowego stosowanego w syntezie UiO-67 miało ułatwić 
koordynację palladu. Dalsze modyfikacje warunków syntezy, pokazały, że tempera- 
tura wpływa na wielkość nanocząstek Pd” i ich umiejscowienie w strukturze MOF. 
Otrzymane katalizatory zawierające Pd” były aktywne w uwodornieniu styrenu do 
etylobenzenu, nawet w temperaturze pokojowej, pod ciśnieniem 1 atm H,. Ponadto, 
łatwo można je wydzielić z mieszaniny reakcyjnej i ponownie użyć. Natomiast ana- 
logiczne kompozyty z Pd” nie wykazały żadnej aktywności katalitycznej w zastoso- 
wanych warunkach. 

Podobny sposób syntezy in situ zastosowano z powodzeniem do otrzymania 
katalizatora bimetalicznego, PANi-in-UiO-67, w którym nanocząstki metali znajdo- 
wały się w porach sieci MOF [28]. Prezentując tę syntezę, autorzy podkreślili ważną 
rolę liganda 2,2 -bipirydyno-5,5 -dikarboksylowego, posiadającego dwa rodzaje 
grup funkcyjnych w rozumieniu teorii miękkich i twardych kwasów i zasad. Grupy 
karboksylanowe, o charakterze twardych zasad silne wiążą się z oksofilnymi katio- 
nami Zr“ tworzącymi strukturę MOF, Z kolei miękkie grupy bipirydylowe wiążą 
prekursory metali, Pd i Ni. W celu zbadania synergicznego efektu katalitycznego 
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działania dwóch metali (Pd, Ni) w reakcji uwodornienia nitrobenzenu diwodorem, 
otrzymano kompozyty różniące się składem. Dla porównania autorzy zbadali rów- 
nież katalizatory mono-metaliczne, jednak katalizator, który zawierał tylko nano- 
cząstki niklu (Ni-in-UiO-67) nie był aktywny w tej reakcji. Natomiast katalizator 
zawierający nanocząstki palladu (Pd-in-UiO-67) pozwalał osiągnąć całkowitą kon- 
wersję nitrobenzenu do aniliny dopiero po 18 godzinach. Katalizatory bimetaliczne 
były znacznie bardziej aktywne, zapewniając całkowitą konwersję nitrobenzenu 
w ciągu 1,5-3,5 h w zależności od składu. Autorzy podkreślili korzystne cechy 
katalizatora i powiązali je z umiejscowieniem nanocząstek wewnątrz sieci MOF. 
Podkreślono wysoką aktywność, a także trwałość katalizatora w kolejnych cyklach. 
Obserwowana aktywność katalityczna jest wyższa niż aktywność podobnych katali- 
zatorów, w których aktywne metale znajdują się na powierzchni nośnika typu MOF. 

Odpowiedni dobór materiału typu MOF jako nośnika katalizatora może mieć 
kluczowe znaczenie dla selektywności reakcji katalitycznej. W szczególności, wiel- 
kość porów może ograniczać dostęp zbyt dużych substratów do centrów aktywnych 
i przez to wpływać na selektywność reakcji. Ilustracją tej zależności są wyniki uzy- 
skane w reakcji uwodornienia alkenów katalizowanej przez Pd@ZIF-L, w którym 
nanocząstki palladu znajdują się wewnątrz porowatej struktury (Rys. 11) [29]. Śred- 
nice kinetyczne 1-heksenu, cykloheksenu i cyklooktenu wynoszą odpowiednio 1,7, 
4,2 i 5,5 A i są mniejsze niż wielkość porów ZIF-L (6,64 A), co umożliwia ich kon- 
wersję do produktów uwodornionych. Natomiast brak reakcji tetrafenyloetylenu 
może być związany z jego większym rozmiarem (6,7 A). 


d 


Konwersja % 


Rysunek 11. Konwersja różnych alkenów w reakcji uwodornienia katalizowanej przez PAQZIF-L 
Figure 11. Conversion of different alkenes in hydrogenation catalyzed by PA@ZIF-L 
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5. ZASTOSOWANIE KATALIZATOROW PD@MOF 
W PRZEMYSLOWYCH REAKCJACH UWODORNIENIA 


W Tabeli 1 przedstawiono przykłady zastosowań katalizatorów Pd@MOF 
w reakcjach uwodornienia różnych substratów, także w procesach uwodornie- 
nia o znaczeniu technologicznym. Dobre wyniki uzyskane w tym obszarze mogą 
w przyszłości wpłynąć na zmianę procesu produkcji i wprowadzenie katalizatorów 
immobilizowanych na materiałach typu MOF do przemysłu. 


Tabela 1. Katalizatory PA@MOF zastosowane w reakcjach uwodornienia 
Table 1. Pd@MOF catalyst used in hydrogenation reaction 
Katalizator Substrat 

Pd@ZIF-8 p-nitrofenol 
Pd@Zn-BDC p-nitrofenol 
Pd@MIL-101-Cr p-nitrofenol 
Pd@MIL-125-Ti-NH, p-nitrofenol 
Ni@Pd@MIL-101 nitrobenzen 
Ni@Pd@UiO-66-NH, p-nitrofenol 
Fe,O,@Pd@MIL-100-Fe p-nitrofenol 
Ni,Pd,@Ui067 nitrobenzen 
Pd@UiO-66 acetofenon/benzofenon 


Pd@UiO-67 


acetofenon/benzofenon 


NiPd@MIL-101 


3-heptanon 


Pd@MIL-101-Al-NH, 


5-hydroksymetylofurfural 


Pd@ZIF-8 alkeny 
PdQZIF-8 1-heksen 
PdQZIF-L 1-heksen 
PdGSiO,QZIF-8 1-heksen 
Pd@ZSM-5@ZIF-8 1-heksen 
Pd@MSS@ZIF-8 1-heksen 
C@Pd@ZIF-8 styren 
Pd@MOF-3 styren 
Pd@MOF-5 styren 
Pd@MIL-101 styren 
Pd@UiO-67 styren 
Pd@MIL-101-Cr-NH, aldehyd cynamonowy 


Pd@ZIF-8 


aldehyd cynamonowy 


236 A.W. AUGUSTYNIAK, A.M. TRZECIAK 


Katalizator Substrat 
Pd-PVP@ZIF-8 1,4-butynodiol 
Pd@MIL-101 fenol 
Pd@MIL-53 fenol 
Pd@MIL-54-Al fenol 


Jedną z reakcji o znaczeniu technologicznym jest uwodornienie 2,3,5-trime- 
tylobenzochinonu do 2,3,5-trimetylohydrobenzochinonu, ważnego substratu do 
produkcji witaminy E (Rys. 12). Katalizatorami stosowanymi w tym procesie są Ni 
Raneya lub metale z grupy platynowców, Pd, Ru, Rh, naniesione na nośniki, nato- 
miast Pd/C jest katalizatorem przemysłowym. Jego stosunkowo niska aktywność 
była inspiracją do poszukiwania innych rozwiązań z wykorzystaniem materiału typu 
MOE MIL-101(Cr), jako nośnika dla nanocząstek Pd” [54]. Immobilizowany katali- 
zator otrzymano metodą impregnacji z roztworu kompleksu Pd(acac),, a następnie 
zredukowano pallad w atmosferze H, uzyskując wysoką dyspersję metalu w porach 
nośnika. Zgodnie z przewidywaniem, aktywność katalityczna w reakcji uwodor- 
nienia wzrastała ze wzrostem zawartości palladu, a dodatkowo zaobserwowano 
pozytywny wpływ nośnika na wydajność reakcji. Może on być związany z hydrofo- 
bowymi właściwościami materiału MOF, które sprzyjają kondensacji substratu na 
powierzchni katalizatora. Ponadto, centra kwasowe Lewisa (Cr'*) nośnika wspoma- 
gają aktywację substratu dzięki oddziaływaniu z fragmentem aromatycznym. Dzięki 
synergicznemu działaniu nośnika i nanocząstek Pd” otrzymano katalizator o wyż- 
szej aktywności i trwałości niż Pd/C. 


o OH 


Pd/MIL-101-Cr 


H, 
o OH 


Rysunek 12. Schemat reakcji uwodornienia 2,3,5-trimetylobenzochinonu 
Figure 12. Scheme of catalytic hydrogenation of 2,3,5-trimethylbenzoquinone 


Dalsze badania procesu uwodornienia 2,3,5-trimetylobenzochinonu prowa- 
dzono z użyciem materiału NH,-UiO-66 funcjonalizowanego grupami aminowymi 
[55]. Wybór tego nośnika był uzasadniony jego wysoką odpornością termiczną 
i chemiczną, obecnością dobrze zdefiniowanych trójwymiarowych kanałów oraz sil- 
nymi właściwościami koordynacyjnymi grup aminowych. Synteza katalizatora pal- 
ladowego wykonana w DME, który był zarówno rozpuszczalnikiem jak i reduktorem 
palladu, doprowadziła do otrzymania materiału o równomiernym rozmieszczeniu 
nanocząstek w nanoporach nośnika. Aktywność katalityczna okazała się wyższa niż 
aktywność komercyjnego Pd/C, a wyjaśnienia tego faktu poszukiwano przy pomocy 
obliczeń teoretycznych dla modelowych układów zawierających różne klastery Pd 
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(do 32 atomów). Wykazano, że obecność grup NH, wzmacnia wiązanie benzochi- 
nonu do katalizatora w porównaniu z niemodyfikowanym nośnikiem UiO-66, co 
wpływa na wyższą aktywność, która nie zmieniała się w kolejnych 8 cyklach. 

Jako nośnik katalizatora palladowego przeznaczonego do uwodornienia 
2,3,5-trimetylobenzochinonu zastosowano także materiał UiO-66(Hf), analog 
stosowanego częściej UiO-66(Zr) [56]. Motywacją do użycia związku hafnu było 
jego wyższe powinowactwo do tlenu, które ilustruje wyższa energia wiązania Hf-O 
(802 kJ.m"') niż Zr-O (776 kJ.m‘'). Zdaniem autorów ta cecha powinna skutkować 
wzmocnieniem charakteru kwasu Broensteda grup u,-OH w podjednostkach budu- 
jących strukturę MOF w porównaniu ze związkiem Zr. Z kolei, jak wykazano wcze- 
śniej, obecność centrów kwasowych korzystnie wspomaga adsorpcję i aktywację 
substratu. Rzeczywiście, kwasowość UiO-66(Hf) okazała się 5-krotnie wyższa niż 
UiO-66(Zr), a badany katalizator wykazał się aktywnością wyższą niż Ni Raneya, 
a także inne katalizatory osadzone na materiałach MOF, PA@MIL-101 i PAGNH,. 
-UiO-66(Zr). 

Uwodornienie fenolu do cykloheksanonu jest ważnym procesem ze względu 
na zastosowanie cykloheksanonu do produkcji kaprolaktamu i kwasu adypinowego. 
Ponadto, uwodornienie fenolu może być wykorzystane w procesie produkcji biopa- 
liw, a właściwie do ich oczyszczania. Biopaliwa zawierają znaczne ilości związków 
aromatycznych z grupami tlenowymi, wśród których ponad 50% stanowią związki 
fenolowe. Ich obecność zwiększa lepkość paliwa i obniża jego trwałość zwiększając 
przy tym właściwości korozyjne. Uwodornienie związków fenolowych do ketonów, 
alkoholi lub węglowodorów podnosi znacząco jakość biopaliwa. 

Do reakcji uwodornienia fenolu użyto katalizatory palladowe naniesione 
metodą impregnacji mokrej na UiO-66 i UiO-66-NH, [57]. Reakcję prowadzono 
w wodzie pod ciśnieniem H, 0,2 MPa, a w produktach identyfikowano cykloheksa- 
nol i cykloheksanon. Przy zastosowaniu Pd@UiO-66 uzyskano całkowitą konwersję 
fenolu do cykloheksanolu w ciągu 2 h w 120°C, co jest bardzo dobrym wynikiem 
w porównaniu ze znanymi katalizatorami. Drugi katalizator, Pd/UiO-66-NH, oka- 
zał się wyraźnie mniej aktywny, zapewniając konwersję fenolu na poziomie ok. 
20%, przy czym ilości cykloheksanolu i cykloheksanonu były podobne. Obserwo- 
waną różnicę aktywności autorzy powiązali z obecnością grup NH,, które wiązały 
pallad na powierzchni nośnika, podczas gdy w UiO-66 pallad mógł znajdować się 
wewnątrz nośnika. W efekcie oddziaływań między Zr i Pd w nośniku UiO-66-NH, 
było słabsze, a sam nośnik ulegał częściowej destrukcji. Różnice strukturalne miały 
wpływ na selektywność reakcji. Katalizator oparty na UiO-66 wykazywał wysoką 
selektywność do cykloheksanolu, a osadzony na UiO-66-NH, był bardziej selek- 
tywny do cykloheksanonu. 

W reakcji uwodornienia fenolu testowano serię związków Pd@AI-MIL-53-X, 
gdzie X= H, OCH,, NH, Cli NO,, a przedstawione wyniki ilustrują wpływ liganda 
tworzącego strukturę MOF na aktywność katalityczną [58]. Do badań wybrano mate- 
riały typu AI-MIL, które są nietoksyczne i odporne na hydrolizę i mogą być uznane za 
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przyjazne dla środowiska, o potencjalnym zastosowaniu przemysłowym. W syntezie 
użyto pochodne kwasu tereftalowego z różnymi podstawnikami. Wszystkie katali- 
zatory były aktywne w uwodornieniu fenolu do cykloheksanonu w 19°C. Konwersja 
fenolu w trakcie 24 h mieściła się w zakresie od 8,4 do 36,3%, z selektywnością do 
cykloheksanonu powyżej 95% we wszystkich przypadkach w kolejności dla odpo- 
wiednich podstawników w ligandzie H (36,3%) > OCH, (30,4%) > NH, (25%) > Cl 
(11,5%) > NO, (8,4%). Ten kierunek zmian okazał się odwrotny do zmian energii 
aktywacji, która była najwyższa dla liganda podstawionego grupą NO,. Obecność 
podstawników zwiększa hydrofobowy charakter powierzchni i powoduje zmniej- 
szenie powierzchni właściwej, co nie jest korzystne dla aktywacji substratu. Warto 
jednak dodać, że w warunkach wyższej temperatury, 90°C, aktywność katalityczna 
była znacznie wyższa, a negatywny wpływ podstawników nie był widoczny. 

Te same katalizatory zastosowano także do uwodornienia fenyloacetylenu do 
styrenu [58]. Ta reakcja ma znaczenie praktyczne przy usuwaniu fenyloacetylenu 
obecnego jako zanieczyszczenie w surowcu styrenowym do procesu polimeryza- 
cji. Istotny jest zatem taki dobór warunków i katalizatora, aby fenyloacetylen ulegał 
uwodornieniu do styrenu, a nie do etylobenzenu. W serii katalizatorów Pd@AI- 
-MIL-53-X z podstawnikami, H, OCH,, NH, uzyskano znacznie lepszą aktywność 
katalityczną niż w obecności Cl i NO,. Katalizator z grupą NH, wykazywał najlepszą 
aktywność katalityczną, z 24,7% konwersją fenyloacetylenu i 98,8% selektywnością 
do styrenu w ciągu 30 minut. W dłuższym czasie konwersja wzrosła do ponad 99%, 
a uzyskane wyniki są lepsze niż wyniki z użyciem katalizatora palladowego osadzo- 
nego na nieorganicznym tlenku. 

Dwa katalizatory hybrydowe, Pd@UiO-66(Hf) i Pd/UiO-66(Hf) zostały 
użyte do częściowego uwodornienia fenyloacetylenu [59]. Katalizator Pd/UiO- 
-66(Hf) otrzymano przez zmieszanie związku MOF z koloidalnymi nanocząstkami 
Pd. Kompozyt Pd@UiO-66(Hf) otrzymano w syntezie solwotrmalnej, w której 
do zawiesiny nanocząstek Pd w DMF dodano substraty potrzebne do utworzenia 
związku MOF, czyli HfCI,.6H,O i ligand łącznikowy. W ten sposób wytworzono 
materiał, w którym związek typu MOF obudował nnaocząstki Pd. Stwierdzono, że 
ten kompozyt wykazywał wyższą selektywność do styrenu (93,2%) podczas gdy 
selektywność drugiego katalizatora była niższa, 84,3%. Warto dodać, że uwodor- 
nienie prowadzono metodą „przeniesienia wodoru”, a źródłem wodoru był borazan 
(H,N.BH,). Kiedy do reakcji użyto gazowy wodór uzyskano całkowitą konwersję 
fenyloacetylenu do etylobenzenu. 


PODSUMOWANIE 


Sieci metaliczno-organiczne MOF to krystaliczne materiały porowate, które 
charakteryzują się wyjątkowymi właściwościami w porównaniu do takich materia- 
łów jak tlenki metali, węgle aktywne czy zeolity. Chemiczna i strukturalna różno- 
rodność, bardzo dobrze rozwinięta powierzchnia właściwa, zdolność do inkorpora- 
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cji kompleksów i nanocząstek metali to niewątpliwie główne zalety odróżniające je 
korzystnie od znanych nośników katalizatorów. W artykule przedstawiono przykład 
związku MOF, [Pd(2-pymo),]_, w którym węzłami są jony Pd(II) wykazujące aktyw- 
ność katalityczną w reakcji uwodornienia. Dla tego związku zaproponowano dwa 
schematy aktywacji wodoru, homolityczny i heterolityczny. W tym drugim przy- 
padku grupa NH, liganda pymo uczestniczyła w reakcji jako akceptor protonu, co 
umożliwiało utworzenie aktywnego kompleksu wodorkowego z wiązaniem Pd-H. 

Przedstawione w niniejszym przeglądzie przykłady kompozytów Pd@MOF 
aktywnych w reakcjach uwodornienia stanowią jedynie niewielką część opisanych 
w literaturze układów katalitycznych tego typu. W artykule omówiono przykłady 
użycia katalizatorów palladowych immobilizowanych w różnych polimerach typu 
MOE, zwracając uwagę na reakcje w których widoczne jest synergiczne działanie 
palladu i nośnika MOF. Najważniejszą korzyścią, którą można uzyskać stosując te 
katalizatory jest wysoka selektywność reakcji uwodornienia. Może ona wynikać 
z relacji między rozmiarem substratu a wielkością porów, która limituje dostęp sub- 
stratu do centrów aktywnych znajdujących się wewnątrz kompozytu. W zależności 
od budowy, materiały typu MOF mogą wpływać na szybkość dyfuzji substratów, 
kontrolując w ten sposób szybkość procesu katalitycznego i jego selektywność. Przy- 
kładem reakcji, w której oddziaływanie substratu z polimerem typu MOF decyduje 
o selektywności jest uwodornienie fenyloacetylenu do styrenu, zamiast do etyloben- 
zenu. Dzięki specyficznym oddziaływaniom z katalizatorem reakcja zatrzymuje się 
na etapie styrenu, który nie ulega dalszemu uwodornieniu. Ta reakcja ma znaczenie 
w przemysłowym procesie oczyszczania surowca do polimeryzacji i selektywnie 
działające katalizatory są szczególnie poszukiwane. Widać, że kompozyty PAU@MOF 
mają potencjał w tym obszarze i mogą w przyszłości znaleźć praktyczne zastosowanie. 
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ABSTRACT 


Aromaticity/aromatic and substituent/substituent effects belong to the most 
commonly used terms in organic chemistry and related fields. They are used for 
more than a century, and so far are the subject of thousands publications a year. 

The quantitative description of the aromaticity of planar z-electron cyclic 
molecules is based on four criteria: (i) they are more stable than their acyclic unsa- 
turated analogues, (ii) bonds have intermediate lengths between those for the single 
and double ones, (iii) external magnetic field induces z-electron ring current, and 
(iv) aromatic systems prefer reactions in which the z-electron structure is preserved. 
conserved. Quantitative characteristics based on these criteria, named as aromati- 
city indices, allow to relate aromaticity to the substituent effect. This latter can be 
described using either traditional Hammett-type substituent constants or characte- 
ristics based on quantum-chemistry. For this purpose, the energies of properly desi- 
gned homodesmotic reactions and electron density distribution are used. In the first 
case, a descriptor named SESE (substituent effect stabilization energy) is obtained, 
while in the second case — cSAR (charge of the substituent active region), which is 
the sum of the charge of the ipso carbon atom and the charge of the substituent. The 
application of these substituent effect descriptors to a set of zr-electron systems, such 
as: benzene, quinones, cyclopenta- and cyclohepta-dienes, as well as some azoles, 
allowed to draw the following conclusions: 

(i) The less aromatic the system, the stronger the substituent influences the 
n-electron system. Highly aromatic systems are resistant to the substituent effect, in 
line with the organic chemistry experience that aromatic compounds dislike reac- 
tions leading to changes in the 7r-electron structure of the ring. 

(ii) Intramolecular charge transfer (resonance effect) is privileged in cases 
where the number of bonds between the electron-attracting and electron-donating 
atoms is even. These effects are much weaker when this number is odd. Classically, it 
may be related to traditional para vs meta substituent effects in benzene derivatives. 
We should note that in electron-accepting groups, such as CN or NO, (and others), 
electron-accepting atoms are second counting from Cipso. 

(iii) In all cases, when the substituent changes number of -electrons in the ring 
in the direction of 4N+2, its aromaticity increases, for example electron-donating 
substituents in exocyclic substituted pentafulvene, or a halogen atom in complexes 
with heptafulvene. 


Keywords: aromaticity, HOMA, NICS, substituent effect, Hammett constants, quan- 
tum chemistry modeling 

Stowa kluczowe: aromatycznos¢, HOMA, NICS, efekt podstawnikowy, stale Ham- 
metta, metody chemii kwantowej 
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ASE 
cc 
CFI 
cSAR 
DRE 
FLU 


HOMA 


HMO 


HOSE 


SESE 


WYKAZ STOSOWANYCH SKROTOW 


aromatyczna energia rezonansu (ang. aromatic stabi- 
lization energy) 

współczynnik korelacji (ang. correlation coefficient) 
indeks przepływu ładunku (ang. charge flow index) 
ładunek aktywnego obszaru podstawnika (ang. charge 
of the substituent active region) 

energia rezonansu Dewara (ang. Dewar resonance 
energy) 

indeks fluktuacji gęstości ładunku (ang. aromatic 
fluctuation index) 

model aromatyczności za pomocą oscylatora harmo- 
nicznego (ang. harmonic oscillator model of aroma- 
ticity) 

teoria orbitali molekularnych Hiickela (ang. Hiickel 
molecular orbital theory) 

model oscylatora harmonicznego w zastosowaniu do 
energii stabilizacji (ang. harmonic oscillator stabiliza- 
tion energy) 

przesunięcie chemiczne niezależne od jądra (ang. 
nucleus independent chemical shift) 

indeks delokalizacji para (ang. para-delocalization 
index) 

indeks donorowo-akceptorowy elektronów 7 (ang. 
melectron donor-acceptor index) 

energia rezonansu (ang. resonance energy) 

energia rezonansu w przeliczeniu na elektron 7 (ang. 
resonance energy per pi-electron) 

efekt podstawnikowy (ang. substituent effect) 

energia stabilizacji wpływu podstawnika (ang. substi- 
tuent effect stabilization energy) 
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WPROWADZENIE 


Pojęcie aromatyczności (aromaticity/aromatic) należy do najczęściej używa- 
nych w chemii i naukach pokrewnych - średnio ukazuje się 35 publikacji dzien- 
nie na ten temat [1]. Efekt podstawnikowy (substituent/substituent effect) jest też 
popularnym terminem, ukazuje się średnio ok. 12 publikacji dziennie. Mimo takiej 
popularności obu terminów, prac łączących te dwa zagadnienia jest znacznie mniej: 
średnio dziennie ~2. Obie koncepcje, aromatyczność i efekt podstawkowy, mają 
długą (ponad stuletnią) historię, natomiast ich ilościowy opis jest trochę młodszy. 
Oba pojęcia (zjawiska) nieustająco fascynują i inspirują do szukania nowych cha- 
rakterystyk, które umożliwią zarówno opis i porównanie właściwości fizycznych, 
chemicznych i biochemicznych związków jak i projektowanie nowych związków, 
wykorzystywanych w chemii, biochemii i dziedzinach pokrewnych. 

Po raz pierwszy chemiczna idea aromatyczności pojawiła się jako koncepcja 
strukturalna. Kekule przypisał ten termin związkom zawierającym pierścień benze- 
nowy [2], natomiast rok później Erlenmayer [3] jako aromatyczne wskazał związki 
o właściwościach podobnych do pochodnych benzenu. 

Powszechnie wiadomo, że benzen jest substancją trującą i kancerogenną. Nato- 
miast jego pochodna — kwas benzoesowy, w postaci soli sodowej lub wapniowej są 
uznanymi konserwantami, zaś kwas acetylosalicylowy to powszechnie stosowany 
lek, znany jako aspiryna [4]. Pierwsze ilościowe ujęcie opisu efektu podstawniko- 
wego zaproponował Louis Plack Hammett [5], na podstawie wyników badań pod- 
stawionych pochodnych kwasów benzoesowych. Jego idee i zastosowania są do dziś 
wykorzystywane w codziennej pracy naukowej chemików. 

Jak widać, aromatyczność i efekt podstawnikowy łączą obiekty badań, zwłaszcza 
gdy pierścień jest transmiterem efektu podstawnika. Ważnym jest zatem spojrzenie 
na ich wzajemne relacje, obserwowane w przypadku podstawionych pochodnych 
benzenu oraz innych układów jednopierścieniowych. Przed ich zaprezentowaniem, 
zostaną przedstawione najczęściej wykorzystywane podejścia do ilościowego opisu 
obu pojęć. 


1. INDEKSY AROMATYCZNOŚCI 


W tradycyjnym rozumieniu [6], jako związki aromatyczne przyjęto określać, 
w latach sześćdziesiątych, płaskie cykliczne układy z-elektronowe wykazujące 
następujące właściwości [7, 8]: 
(i) są one bardziej trwałe od ich acyklicznych analogów [9, 10], 
(ii) mają długości wiązań pośrednie między długościami wiązań pojedyn- 
czych i podwójnych, 
(iii) zewnętrzne pole magnetyczne indukuje z-elektronowy prąd kołowy, 
powodujący wzrost podatności diamagnetycznej [7, 11], 
(iv) chętniej ulegają reakcji substytucji niż addycji [12]. 
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Powyższe kryteria mają swoje ilościowe deskryptory, zwane czasem indeksami 
aromatyczności. 

ad (i) Koncepcja energii rezonansu (ang. resonance energy, RE) [13, 14] jako 
różnicy energii zdelokalizowanego układu m-elektronowego i jego zlokalizowa- 
nego analogu pozwala oszacować ilościowo aromatyczność układu. Energie mogą 
pochodzić zarówno z pomiarów kalorymetrycznych [15], jak i można je wyznaczać 
metodami chemii kwantowej [16-18], tzw. DRE (ang. Dewar resonance energy) [19]. 
Zastosowanie HMO (ang. Hiickel Molecular Orbital theory) pozwoliło wyznaczyć 
RE dla ogromnej liczby karbo- i hetero- cyklicznych cząsteczek m-elektronowych. 
Przeliczone wartości RE na elektron 7, REPE (ang. Resonance energy per n-electron) 
[20-23] stały się wielkościami intensywnymi, pozwalającymi na porównanie właści- 
wości aromatycznych układów o różnych rozmiarach [24]. Im większa wartość 
REPE tym układ jest bardziej aromatyczny. Aktualnie termin RE zastępuje się często 
pojęciem ASE (ang. Aromatic Stabilization Energy) [25], a liczba modyfikacji ASE 
jest bardzo duża, zależy nie tylko od metody obliczeniowej ale też od wybranego 
układu odniesienia. W artykule przeglądowym [25] w przypadku benzenu przed- 
stawionych jest aż 15 takich układów. Interesującym podejściem energetycznym jest 
koncepcja HOSE (ang. Harmonic Oscillator Stabilization Energy), pozwalająca nie 
tylko na ocenę energetyczną układu m-elektronowego [20, 22] ale także na oszaco- 
wanie udziału struktur rezonansowych go opisujących [26], z bardzo dobrą zgodno- 
ścią z innymi tego typu opracowaniami [27, 28]. 

ad (ii) Geometryczne indeksy aromatyczności. Ilościowe wykorzystanie wyrów- 
nania długości wiązań jako deskryptora aromatyczności A, zaproponowali Julg 
i Francois [29], w postaci znormalizowanej funkcji wariancji długości wiązań peryfe- 
ryjnych, tak aby dla benzenu A, = 1 a dla butadienu A, = 0. Rok później został zapro- 
ponowany podobny indeks, oparty na rzędach wiązania [30], obliczanych metodą 
HMO. Kolejnym krokiem była koncepcja odniesienia wartości długości wiązań, d, 
nie do średniej długości wiązań ale do empirycznej wielkości nazwanej optymalną 
długością wiązania, d. tj. taką jaką dane wiązanie (są tu rozważane także układy 
z heteroatomami) miało by w układzie idealnie aromatycznym [31, 32]. Indeks ten 
nazwano HOMA (ang. Harmonic Oscillator Model of Aromaticity), w przypadku 
węglowodorów ma postać (1): 


HOMA =1--)a/d 


H 


opt —d, y (1) 
gdzie n jest liczbą wiązań CC, a = 257,7 jest stałą empiryczną taką aby HOMA = 0 
dla układów niearomatycznych i HOMA = 1 dla układów, w których Św = 1,388 A, 
zaś d, są długościami wszystkich wiązań w danym układzie. 

Wartości parametrów empirycznych dla układów z heteroatomami są podane 
w publikacjach [32-37]. Bardzo ważną zaletą indeksu HOMA jest to, że można go 
stosować do dowolnego fragmentu 7-elektronowego cząsteczki, znajdując w ten 
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sposób informację o delokalizacji jego elektronów m. Indeks HOMA ma kilka 
modyfikacji [38-40], które jednak nie zmieniają jego zasadniczej idei. 

ad (iii) Magnetyczne indeksy aromatyczności. Początkowo aromatyczność 
łączono z takimi właściwościami magnetycznymi jak wzbudzany zewnętrznym 
polem magnetycznym prąd pierścieniowy (ang. ring current) [41] i związana z nim 
egzaltacja podatności magnetycznej [42-44]. Warto zauważyć, że ta ostatnia wiel- 
kość opisuje cząsteczkę jako całość. Ważnym przełomem było zaproponowanie 
przez Schleyera [45] ilościowej miary aromatyczności, w postaci ujemnej wartości 
absolutnej wielkości przesłaniania obliczanej w geometrycznym środku cząsteczki 
(pierścienia), nazwanej NICS (ang. Nucleus Independent Chemical Shift). NICS 
zależy explicite od wartości prądu kołowego pierścienia i może służyć do oceny delo- 
kalizacji jego elektronów m. Aktualnie NICS jest obliczany także w innych miejscach 
cząsteczki, lub określany jako jego składowa tensora prostopadła do płaszczyzny 
cząsteczki, tj. NICS(1)zz [46, 47]. Pewną wadą NICS jest to, że jego wartości zależą 
od wielkości pierścienia. 

ad (iv) Jedynym ilościowym, a równocześnie ogólnym, ujęciem charakteru aro- 
matycznego cząsteczek pi-elektronowych opartym na braku ich skłonności do ule- 
gania reakcji addycji, gdy możliwa jest reakcja konkurencyjna czyli substytucja, jest 
indeks KK [48]. Jest on zdefiniowany jako ilość energii z-elektronowej (obliczana 
metodą HMO) jaką cząsteczka traci w wyniku reakcji addycji w pozycji r i s, gdy 
w tych pozycjach następuje zmiana hybrydyzacji z sp” na sp”. Im większa jest war- 
tość KK, tj. większa strata energii powodowanej przez reakcję addycji, tym trudniej 
zachodzi ta reakcja i tym bardziej aromatyczna jest cząsteczka. Dzięki zastosowaniu 
indeksu KK możliwa była klasyfikacja układów zr-elektronowych na aromatyczne 
(wg reguły Hiickela zawierające 4N+2 elektrony 7) oraz anty-aromatyczne (zawiera- 
jące 4N elektronów 7), co w przypadku annulenów dokonane było przy okazji zasto- 
sowania różnie zdefiniowanych energii rezonansu [49, 50]. Dodatkowo, indeks KK 
pozwolił do takiej klasyfikacji włączyć inne układy 7-elektronowe, jak pokazuje to 
Rys. 1 [51]. Przedstawione krzywe opisują anuleny 4N+2 (górna, przerywana) i 4N 
(dolna, ciągła), pozostałe cząsteczki (również nie pierścieniowe) mają oznaczenia 
indywidualne. Uogólnienie w/w podejścia przedstawili Krygowski i Kruszewski [52]. 

Oprócz przedstawionych wyżej charakterystyk aromatyczności, są jeszcze 
podejścia oparte na właściwościach elektronowych cząsteczek, w szczególności na 
wartości ładunków z-elektronowych na atomach. Już w 1970 r. zaproponowano 
indeks aromatyczności D_ zdefiniowany jako suma wartości bezwzględnych różnic 
między wartością średnią ładunku 7-elektronowego (wyznaczanego metodą HMO) 
a poszczególnymi ładunkami na atomach węgla w cząsteczce [53]. Im większa jest 
wartość D, tym mniejsza aromatyczność układu. Bardziej nowoczesne są podejścia 
oceniające fluktuacje elektronów między sąsiadującymi atomami w układzie, FLU 
[54], albo indeks delokalizacji elektronów w pozycji para, PDI [55]. Więcej infor- 
macji o nowoczesnych elektronowych indeksach aromatyczności można znaleźć 
w pracy przeglądowej [56]. 
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Rysunek 1. Zależność KK od liczby elektronów pi cząsteczki [51] 
Figure 1. Dependence of KK index on the number of pi-electrons in molecules. Reprinted (adapted) from 
Tetrahedron Lett., 1970, 11, 1311 [51]. Copyright (1970), with permission from Elsevier 


W świetle takiej różnorodności metod ilościowego opisu aromatyczności, 
warto postawić pytanie; w jakim stopniu są one równoważne? Problem ten podjęło 
wielu autorów [57-62]. Odpowiedzią jest stwierdzenie, że w sensie przybliżonym 
są one porównywalne [63] i to tym bardziej im w większym stopniu rozpatrywane 
układy są do siebie podobne strukturalnie. 


2. CHARAKTERYSTYKI EFEKTU PODSTAWNIKOWEGO 


Od dawna wiadomo, że wprowadzenie podstawnika i jego zmiany mogą 
w decydujący sposób zmienić właściwości chemiczne, fizyczne a nawet bioche- 
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miczne cząsteczek chemicznych. Hammett, jako pierwszy, zaproponował ilościowy 
opis efektu podstawnikowego [5, 64], wprowadzając pojęcie stałej podstawnikowej 
o zdefiniowanej jako 


o = lg K(X) - Ig(H) (2) 


gdzie K(X) i K(H) są stałymi równowagi dysocjacji kwasowo-zasadowej podstawio- 
nego i niepodstawionego kwasu benzoesowego w warunkach normalnych w wodzie. 

Podstawowym założeniem koncepcji Hammetta było, że zmiany właściwości 
chemicznych (co potem rozszerzono na inne właściwości, P(X)), innych podsta- 
wionych układów będą podobne do tych obserwowanych w kwasach benzoesowych 
i zatem będą skorelowane ze stałymi podstawnikowymi o(X) w postaci równania 
(3), zwanego równaniem Hammetta: 


P(X) = po(X) + const (3) 


gdzie p jest tzw. stałą reakcji i zależy od rodzaju reakcji/procesu oraz warunków jej 
przeprowadzenia i opisuje czułość reakcji/procesu na efekt podstawnikowy. 
Ostatnio, Sadlej-Sosnowska [65, 66] zaproponowała metodę oceny efektu pod- 
stawnikowego w oparciu o zastosowanie metod obliczeniowych chemii kwantowej, 
wprowadzając pojęcie „charge of the substituent active region” , cSAR(X), zdefinio- 
wanego jako suma ładunków atomowych podstawnika X oraz ładunku na podsta- 
wionym atomie węgla, je cSAR(X) = q(X) + a Wartości cSAR(X) dobrze 
korelują ze stałymi podstawnikowymi [65, 66], a w dwu-podstawionych pochod- 
nych X-Ph-Y (Y - miejsce reakcji) pozwalają oszacować przeniesienie ładunku (ang. 
charge transfer) z X do Y [67]. W odróżnieniu od ładunku na podstawniku, q(X), 
niezależnie od metody obliczania ładunków atomowych wartości cSAR(X) dobrze 
korelują ze sobą [68]. Wiąże się to z faktem, że wiązania C-X mogą być bardzo 
polarne, podczas gdy wiązania CC są wolne od tej cechy, co pokazano na Schema- 
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Schemat 1. Graficzna prezentacja definicji q(X) (a) i cSAR(X) (b) 
Scheme 1. Graphical presentation of q(X) (a) and cSAR(X) (b) definitions 


Użytecznym sposobem opisu efektu podstawnikowego jest zaproponowana 
wcześniej reakcja homodesmotyczna [69, 70] 
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X-R-Y + R = X-R + R-Y 


Energia tej reakcji, nazwana SESE (ang. Substituent Effect Stabilization Energy), 
jest energetycznym deskryptorem efektu podstawnikowego, obliczanym zgodnie 
z równaniem (4) 


SESE = E(X-R) + (E(R-Y) - [E(X-R-Y) + E(R)] (4) 


Wartości SESE są ona na ogół dobrze skorelowane ze stałymi podstawnikowymi 
[71] Dodatnie wartości SESE (SESE > 0) oznaczają, że oddziaływania wewnatrzcza- 
steczkowe X--:Y mają charakter stabilizujący cząsteczkę X-R-Y. 

Szersze omówienie kwantowo-chemicznych modeli efektu podstawnikowego 
przedstawiliśmy dwa lata temu [4]. 


3. UKŁADY JEDNOPIERŚCIENIOWE 


3.1. POCHODNE BENZENU I INNYCH PIERŚCIENI WĘGLOWYCH 


Pierwszą ilościową próbą określenia zależności aromatyczności cząsteczki od 
efektu podstawnikowego było zastosowanie indeksu D, do 10 monopodstawionych 
pochodnych benzenu [53]. Równanie (5) przedstawia równanie regresji liniowej D, 
względem wartości bezwzględnej stałych podstawnikowych Hammetta o. 


D, = -0,915 |o,| + 0,084 (5) 


Współczynnik korelacji wynosi —0,946. 

Powyższa zależność oznacza, że im silniej oddziałuje podstawnik z pierścieniem 
benzenowym, tym większe są zmiany ładunków n-elektronowych na atomach węgla 
pierścienia i tym większe są wartości D, i tym mniejsza aromatyczność pierścienia. 

Zastosowanie bardziej nowoczesnych indeksów aromatyczności do takich 
samych układów miało miejsce ponad trzydzieści lat później [72], wykorzystano 
następujące indeksy: ASE, HOMA [32], NICS% [47] i PDI [55]. W badaniach uwz- 
ględniono 19 podstawników, a do ich opisu użyto stałych o. W rezultacie, z wyjąt- 
kiem ASE, nie potwierdziły się wyniki równania (5); wszystkie pozostałe indeksy 
wskazały na dużą odporność pierścienia benzenowego na zaburzenie jego struktury 
n-elektronowej przez wprowadzenie podstawnika. Jest to w zgodzie ze znaną regułą 
chemii organicznej, że struktura elektronowa układów aromatycznych jest odporna 
na działanie czynników zewnętrznych (substytucja vs addycja) [12]. 

Naturalną konsekwencją powyższej konkluzji było sprawdzenie wpływu pod- 
stawników na układy n-elektronowe o niskiej aromatyczności. Klasycznymi przy- 
kładami takich układów są pentafulwen i heptafulwen (Schemat 2). Należy spodzie- 
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wać się, że elektrono-donorowy podstawnik podnosi aromatyczność pentafulwenu, 
natomiast elektrono-akceptorowy podnosi ją w przypadku heptafulwenu. W obu 
przypadkach liczba elektronów 7 w pierścieniach zbliża się do 4N+2 = 6. 
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Schemat 2. Pentafulwen (a) i heptafulwen (b) 
Scheme 2.  Pentafulvene (a) and heptafulvene (b) 
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Wpływ podstawników na delokalizację elektronów 7 w pierścieniu (tj. aroma- 
tyczność) w egzocyklicznie podstawionych pochodnych fulwenu opracowano wyko- 
rzystując dane doświadczalne geometrii cząsteczek [73], zastosowane do wyzna- 
czenia indeksu HOMA [74]. Znaleziono bardzo duży zakres zmienności wartości 
HOMA, od -0,106 dla 6-(4-dimetylo-aminofenylo)fulwenu do 0,702 dla 6-dime- 
tyloaminopiperidyno-fulwenu. Natomiast największą wartość, HOMA = 0,986, 
wyznaczono w przypadku soli wapniowej di-cyclopentadienylowej, gdzie pierścień 
pięcioczłonowy jest częścią anionu i ma już 6 elektronów 7, a zatem spełnia regułę 
4N+2 Hiickela. Podobne wyniki uzyskano dla kompleksu pentafulwenu z litem [75], 
a następnie jego kompleksów z innymi metalami alkalicznymi [76]: Li, Na, K, Rb 
oraz Cs, stosując różnorodne indeksy aromatyczności (NICS, HOMA, pEDA). W 
przypadku wolnego pentafulwenu HOMA = -0,297, natomiast w badanych ukła- 
dach HOMA ~ 0,56, wartości te (zmiany aromatyczności) są zgodnie ze wskaza- 
niami pozostałych indeksów aromatyczności. 

Badano również 29 egzocyklicznie podstawionych pochodnych pentafulwenu 
(Schemat 2a), analizując wpływ podstawnika na delokalizację 2-elektronowa w pier- 
ścieniu [77], do opisu której wykorzystano prąd pierścieniowy, HOMA i pEDA 
(suma obsadzeń orbitali 7-elektronowych pierścienia) [78]. Stwierdzono bardzo 
dobrą korelację pomiędzy HOMA i pEDA (R° = 0,988), zgodnie z poniższym rów- 
naniem: 


HOMA = 1,502 pEDA - 0,239 (6) 


W zależności od elektrono-donorowych/akceptorowych zdolności podstawni- 
ków, wartości HOMA zmieniają się w bardzo dużym zakresie (od ~-0,5 do ~0,7). 

Hipotezę, że elektrono-akceptorowy podstawnik powinien spowodować wzrost 
aromatyczności pierścienia heptafulwenu potwierdziły wyniki badań jego kom- 
pleksów z atomami chlorowców (F, Cl, Br, I i At) [79]. Do opisu aromatyczności 
wykorzystano indeksy: HOMA, pEDA i NICS. W pierścieniu heptafulweniowym 
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znajduje sie 7 mr-elektronów, a ich oddziaływanie z atomami chlorowca prowadzi 
do przeniesienia ładunku na ten atom i w konsekwencji powstaje anion chlorow- 
cowy. Proces ten powoduje zmianę wartości HOMA z 0,165 dla wolnej cząsteczki do 
0,640 dla soli fluorowej. Obserwowana zmiana jest tym większa im bardziej elektro- 
ujemny jest chlorowiec. Co więcej, tu również stwierdzono bardzo dobrą korelację 
pomiędzy HOMA i pEDA (cc = —0,999) oraz pomiędzy energią oddziaływań i NICS 
(cc = —0,995). 

Spektakularnym przykładem wpływu podstawnika na aromatyczność pierście- 
nia benzenowego jest zależność HOMA pierścienia od podstawników w przypadku 
wielopodstawionych pochodnych fenolu [80], w których „moc wiązań wodoro- 
wych” jest charakteryzowana przez długość wiązania CO. Zależność, uzyskaną na 
podstawie wyników badań krystalograficznych [73], przedstawia Rysunek 2. 

Symulacja podobnego typu oddziaływań metodami chemii kwantowej, wyko- 
rzystując prosty model p-podstawionych fenoli oddziaływujących z jonem fluorko- 
wym, doprowadziła do bardzo podobnego obrazu [81, 82]. Dodatkowo, udokumen- 
towano dobrą korelację pomiędzy indeksami HOMA i NICS, wykorzystanymi do 
opisu aromatyczności pierścienia. 
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HOMA 
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Cé 
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ry 
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Rysunek 2. Zależność HOMA od długości wiązania C—O, d. dla różnie podstawionych pochodnych fenolu 
[80] 

Figure 2. Relationship of HOMA on C-O bond length, d.,, for variously substituted phenols interac- 
ting with bases. Reprinted with permission from J. Chem. Inf. Comput. Sci., 2004, 44, 2077 [80]. 
Copyright (2004) American Chemical Society 


Przykładami słabo- albo nawet niearomatycznych układów z-elektronowych 
są chinony, a dwie grupy karbonylowe mają dwie różne ścieżki (drogi) oddziały- 
wania z podstawnikiem. Zastosowanie indeksów: HOMA, MCI, DI oraz FLU, do 
jedenastu 4-podstawionych pochodnych 1,2-benzo-chinonu (Schemat 3) pokazało 
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różnicę w delokalizacji w zależności od drogi oddziaływania [83]. Bardzo wyraźnie 
widać to na wykresie zależności długości wiązania CO od stałej podstawnikowej 
o, (o, dla podstawników elektrono- donorowych), przedstawionym na Rysunku 3. 
Im bardziej elektrono-donorowy podstawnik, tym większy charge transfer do grupy 
karbonylowej i tym dłuższe wiązanie, co więcej — jedna z grup karbonylowych jest 
znacząco bardziej wrażliwa na efekt podstawnikowy. 


O 


X 


Schemat 3. 4-X-1,2-benzochinony, X = NO, NO,, CN, CHO, H, Me, OMe, OH, NH,, NHMe i NMe, 
Scheme 3. 4-X-1,2-benzoquinones, X = NO, NO,, CN, CHO, H, Me, OMe, OH, NH,, NHMe and NMe, 


doo = -0.0046 o + 1.212 


dc=o = -0.0005 o + 1.210 


-2.0 10 6 0.0 1.0 


Rysunek 3. Zależność d. — od stałej podstawnikowej o (o, dla podstawników elektrono-donorowych, dla 


pozostałych 0,) [83] 

Figure 3. Dependence of CO bond lengths, d._„, on substituent constants o (o; for electron-donating sub- 
stituents and for others o). Reprinted with permission from J. Org. Chem., 2011, 76, 550 [83]. 
Copyright (2011) American Chemical Society 


Pdobne zależności otrzymano w przypadku relacji pomiędzy HOMA oraz MCI 
a stałymi podstawnikowymi. Z uzyskanych równań liniowych wynika, że większe 
są zmiany nachylenia dla oddziaływań podstawnika X w pozycji 4 z grupą karbo- 
nylowa C20 niż z CIO, czyli silniejsze są oddziaływania na drodze OC2C3C4 niż 
OC1C6C5C4. Wyniki te pozwalają na sformułowanie prostej reguły [83]: jeżeli 
między atomem elektrono-akceptorowym i elektrono-donorowym liczba wiązań 
jest parzysta, wówczas uprzywilejowany jest efekt rezonansowy, tj. przeniesienie 
ładunku. 
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Przeprowadzone zostały również systematyczne badania wpływu podstawnika 
na delokalizacje 2-elektronowa i właściwości transmisyjne 3-X i 4-X di-podsta- 
wionych pochodnych benzenu (BEN) i cykloheksa-1,3-dienu (CHD), X = NMe,, 
NH, OH, OMe, CH, H, F, Cl, CF, CN, CHO, COMe, CONH,, COOH, NO,, NO. 
Właściwości podstawników charakteryzowano wykorzystując zarówno deskryptory 
oparte na chemii kwantowej: cSAR i SESE, jak również tradycyjne stałe podstaw- 
ników Hammetta, natomiast delokalizację mr-elektronową oszacowano przy użyciu 
indeksu HOMA. Wyniki dla 4-X podstawionych 1-nitro i 1-hydroksy pochodnych 
CHD [84] przedstawiono na Rysunku 4. W obu seriach (Y = NO, i OH) stwier- 
dzono, że gdy właściwości elektronowe podstawników są przeciwne do właściwości 
ustalonej grupy (Y), wartości HOMA dobrze korelują z SESE, ze znaczną warto- 
ścią współczynnika kierunkowego. Efekt ten nie jest obserwowany, jeśli X i Y mają 
podobne właściwości, czyli oba są elektrono-donorowe lub elektrono-akceptorowe. 


= Y=NO; 
=OH 
y = 0.125x+ 0.236 e Hei 
R=0899 e 
GË 
y = -0.016x + 0.280 0.6 PE 
$ "w" (Zi ©0<0 
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= ve 297 E  v=0.083x+0.305 
"HM O.O... OOS, "o R? = 0.982 0620 
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S e? Z 8.0<0 
y=-0.019x+0.135 @ 38 
R? =0.086 ri 
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Rysunek 4. Zależności HOMA od SESE dla 4-X-CHD-NO, i 4-X-CHD-OH pochodnych 

Figure 4. Dependences of HOMA on SESE for 4-X-CHD-NO2 and 4-X-CHD-OH series. Reused from [84], 
this work is licensed under the Creative Commons Attribution 4.0 International License (http:// 
creativecommons.org/licenses/by/4.0/) 


Koncepcja cSAR umożliwia opis właściwości transmisyjnych podstawnika 
z danej pozycji do centrum reakcji. Wskaźnik przepływu ładunku (ang. charge flow 
index, CFI) zdefiniowano jako: 


CFI = cSAR(Y) - cSAR(X) (7) 


Wyznaczone zależności liniowe CFI, , zap VS CFL 4 xcup.y Oraz CFL a przy 
vs CFI, , weu Umożliwiły porównanie transmisji efektu podstawnikowego w obu 
seriach, w każdej z nich uwzględniono sześć ustalonych grup (Y = NO,, COOH, 
COO, OH, O i NH,). Wartości współczynników kierunkowych wskazują na 
słabszą transmisję z pozycji meta (3-) niż z pozycji para (4-). Co więcej, efekt ten 
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(słabsza transmisja) jest znacznie silniejszy w przypadku pochodnych cykloheksa- 
-1,3-dienu niż benzenu, średnie wartości współczynników kierunkowych wynoszą 
odpowiednio 0,532 i 0,735. Należy dodać, że współczynniki determinacji wyznaczo- 
nych korelacji były większe niż 0,92. 


3.2. UKŁADY HETEROCYKLICZNE 


Aromatyczność pięcioczłonowych pierścieni zawierających atom(y) azotu, 
nazywanych azolami, także zależy od obecności podstawników. Tetrazol (Sche- 
mat 4) występuje w dwóch formach tautomerycznych (1H i 2H) i podobnie jak ben- 
zen ma 6 elektronów 7. Badania C5-podsatwionych pochodnych (X = AIH,, BeH, 
DH. CCH, CF,, CHO, Cl, CN, F, Me, NH,, NO, NO, OH, SH, SiH,) pokazały [85], 
że zmiana struktury 7-elektronowej, opisana przez pEDA, zależy wyraźnie od sta- 
łych podstawnikowych: 


1H tautomer pEDA = -0,146 0, - 0,034 cc = —0,905 
2H tautomer pEDA = —0,141 0, - 0,031 cc = —0,924 
2 3 H 
N=N NN 
ee i 
X X 
a b 


Schemat 4. C5-X podstawione pochodne 1H (a) i 2H (b) tautomerów tetrazolu 
Scheme 4. C5-X substituted 1H (a) and 2H (b) tetrazoles 


Ponadto, obsadzenie orbitali 2pz wszystkich atomów azotu, z wyjątkiem N3, 
dobrze koreluje z wartościami pEDA (cc > 0,95). Brak korelacji w przypadku atomu 
N3 może sugerować, że ta pozycja w 1H i 2H tetrazolach przypomina do pewnego 
stopnia pozycję meta w benzenie. Potwierdziły to analogiczne badania mono-pod- 
stawionych pochodnych benzenu. Analogiczne zależności obsadzenia orbitali 2pz 
atomów węgla w pozycji o- i p- od pEDA mają wysoki poziom korelacji (odpowied- 
nio cc = 0,971 i 0,968) i dodatnie wartości współczynników nachylenia, natomiast 
w przypadku pozycji m- cc = -0,781, z nachyleniem ujemnym. Ten efekt jest jakby 
dodatkowym obrazem przyczyn dla których wpływ podstawnika w o- i p- pochod- 
nych benzenu tak się różni od takich oddziaływań dla m-pochodnych. 

Różnice między wpływem podstawnika na aromatyczność w mono-podstawio- 
nych pochodnych benzenu i pirazolu oraz imidazolu, w każdym przypadku pier- 
ścień zawiera 6 elektronów 7, jest doskonale pokazana na Rysunku 5 [86]. W prze- 
ciwieństwie do powyżej omawianych tetrazoli, w badanych pochodnych pirazolu 
i imidazolu podstawnik jest dołączony do atomu azotu. Do opisu aromatyczności 
pierścieni wykorzystano indeks HOMA. W przypadku podstawników elektrono- 
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-donorowych (o < 0) wartości HOMA są większe niż 0,8, natomiast dla podstaw- 
ników elektrono-akceptorowych — mniejsze niż 0,8. W tym ostatnim przypadku 
podstawniki wyciągają elektrony 7 z pierścienia, prowadząc do utworzenia układów 
nie spełniających zasady 4N + 2. 


HOMA(Bz) = 0.0070} + 0.988 


Pie 1.000 
AR a 0.950 
UN 
N 
| 0.900 
Pz X e 
0.850 Š 
= 
N 0.800 
W \ HOMA(PZ) = -0.1620 + 0.836 
r=-0.87 
Im y HOMA(Im) = -0.166c, + 0.776 9.220 
x = -0.91 
0.700 


-0.80 -0.60 -0.40 -0.20 0.00 0.20 0.40 


Rysunek 5. Korelacde pomiędzy indeksem HOMA a rezonansową stałą podstawnika (o,) dla pochodnych 
benzenu (Bz), pirazolu (Pz) i imidazolu (Im) [86] 

Figure 5. Correlations between HOMA aromaticity index and resonance substituent constant (0,) for sub- 
stituted benzene (Bz), pyrazole (Pz) and imidazole (Im) derivatives. Adapted with permission 
from J. Phys. Chem. A, 2011, 115, 8571 [86]. Copyright 2011 American Chemical Society 


PODSUMOWANIE 


Zjawiska znane jako aromatyczność i efekt podstawnika są jednym z najważ- 
niejszych zagadnień w chemii, biochemii i dziedzinach pokrewnych. Do ilościo- 
wego opisu aromatyczności najczęściej używane są indeksy oparte na energii, 
geometrii (na przykład HOMA), właściwościach magnetycznych (na przykład 
NICS) i strukturze elektronowej (np. FLU). Efekt podstawnikowy najczęściej jest 
charakteryzowany przez stałe podstawnika (np. stale Hammetta). Jednakże, rozwój 
metod obliczeniowych umożliwił zastosowanie deskryptorów efektu podstawnika 
wyznaczanych w oparciu o metody chemii kwantowej. W tym celu wykorzystywane 
są energie odpowiednio zaprojektowanych reakcji homodesmotycznych, rozkład 
gęstości elektronowej lub potencjał elektrostatyczny. Oczywiście, możliwość ich 
zastosowania jest weryfikowana przez porównanie „nowych” deskryptorów z otrzy- 
manymi przy użyciu opisu klasycznego (czyli stałe typu Hammetta). Wśród nowych 
koncepcji fizykochemicznych opisu efektu podstawnikowego najbardziej obiecują- 


O ZALEŻNOŚCIACH POMIĘDZY AROMATYCZNOSCIA I EFEKTEM PODSTAWNIKOWYM 259 


cym podejściem jest cSAR, który pozwala badać zarówno klasyczne jak i odwrotne 
efekty podstawników. 


Wzajemne relacje między aromatycznością a efektem podstawnikowym można 


podsumować następująco: 


(i) Im mniej aromatyczny jest układ, tym bardziej wrażliwy na efekt pod- 
stawnika. Silnie aromatyczne cząsteczki są odporne na działanie podstaw- 
ników. 

(ii) Wewnątrzcząsteczkowego przeniesienie ładunku (efekt rezonansowy) jest 
znacznie bardziej efektywne w przypadku parzystej liczby wiązań pomię- 
dzy atomami elektrono-akceptorowymi i elektrono-donorowymi grup 
w układzie 7-elektronowym. 

(iii) Wewnątrzcząsteczkowe oddziaływania — transmisja efektu podstawnika 
— są znacznie słabsze w układach podstawionych w pozycjach 1,3 (czyli 
typu meta) niż w pozycjach 1,4 (typu para). 
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ABSTRACT 


General anesthesia is defined as impairment of the central nervous system 
(UON) caused by intravenous or volatile anesthetics. The state of loss of conscio- 
usness or even amnesia and the disappearance of perception into external stimuli 
is achieved by the use of a large group of chemical compounds. The use of nitrous 
oxide in 1844 revolutionized surgery and medicine at that time. From that moment, 
anesthesiology develops dynamically, allowing more and more complex procedures. 
Despite more than 170 years of history of anesthesia, understanding the mecha- 
nism of reversible loss of awareness and sensitivity to pain caused by the action of 
general anesthetics is one of the greatest challenges of modern pharmacology and 
neuroscience. Incredibly high diversity of anesthetics, including both noble gases 
and complex steroids, combined with human sensation makes the above problem 
extremely difficult to solve. 

The reversibility of the anesthesia phenomenon suggests that the analyzed 
phenomenon is based on disturbance of weak intermolecular interactions, such as 
hydrogen bond or van der Walls forces. Anesthetic molecules may bind directly 
to the hydrophobic region of protein, which causes its conformational changes or 
disturb ion channel activity by anesthetic-induced perturbations of lipid bilayers. 
The mechanism of anesthesia is thus very often attributed to both protein and lipid 
membrane targets. The influence of anesthetic molecules on biomolecular systems 
can be studied successfully using many different physico-chemical methods, such 
as, infrared, fluorescence or nuclear magnetic resonance spectroscopy. Vibrational 
circular dichroism as well as differential scanning calorimetry can also be used. 


Keywords: inhalation anesthetics, anesthesia, weak intermolecular interactions, 
lipid membrane 

Słowa kluczowe: anestetyki wziewne, anestezja, słabe oddziaływania międzyczą- 
steczkowe, błona lipidowa 
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WYKAZ STOSOWANYCH SKROTOW 


P-helisa zawierająca 6,3 reszt aminokwasowych na 
skret 

dipalmitoilofosfatydylocholina (ang. dipalmitoilpho- 
sphatidylcholine) 

różnicowa kalorymetria skaningowa (ang. differential 
scanning calorimetry) 

kwas gamma-aminomasłowy (ang. gamma aminobu- 
tyric acid) 

lamelarna faza ciekłokrystaliczna (ang. lamellar liquid 
crystalline) 

lamelarna faza zelowa (ang. lamellar gel phase) 
lamelarna faza przemienna (ang. lamellar interdigita- 
ted phase) 

lamelarna faza krystaliczna (ang. lamellar crystalline 
phase) 

minimalne stezenie pecherzykowe (ang. minimum 
alveolar concentration) 

magnetyczny rezonans jadrowy (ang. nuclear magne- 
tic resonance) 

pofalowana faza zelowa (ang. rippled gel phase) 
temperatura głównego przejścia fazowego (ang. 
the temperature of the main phase transition) 
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WPROWADZENIE 


Znieczulenie ogólne jest definiowane jako upośledzenie ośrodkowego układu 
nerwowego (UON) wywołane przez dożylne bądź lotne anestetyki [1, 2]. Stan utraty 
świadomości czy nawet amnezji oraz zaniku percepcji na zewnętrzne bodźce osią- 
gany jest przez zastosowanie licznej grupy związków chemicznych. Wykorzysta- 
nie w 1844 roku, po raz pierwszy w historii, podtlenku azotu zrewolucjonizowało 
ówczesną chirurgię i medycynę. Od tego momentu anestezjologia rozwija się pręż- 
nie, umożliwiając przeprowadzanie coraz bardziej złożonych zabiegów [3]. Mimo 
ponad 170-letniej historii, zrozumienie mechanizmu odwracalnej utraty świado- 
mości i czułości na ból, wywołanej działaniem środków znieczulenia ogólnego, jest 
jednym z największych wyzwań współczesnej farmakologii i neuronauki. Niewia- 
rygodnie duża różnorodność anestetyków, obejmująca zarówno gazy szlachetne jak 
i złożone steroidy, w połączeniu z ludzkim odczuwaniem czyni powyższy problem 
niezwykle trudny do rozwiązania. 


1. ANESTETYKI WZIEWNE 


Pierwsze środki znieczulenia ogólnego, jak eter dietylowy czy chloroform, cha- 
rakteryzowały się łatwopalnością i wysoką toksycznością (Rys. 1). W celu zwięk- 
szenia bezpieczeństwa operowanego pacjenta dążono do redukcji niepożądanych 
efektów, m.in. przez podstawienie atomów wodoru atomami fluorowców. 


cl 


cl cl 


CHLOROFORM 


Rysunek 1. Historia rozwoju anestezjologii 
Figure 1. History of development of anesthezia 


W roku 1960 wprowadzono do użytku klinicznego metoksyfluran, będący 
pochodną eterową. Mimo eliminacji wielu efektów ubocznych, wywołanie ane- 
stezji przez ten związek trwało zbyt długo. Równie trudne okazało się wybudza- 
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nie pacjenta. W niewielkim odstępie czasowym wdrożono dwa izomery: enfluran 
i izofluran, spośród których izofluran charakteryzował się większą odpornością na 
metabolizm, a tym samym nie wywierał szkodliwego wpływu na organy. Dążąc do 
uzyskania farmakologicznie perfekcyjnego środka znieczulającego, nie poprzestano 
na izofluranie. Usunięcie wszystkich halogenów, za wyjątkiem fluoru poskutkowało 
otrzymaniem niepalnej cieczy o niezwykłej wytrzymałości na metabolizm. Desflu- 
ran, eter etylowo-metylowy, w którym wszystkie wodory zostały zastąpione atomami 
fluoru, wprowadzono w 1992 roku. Dwa lata później zsyntezowano pochodną eteru 
metylowo-izopropylowego, sewofluran. Spełnienie niemal wszystkich wymagań 
stawianych idealnemu środkowi znieczulającemu pozwoliło na stosowanie zarówno 
sewofluranu jak i desfluranu we współczesnej anestezjologii [2-6]. 


1.1. WŁAŚCIWOŚCI FIZYKOCHEMICZNE WZIEWNYCH 
SUBSTANCJI ZNIECZULAJĄCYCH 


Grupa środków znieczulenia ogólnego obejmuje szereg strukturalnie różno- 
rodnych związków, które mogą być podawane na sposób wziewny lub iniekcyjny. 
Anestetyki inhalacyjne należą do trzech rozległych klas: eterów, alkanów oraz 
substancji, które w temperaturze i ciśnieniu standardowym występują w gazowym 
stanie skupienia. Nowoczesne środki znieczulające to w głównej mierze pochodne 
eterów, których cząsteczki w swej strukturze posiadają atomy fluoru oraz chloru, 
zastępujące atomy wodoru. Zabieg halogenacji prowadzi do redukcji toksyczności 
leków oraz do wzrostu oporności na metabolizm w organizmie pacjenta. Analogi 
eterów i alkanów w warunkach normalnych występują jako lotne ciecze, których 
temperatury wrzenia sięgają 50-60°C. Nowoczesne anestetyki wziewne są niepalne, 
jak również nie posiadają właściwości wybuchowych [7-11]. 


Tabela 1. Właściwości fizykochemiczne wybranych anestetyków wziewnych [2] 
Table 1. Physicochemical properties of selected inhaled anesthetics [2] 
Właściwości N,O Enfluran Izofluran Sewofluran Desfluran 
Masa czasteczkowa [g] 44 184,5 184,5 200 168 
Gęstość - 1,520 1,502 1,520 1,465 
Temperatura wrzenia [°C] —88,5 56,5 48,5 58,5 22,8 
Prężność pary 
39,8 172 240 170 669 
[mm Hg; 20°C] 
Współczynnik podziału 0,46 1,90 1,46 0,65 0,42 
krew:gaz 
MAC [%] 104,00 1,68 1,14 2,05 6,00 


Transport anestetyków wziewnych w organizmie, w odróżnieniu od leków 
podawanych dożylnie, zależy od efektywności poboru leku przez płuca, a następ- 
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nie dostarczenia do miejsca ich działania w ośrodkowym układzie nerwowym [7]. 
Wydajność dystrybucji leku jest związana z ich właściwościami fizycznymi, roz- 
puszczalnością we krwi oraz prężnością pary (Tab. 1). 

Rozpuszczalność środków znieczulenia ogólnego we krwi jest charaktery- 
zowana przez współczynnik podziału krew/gaz. To stosunek stężeń anestetyku 
wziewnego we krwi oraz w fazie lotnej. Wysoki współczynnik wiąże się z silnym 
powinowactwem leku do krwi. Przeciwnie jest w przypadku niskiej wartości tego 
parametru. Innymi słowy anestetyki o niewielkich współczynnikach charakteryzują 
się szybkim wywołaniem znieczulenia z jednoczesnym krótkim czasem wybudzenia. 
Obok desfluranu najniższą wartość stosunku krew/gaz posiada sewofluran (0,65). 
Natomiast rzadziej stosowane we współczesnej anestzjologii, izofluran i enfluran, 
posiadają współczynniki znacznie wyższe [9]. 

Siła działania anestetyków wziewnych jest powszechnie wyrażana przez mini- 
malne stężenie pęcherzykowe (ang. minimum alveolar concentration, MAC). Inter- 
pretowane jest jako stężenie anestetyku wziewnego konieczne do powstrzymania 
reakcji ruchowych u 50% populacji ludzi dorosłych w odpowiedzi na określone 
bodźce. Wysoka wartość parametru MAC wiąże się ze wzmożonym poborem gazu 
znieczulającego, co oznacza, że siła działania anestetyku jest niska. Wśród anali- 
zowanych leków najsilniejszy jest izofluran, którego minimalne stężenie pęcherzy- 
kowe wynosi 1,17% [7]. 

W warunkach klinicznych wykorzystywane są mieszaniny racemiczne leków. 
Mimo to stereoselektywność potencjału znieczulającego halogenoeterów jest 
tematem licznych rozważań. Halotan, enfluran, izofluran oraz desfluran posia- 
dają czynny optycznie chiralny atom węgla. Optyczne izomery izofluranu wyka- 
zują różnicę w oddziaływaniu z receptorami białkowymi np. z receptorem kwasu 
y-aminomasłowego (GABA). Ponieważ współczynnik bezpieczeństwa dla aneste- 
tyków wziewnych jest bardzo niski, sugeruje się że nawet stosunkowo niewielki sto- 
pień stereoselektywności może spowodować znaczącą kliniczną przewagę jednego 
z enancjomerów nad mieszaniną racemiczną [4, 12]. 


OF F Fa F 
H-0-0-0-0- F-0-0-0-0-H 
FF F FH F 
Enfluran Izofluran 


Rysunek 2. Struktura chemiczna enfluranu i izofluranu 
Figure 2. The structure of enflurne and isoflurane 


Enfluran i izofluran (Rys. 2) to izomery strukturalne, posiadające jednakową 
masę oraz objętość cząsteczkową, ok. 110 A”. Izofluran posiada niemal dwukrot- 
nie większy moment dipolowy niż enfluran. Ich wartości wynoszą odpowiednio: 
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2,0 D dla izofluranu i 0,7 D dla enfluranu. W tej rozbieżności upatruje się różnicy 
w oddziaływaniu anestetyków z białkami, a tym samym odchyleń w toksyczności 
związków [13]. 


1.2. MECHANIZM DZIAŁANIA ANESTETYKÓW WZIEWNYCH 


Zrozumienie mechanizmu działania środków znieczulenia ogólnego, prowa- 
dzącego do odwracalnej utraty świadomości i czułości na ból, jest jednym z najwięk- 
szych wyzwań farmakologii. Niezwykła różnorodność strukturalna anestetyków 
sprawia, że jest to szczególnie trudny problem do rozwiązania. 

Na poziomie komórkowym, anestetyki wpływają na funkcjonowanie komó- 
rek nerwowych, upośledzając aktywność synaps, a tym samym pracę kanałów 
jonowych. W standardowych warunkach te wyspecjalizowane białka błonowe są 
aktywowane poprzez sygnały chemiczne bądź modyfikacje otoczenia błony komór- 
kowej [14-16]. Po aktywacji, kanały zmieniają pobudliwość elektryczną neuronów, 
kontrolując przepływ przez błonę jonów pobudzających lub hamujących. Działanie 
środków znieczulenia ogólnego opiera się przede wszystkim na wzmocnieniu ruchu 
jonów hamujących bądź blokowanie sygnałów pobudzających. W odniesieniu do 
tych rozważań istnieją dwie teorie: lipidowa i białkowa [17, 18]. 


As 


stężenie pęcherzykowe (MAC) 
D s 


Minimalne 


0,01 


0,01 0,1 100 1000 


1 10 
Współczynnik podziału oliwa/gaz 


Rysunek 3. (A) Wykres zależności Mayera — Overtona [4] (B) Schemat obrazujący mechanizm znieczulenia 
oparty na zmianie ciśnienia lateralnego w błonach [29] 

Figure 3. (A) Mayer — Overton correlation [4] (B) The anesthesia mechanism based on the membrane late- 
ral pressure change [29] 


Wraz z odkryciem, że błona biologiczna jest zbudowana głównie z lipidów 
wysnuto hipotezę, że działanie anestetyków wziewnych ukierunkowane jest na 
dwuwarstwę lipidową [19]. Fundamentem tej teorii są prace Mayera i Overtona 
[20, 21] przedstawiająca korelację pomiędzy siłą działania środków znieczulających, 
a ich rozpuszczalnością w tłuszczach (Rys. 3A). Cząsteczki anestetyku wnikając do 
wnętrza błony komórkowej, w chwili przekroczenia stężenia krytycznego, zmieniają 
właściwości dwuwarstwy lipidowej. Transformacje błony towarzyszące pojawieniu 
się cząsteczek anestetyków objawiają się m.in. powiększeniem jej objętości [22], dez- 
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organizacją ułożenia lipidów lub modyfikacją profilu ciśnienia lateralnego w błonie 
[23]. Większość teorii opartych na zależności Mayera i Overtona [24, 25] postuluje, 
że zaburzenia struktury dwuwarstwy lipidowej, wywołane obecnością anestetyków 
w układzie, wpływają na funkcje białek błonowych obecnych w komórkach nerwo- 
wych. 

W 1973 roku Miller i współpracownicy [26] zaproponowali hipotezę krytycz- 
nej objętości błony lipidowej, w której to wysnuto przypuszczenia, że hydrofobowe 
cząsteczki anestetyków gromadzą się w neuronowych dwuwarstwach lipidowych 
wpływając na profil ciśnienia lateralnego w poprzek błony i powodując jej rozsze- 
rzenie. Zwiększenie objętości błony prowadzi do odwracalnych zmian w funkcjono- 
waniu kanałów jonowych, wywołując ostatecznie stan znieczulenia. 

Nieco inne podejście przedstawił Mohr [27], który sugerował, że działanie 
środków znieczulenia ogólnego związane jest z wyraźną zmianą ciśnienia lateral- 
nego w błonach. Każda błona lipidowa posiada odrębny profil rozkładu ciśnień 
lateralnych, które charakteryzują się wysokimi wartościami (Rys. 3B). Modyfikacje 
tego parametru zaburzają równowagę konformacyjną białek błonowych, wykazu- 
jących wysoką czułość na tego typu zjawiska. Wniknięcie do wnętrza błony cząste- 
czek amfifilowych, takich jak anestetyki, zwiększa selektywnie ciśnienie lateralne 
w pobliżu powierzchni dwuwarstwy. Następnie jest ono wyrównywane w kierunku 
środka błony. Gdy w odpowiedzi na impuls nerwowy kanał jonowy próbuje się 
otworzyć, powiększa się jego powierzchnia w pobliżu zetknięcia się z wodą. Jakie- 
kolwiek zmiany hamują proces otwarcia, uniemożliwiając transport neurotransmi- 
terów [28, 29]. 

Teoria lipidowa oparta na regule Mayera i Overtona [24, 25] posiada pewne 
ograniczenia. Koblin i Eger [30] udowodnili, że istnieje szeroka gama związków, 
m.in. polihalogenowe alkany, które mimo dobrej rozpuszczalności w tłuszczach nie 
wykazują sugerowanego działania znieczulającego. Ponadto przytoczone zmiany 
zachodzące w błonie lipidowej, wymagają użycia anestetyków w stężeniach znacz- 
nie przewyższających te niezbędne do wywołania znieczulenia [31]. Warto jednak 
zauważyć, że rozmieszczenie leków w błonie komórkowej nie jest jednorodne. 
W wielu przypadkach cząsteczki substancji aktywnych kumulują się w obszarach 
łączenia białka z lipidami, co skutkuje wysokim stężeniem lokalnym specyfiku. Dla- 
tego nie wyklucza się, że zastosowanie stężeń klinicznych anestetyków jest niewy- 
starczające do wywołania zmian strukturalnych dwuwarstwy lipidowej w komórce 
nerwowej [25]. 

Równocześnie z pojawieniem się pierwszych wątpliwości co do słuszności 
teorii lipidowej, zaczęto badać struktury białkowe (neurotransmitery, receptory 
czy kanały jonowe) pod kątem potencjalnych celów działania anestetyków. Franks 
i Lieb [32] zaprezentowali wyniki ukazujące zdolność środków znieczulających do 
zahamowania bioluminescencyjnych funkcji lucyferazy. Eksperymenty przeprowa- 
dzone in vitro potwierdziły wpływ środków znieczulenia ogólnego na funkcje bia- 
łek. Dlatego uważa się, że stan znieczulenia może wiązać się z zakłóceniem procesu 
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wiązania liganda przez bioczasteczke lub zaburzeniem zdolności białka do zmian 
konformacyjnych istotnych dla jego funkcjonowania [33-35]. 
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Rysunek 4. Schemat działania anestetyków na receptory GABA, i glicynowe [4] 
Figure 4. The scheme of anesthetic actions on GABA, and glycine receptors [4] 


W strukturze białkowej występuje szereg aminokwasów tworzących kieszenie 
hydrofobowe, które mogą być zajmowane przez anestetyki. Ponadto hydrofobowe 
aminokwasy są głównymi składnikami a-helis, a ich łańcuchy boczne są wyekspo- 
nowane na powierzchni białka, która styka się z błoną lipidową. Cząsteczki środków 
znieczulających, wnikając w przestrzeń łączącą białko z lipidami, mogą zakłócić ich 
oddziaływanie i bezpośrednio wpłynąć na konformację białka [2]. 

W ciągu ostatnich lat, badania nad mechanizmem działania anestetyków 
skupiły się głównie na kanałach jonowych, które znajdują się w błonie komórek 
nerwowych. W stężeniach klinicznych środki znieczulenia ogólnego wpływają na 
właściwości wydzielnicze neuronów, pozostawiając w ten sposób aksonalne prze- 
wodnictwo potencjału czynnościowego, w znacznym stopniu nienaruszone. Jednym 
z najczęściej badanych jest receptor GABA, należący do rodziny kanałów jonowych 
bramkowanych napięciem. Wraz z receptorem glicynowym stanowią grupę recep- 
torów hamujących (Rys. 4). Kwas y-aminomasłowy (GABA) pobudza receptory 
GABA,, które zwiększają przepuszczalność kanałów chlorkowych oraz wywołują 
hiperpolaryzację błony komórkowej. Większość anestetyków wziewnych, wliczając 
w to wszystkie pochodne eterów i wybrane alkany wzmacniają funkcje receptorów 
[35, 36]. 


2. WŁASNOŚCI BŁONY LIPIDOWEJ W KONTEKŚCIE 
INTERAKCJI ZE ŚRODKAMI ZNIECZULAJĄCYMI 


Błony stanowią ośrodek wielu procesów zachodzących w komórkach organi- 
zmów żywych. Pośredniczą w ustalaniu i utrzymywaniu gradientów transbłono- 
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wych, regulują transport z i do komórki, biorą udział w komunikacji międzyko- 
mórkowej oraz tworzą środowisko dla licznych reakcji biochemicznych [37-40]. 
Biorąc pod uwagę mnogość funkcji błon jest oczywiste, że wszelkie modyfikacje 
strukturalne mogą wpływać na fizjologię organizmu [41]. 

Lipidy i białka stanowią integralną część błon biologicznych otaczających 
komórki organizmów żywych. Skład lipidowo-białkowy błony oraz występujące 
pomiędzy nimi oddziaływania ściśle koreluje z jej funkcjonalnością [42]. W wielu 
przypadkach, rodzaj lipidów, które oddziałują z aminokwasami jest w pewnym stop- 
niu określony przez typ białka, z którym sąsiaduje. Białka transbłonowe wpływają 
także na właściwości strukturalne dwuwarstwy lipidowej [43, 44]. Oddziaływania 
białkowo-lipidowe są więc istotnym czynnikiem warunkującym strukturę i funkcję 
struktur komórkowych. 


2.1. STRUKTURA DWUWARSTWY LIPIDOWEJ 


Podstawą strukturalną wszystkich błon biologicznych jest dwuwarstwa lipi- 
dowa. Zbudowana z dwóch przeciwstawnych warstw cząsteczek lipidów, których 
polarne głowy skierowane są w stronę fazy wodnej, a hydrofobowe łańcuchy węglo- 
wodorowe tworzą wnętrze agregatu (Rys. 5) [39, 41]. 


Hydrofilowa głowa 


Błona lipidowa 


Hydrofobowy ogon 


Rysunek 5. Schemat budowy cząsteczki lipidu będącego elementem dwuwarstwy lipidowej [45] 
Figure 5. Structure of lipid molecule as the element of lipid bilayer [45] 
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Lipidy zorganizowane w dwuwarstwy lipidowe należą do cząsteczek amfifilo- 
wych, których hydrofilowe głowy wchodzą w skład kilku klas chemicznych [46]. 
Ponadto kwasy tłuszczowe tworzące ogon lipidowy mogą należeć do związków 
nasyconych bądź nienasyconych o różnej długości łańcucha. Wyodrębniono ponad 
100 różnych lipidów, których zawartość w komórce waha się od 20 do 80% wszyst- 
kich składników [47]. Ze względu na budowę chemiczną lipidy można podzielić na 
cztery klasy [46, 48]: 

— fosfolipidy (diacylofosfoglicerydy) 

-  glikolipidy (diacyloglikoglicerydy) 

-  sfingolipidy 

- sterole. 

Fosfolipidy dominują pod względem zawartości w błonach biologicznych, 
a wśród nich fosfatydylocholina i fosfatydyloetanoloamina (Rys. 6.A). Oprócz łań- 
cuchów węglowodorowych w swej strukturze posiadają resztę glicerolową i fosfora- 
nową, a także alkohol, definiujący ich rodzaj (Rys. 6B) [48, 49]. 


A. Skład lipidowy czerwonych krwinek 
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Rysunek 6. (A) Skład lipidowy błon ludzkich czerwonych krwinek [49] (B) Schemat fosfolipidu DPPC [50] 
Figure 6. (A) Lipid composition of human red blood cells [49] (B) DPPC structure [50] 


Obecność części hydrofobowej w cząsteczce lipidu sprawia, że większość z nich 
jest słabo rozpuszczalna w wodzie i tworzy agregaty lipidowe, które skupiają grupy 
hydrofobowe wewnątrz klasteru oraz eksponują hydrofilowe rejony na otaczający 
polarny rozpuszczalnik [51]. W zależności od środowiska oraz natury lipidów 
można wyodrębnić m.in. struktury micelarne, dwuwarstwowe tworzące liposomy. 
Dla cząsteczek o główkach polarnych zbliżonych rozmiarem do przekroju łańcu- 
chów węglowodorowych i przypominających kształtem walec, preferowane są 
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struktury dwuwarstwy lipidowej, która spontanicznie zwija się tworząc pęcherzyki 
zwane liposomami. Zamknięta powierzchnia liposomów eliminuje niekorzystną 
ekspozycję części hydrofobowych na działanie wody [49, 51, 52]. 

Błony lipidowe są układami polimorficznymi, których struktura nie zależy 
jedynie od budowy samej cząsteczki lipidu czy od jej stopnia uwodnienie, ale także 
od czynników zewnętrznych takich jak temperatura, ciśnienie, pH i siła jonowa 
[51, 53]. W warunkach fizjologicznych większość błon lipidowych występuje w fazie 
ciekłokrystalicznej (L,) bądź w fazie żelowej (L,). Unikalną cechą dwuwarstwy jest 
zdolność do przemiany fazowej będącej odwracalnym i wysoce kooperatywnym 
procesem, opierającym się na topnieniu lipidowych łańcuchów acylowych. Kon- 
wersja uporządkowanej fazy żelowej, w której sztywne łańcuchy węglowodorowe 
są bogate w konformery trans, do nieuporzadkowanej fazy ciekłokrystalicznej, cha- 
rakteryzującej się ruchliwymi łańcuchami acylowymi o konformacji typu gauche, 
nazywana jest głównym przejściem fazowym (Rys. 7). Objawia się zmniejszeniem 
grubości dwuwarstwy i jednoczesnym wzrostem całkowitej objętości zajmowanej 
przez cząsteczki lipidów. Topnienie łańcuchów węglowodorowych zwiększa obszar 
hydrofilowy dostępny dla penetrujących go cząsteczek wody [53-55]. 
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Rysunek 7. Model głównego przejścia fazowego wywołany wzrostem temperatury oraz odpowiadająca mu 
krzywa DSC [38] 

Figure 7. The main phase transition caused by an increase in temperature and corresponding DSC curve 
[38] 


Główne przejście fazowe może być wywołane zmianami temperatury otoczenia. 
Przemianę fazową osiąga się po przekroczeniu tzw. temperatury głównego przejścia, 
która jest charakterystyczna dla danego rodzaju lipidów. Dla dipalmitoilofosfatydy- 
locholiny (DPPC) przemianę fazową trans — gauche obserwuje się w temperaturze 
41°C [56]. 


276 M. HOROCHOWSKA, K. CIESLIK-BOCZULA, M. ROSPENK 


wnt AN ; 
i = Aye 


Faza krystaliczna L, Faza żelowa Lg Faza żelowa P; Faza ciekłokrystaliczna L, 


Rysunek 8. ` Przemiany fazowe zachodzące w uwodnionej dwuwarstwie DPPC [51] 
Figure 8. Phase transition processes of hydrated DPPC bilayer [51] 


Fosfatydylocholiny, które należą do grupy najczęściej występujących lipidów 
mogą ulegać dodatkowo dwóm innym przemianom fazowym (Rys. 8).W najniższych 
temperaturach agregaty lipidowe występują w fazie krystalicznej (L ). Ze wzrostem 
temperatury oddziaływania van der Waalsa stabilizujące L ulegają osłabieniu, pro- 
wadząc do otrzymania fazy żelowej (L,) [51]. Proces nazywany jest subprzejściem 
dla błon DPPC [57] i pojawia się przy 18 °C [58]. W fazie żelowej lipidy są słabiej 
upakowane, a przestrzeń pomiędzy główkami jest silniej uwodniona [51, 59]. Dalsze 
grzanie prowadzi do wzrostu ruchliwości łańcuchów węglowodorowych pozostają- 
cych nadal w konformacji trans. Po osiągnięciu temperatury zwanej temperaturą 
przedprzejścia, faza żelowa L, przechodzi w fazę pofałdowaną P g [51, 57, 60]. Dla 
DPPC przemiana pojawia się przy temperaturze 35°C [58] i nosi nazwę przedprzej- 
ścia. Faza P, charakteryzuje się periodycznie pofalowaną powierzchnią błony wyni- 
kającą ze wzrostu uwodnienia przestrzeni między głowami. W fazie żelowej pofał- 
dowanej dominującą konformacją łańcuchów węglowodorowych jest nadal trans, 
ale rośnie ich swobodna rotacja [60-62]. 


2.2. MECHANIZM ODDZIAŁYWANIA ANESTETYKU Z BŁONĄ LIPIDOWĄ 


Wszelkie procesy związane z obecnością leków w organizmie, z uwagi na 
ich dynamiczną naturę są wysoce złożone i trudne do zrozumienia. Badania nad 
oddziaływaniem lek — błona lipidowa pozwalają przewidywać właściwości farma- 
kokinetyczne substancji leczniczych. Ponadto mogą zostać wykorzystane do projek- 
towania i rozwijania skutecznych systemów dostarczania leków. 

Mechanizmy działania wielu leków opierają się na zmianie właściwości fizycz- 
nych błony lipidowej wynikających z bezpośredniego wiązania się z cząsteczkami 
lipidów [63-65]. Siła działania związków chemicznych, w tym anestetyków, zależy 
więc od zdolności do modyfikacji struktury i dynamiki dwuwarstwy lipidowej. 
Zmiany na poziomie cząsteczkowym, wywołane występowaniem leków w układzie, 
rzutują na właściwości makroskopowe błony takie jak: płynność [66], objętość [67] 
uporządkowanie [68], zwiększona przepuszczalność [69] czy przesunięcie tem- 
peratury głównego przejścia fazowego [70]. Różnicowa kalorymetria skaningowa 
(DSC) umożliwia monitorowanie przemian fazowych zachodzących w układach 
lipidowych [71]. W konsekwencji obserwowalny jest wpływ obecności anestetyków 
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na proces przejścia fazowego błony z fazy żelowej do ciekłokrystalicznej. Ponadto 
odnotowano zanik procesu przedprzejścia wraz ze wzrostem stężenia środków 
znieczulających w strukturze błony [72, 73]. Zbliżone działanie wykazują alkohole 
krótkołańcuchowe będące, podobnie jak wiele anestetyków, cząsteczkami amfifilo- 
wymi [73]. Modulacja temperatury głównego przejścia fazowego, któremu towarzy- 
szy redukcja grubości dwuwarstwy, jest wynikiem tworzenia się fazy przemiennej 
(ang. interdigitated), charakteryzującej się wzajemnym przenikaniem łańcuchów 
węglowodorowych lipidów z przeciwnych monowarstw [73-77]. 
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Rysunek 9. Diagram fazowy DPPC przy różnym stosunku woda/etanol [79] 
Figure 9. DPPC phase diagram at different water/ethanol ratio [79] 


Najdokładniej przebadanym związkiem wywołującym powstawanie fazy 
naprzemiennej jest etanol (Rys. 9). W błonach lipidowych w fazie żelowej wyka- 
zuje skłonności adsorbowania do grup obecnych w hydrofilowej głowie cząsteczek 
lipidu, zwłaszcza w pobliżu łańcuchów węglowodorowych. Uważa się, że cząsteczki 
amfifilowe tworzą wiązania wodorowe z grupami P=O oraz C=O reszty glicerolowej 
lipidu, powodując ich dehydratację. Na skutek oddziaływania grup OH alkoholu 
z resztami fosforanowymi formują się przestrzenie pomiędzy główkami lipidów, co 
w następstwie prowadzi do pojawiania się luk w rdzeniu hydrofobowym. Tworzone 
w regionie hydrofobowym wolne przestrzenie są niekorzystne energetycznie co pro- 
wadzi do utworzenie fazy przemiennej. Dzięki zabiegowi przenikania się łańcuchów 
wzmacniają się oddziaływania van der Waalsa, a entropia określana przez uporząd- 
kowane cząsteczki wody wzrasta co minimalizuje energię układu [77-81]. 

Jednym z głównych następstw tworzenia się fazy przemiennej jest efekt dwu- 
fazowy. Objawia się początkowym spadkiem wartości T „ a po osiągnięciu stężenia 
progowego alkoholu, stopniowym wzrostem T . Zjawisko związane jest z pojawie- 
niem się dwóch niezależnych rodzajów oddziaływań alkohol-lipid. Przy stosun- 
kowo niskich stężeniach, uprzywilejowane jest wbudowywanie się etanolu w błonę 
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lipidową w fazie ciekłej. Natomiast powyżej stężenia progowego, etanol oddziałuje 
z grupami chemicznymi obecnymi w główkach lipidów co powoduje stabilizację 
fazy żelowej, a w konsekwencji podwyższenie T [82, 83]. 

Obecność cząsteczek anestetyków w układzie wpływa na ciśnienie lateralne 
błony lipidowej. Wyniki oparte na obliczeniach teoretycznych dowodzą, że cząs- 
teczki o właściwościach znieczulających lokują się zarówno w części hydrofilowej 
jak i interfazie błony. Mechanizmy zapoczątkowane przez konkurencyjne tworzenie 
się wiązań wodorowych prowadzą do zmian w profilu ciśnienia lateralnego. Mody- 
fikacje właściwości błony przyczyniają się do zmian konformacyjnych białek błono- 
wych [23, 84-86]. 


3. GRAMICYDYNA JAKO MODEL BIAŁKA TRANSBLONOWEGO 


Zaproponowana przez Franks’ i Lieb’a [32, 87] teoria białkowo-zależna aneste- 
zji przewiduje, że wskutek bezpośredniego oddziaływania z białkami błonowymi, 
anestetyki na poziomie komórkowym wpływają na funkcję neuronów. Przewod- 
nictwo w aksonach oraz aktywność synaps są odpowiednio zależne od białkowych 
kanałów jonowych bramkowanych napięciem lub ligandem. Liczne badania wyka- 
zały, że zarówno kanały sodowe jak i potasowe, zaliczane do grupy zależnej od 
potencjału, nie wykazują szczególnej wrażliwości na działanie środków znieczulają- 
cych [33, 88]. Natomiast kanały aktywowane przez ligandy, a w szczególności przez 
glutaminian oraz GABA, zaliczane są do głównych celów wpływu anestetyków. 

Kwas y-aminomasłowy (GABA) jest najważniejszym hamującym neuro- 
przekaźnikiem w ośrodkowym układzie nerwowym. Pośredniczy w przepływie 
jonów chlorkowych przez błony synaptyczne, prowadząc do jej hiperpolaryzacji 
(Rys. 10).W efekcie zostaje zmniejszona pobudliwość neuronów oraz zahamowane 
przewodzenie impulsów [33, 89]. Receptory GABA mogą być modulowane przez 
liczną grupę substancji chemicznych [33, 35, 90, 91]. Wiele z tych związków zali- 
czanych jest do cząsteczek amfifilowych, których właściwości mogą również wpły- 
wać na strukturę i dynamikę otaczającej błony lipidowej. Modyfikacje strukturalne 
zachodzące w sąsiadujących z białkiem cząsteczkach lipidów często wpływają na 
jego aktywność. Pojawiające się oddziaływania hydrofobowe między regionem 
transbłonowym białka, a rdzeniem dwuwarstwy, wywołują miejscowe odkształcenia 
błony, które jednocześnie mogą wymusić zmiany konformacyjne łańcucha polipep- 
tydowego [92]. 
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Rysunek 10. (A) Schemat budowy receptora GABA, z podziałem na podjednostki [93]. (B) Struktura krysta- 
liczna ludzkiego receptora GABA [opracowane na podstawie danych bazy PDB] 

Figure 10. (A) GABAA receptor model [93]. (B) Crystal structure of GABA receptor [on the base on PDB 
data] 


Zrozumienie pełnego mechanizmu anestezji na poziomie molekularnym 
wymaga uwzględnienia w badaniach informacji o modyfikacjach strukturalnych 
kanałów jonowych. Mimo postępów w dziedzinie obrazowania biocząsteczek [94], 
analiza oddziaływań receptorów z anestetykami stanowi duże wyzwanie [95]. Roz- 
wiązaniem tego problemu jest zastosowanie dobrze zdefiniowanych układów mode- 
lowych jak albumina surowicy bydlęcej [96], albumina surowicy ludzkiej [97] czy 
lucyferaza muchy [98]. Rozpuszczalne w wodzie białka globularne dostarczyły wiele 
istotnych danych na temat oddziaływań białko-anestetyk. Jednak, aby odnieść się 
do działań środków znieczulających na strukturę białek błonowych wykorzystuje się 
modelowy kanał jonowy jakim jest Gramicydyna [99-101]. 


3.1. ZMIANY KONFORMACYJNE GRAMICYDYNY 
ZALEŻNE OD RODZAJU ROZPUSZCZALNIKA 


Gramicydyna to liniowy peptyd o właściwościach antybiotycznych. Wytwarzany 
przez bakterie glebowe Bacilus brevis, zbudowany jest z 15 reszt aminokwasowych 
L- i D-ułożonych naprzemiennie [102, 103]. W warunkach naturalnych Gramicy- 
dyna występuje w formie mieszaniny, której głównym składnikiem jest Gramicy- 
dyna A (85%). Ponadto obserwuje się obecność Gramicydyny B oraz C różniących 
się naturą aromatycznej reszty aminokwasowej w pozycji 11 [104] (Rys. 11). 
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Rysunek 11. Sekwencja aminokwasów w Gramicydynie A, B i C [105] 
Figure 11. Amino acid sequence of Gramicidin A, B and C [105] 


Unikalna sekwencja L- i D- aminokwasów sprawia, ze Gramicydyna wykazuje 
silną wrażliwość na działanie otaczającego środowiska. Obserwowalne jest przyjmo- 
wanie przez peptyd szerokiego spektrum zależnych od otoczenia konformacji. Stąd 
niejednokrotnie wykorzystywany jest jako model realistycznego określenia kon- 
formacyjnej preferencji białek w środowisku błony. Zidentyfikowano dwa główne 
motywy fałdowania Gramicydyny [102]: 

a) „głowa do głowy” 

b) podwójna helisa 

Konformacja podwójnej helisy dominuje w rozpuszczalnikach organicznych. 
Natomiast konformacja jaką przyjmuje Gramicydyna wbudowująca się do wnętrza 
błony lipidowej jest zależna od natury rozpuszczalnika, w jakim została wstępnie 
rozpuszczona (Rys. 12) [104]. 

Motyw podwójnej helisy, formuje się na skutek tworzenia się pomiędzy dwoma 
łańcuchami polipeptydowymi licznych, międzycząsteczkowych wiązań wodoro- 
wych. Następnie łańcuchy aminokwasów zwijają się w helisy. Tego typu struktura 
często określana jest symbolicznie jako struktura 27° lub zs", gdzie górny indeks 
wskazuje liczbę reszt przypadających na jeden skręt w helikalnym monomerze [102, 
104, 106]. Jak już zostało nadmienione, początkowa konformacja jaką przyjmuje 
Gramicydyna wbudowująca się w błonę zależy od rozpuszczalnika. Rozpuszczenie 
peptydu w roztworach metanol/chloroform, benzen/metanol lub etanolu prowadzi 
do uzyskania motywu podwójnej helisy uznawanej za formę nie-kanałową [104, 
107, 108]. 

Użycie trifluoroetanolu (TFE) jako rozpuszczalnika prowadzi do uzyskania 
formy kanałowej Gramicydyny po wbudowaniu się w strukturę błony lipidowej. 
W helikalnym dimerze, każdy z łańcuchów polipeptydowych tworzy serię wiązań 
wodorowych wewnątrzcząsteczkowych. Reszty aminokwasowe na N-końcu każdej 
helisy tworzą natomiast międzycząsteczkowe wiązania wodorowe dzięki czemu uzy- 
skuje się konformację dimeru „głowa do głowy” (B°"). Wnętrze kanału jest zbudo- 
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wane z polarnego szkieletu peptydowego. Na zewnątrz wyeksponowane są łańcuchy 
boczne aminokwasów, które są w kontakcie z łańcuchami acylowymi sąsiadujących 
lipidów oraz pomagają w regulowaniu przewodnictwa kanału. Rozmiar porów 
kanału jest wystarczająco duży by umożliwić przejście jednododatnim kationom 
[104, 108]. 


GRAMICYDYNA 


FORMA KANAŁOWA FORMA NIEKANAŁOWA 


Rysunek 12. Konformacje dimerów peptydowych Gramicydyny zależne od rodzaju rozpuszczalnika 
Figure 12. Solvent dependent Gramicidin dimer conformations 


3.2. WPŁYW ANESTETYKÓW NA STRUKTURĘ GRAMIDYCYNY 


Działanie środków znieczulających na białka błonowe może opierać się na bez- 
pośrednim wiązaniu się anestetyku do kanału jonowego lub na pośrednim wpływie 
na jego aktywność przez zaburzenie struktury i dynamiki błony lipidowej [109]. 

Próby identyfikacji miejsca wiązania się anestetyków do modelowego białka 
błonowego, jakim jest Gramicydyna, zostały podjęte przez Tanga i współpracow- 
ników, przy wykorzystaniu technik NMR [99], znakowania fotopowinnowactwa 
[101] czy modelowania molekularnego [110]. Silna zależność częstotliwości rezo- 
nansowej protonu amidowego w pierścieniu indolowym od stężenia halotanu 
pozwoliła ustalić, że strukturalnie selektywne oddziaływania pomiędzy Gramicy- 
dyną i środkiem znieczulającym występują w pobliżu reszt tryptofanowych pep- 
tydu przyjmującego konformację kanałową. Natomiast, reorganizacja dimeru do 
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formy podwójnej helisy prowadzi do niemal całkowitego zaniku tego zjawiska. Na 
tej podstawie wysnuto wniosek, że reszty tryptofanowe nie są jedynym czynnikiem 
determinującym oddziaływanie anestetyk — Gramicydyna, ale istotne jest również 
rozmieszczenie łańcuchów bocznych aminokwasów w dimerze peptydowym. Uni- 
kalna orientacja pierścieni indolowych względem otoczenia lipidowego, zwłaszcza 
Trp-9 i Trp-15 pozwala na utworzenie kieszeni amfifilowych preferowanych przez 
środki znieczulenia ogólnego. Ponadto przypuszcza się, że oddziaływanie Trp-9 
z wodą, utrudnione przez bliską obecność łańcucha bocznego Trp-15, jest możliwe 
po wprowadzeniu do układu cząsteczek anestetyków [99, 101]. 


UWAGI KOŃCOWE 


Najważniejszym wyzwaniem współczesnej anestezjologii jest optymalizacja 
dawkowania środków znieczulenia ogólnego. Powszechnie praktykowana jest tzw. 
anestezja złożona polegająca na połączeniu kilku specjalistycznych leków umożli- 
wiających ograniczenie niepożądanego wpływu anestetyków, przy jednoczesnym 
skutecznym znieczuleniu. Sprecyzowanie molekularnego mechanizmu działania 
środków znieczulenia ogólnego pozwoliłoby uzyskać adekwatną, dostosowaną do 
potrzeb pacjenta kontrolowaną utratą świadomości. Odwracalność zjawiska aneste- 
zji sugeruje, że analizowane zjawisko opiera się na zaburzeniu słabych oddziaływań 
międzycząsteczkowych, jak wiązanie wodorowe i siły van der Wallsa. Cząsteczki 
anestetyków mogą więc wiązać się bezpośrednio z hydrofobowymi regionami białka, 
wywołując jego zmiany konformacyjne bądź wywierać wpływ na lipidy otaczające 
receptory prowadzące do modyfikacji strukturalnych oraz organizacyjnych błony, 
które zaburzą aktywność białek. 

Wyniki uzyskane do tej pory wskazują na dwutorowy mechanizm znieczulenia 
ogólnego, w którym to anestetyki wnikają w obszar interfazy błony lipidowej, gdzie 
za pomocą obecnych w ich strukturze kwasowych wodorów, oddziałują zarówno z 
grupami karbonylowymi lipidów, jak i resztami tryptofanowymi białka. 
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ABSTRACT 


The question about the harmfulness of electromagnetic radiation is ever more 
often raised in scientific discussions it is also interesting for the public. Unfortuna- 
tely, there are no easy answers and no settled opinion. Followers of the catastrophic 
vision of the future argue, that the electromagnetic smog never grew as fast as today, 
and our bodies are probably not prepared for such conditions. Opponents, being 
accused of favoring telecommunication companies, claim that there are (or were) 
enough sources of electromagnetic radiation in the natural environment for organi- 
sms to learn how to self-protect, either by adaptation or development of specific sen- 
sory mechanisms. Questions are raised not only about health aspects, but also about 
security of data transmission and possible military applications, such as protection 
against radar localization. We want to present our recent investigations about the 
uses of natural materials for suppressing electromagnetic waves. The materials in 
question are ordinary clays, which role in Nature is barely now being discovered. 
Clay minerals, the basic ingredient of clay, are unusual materials, formed by Nature 
into nanostructures. Just recently we have begun to understand their role, behavior 
and possible applications. Those minerals, especially kaolinite and halloysite, are 
widespread, but in pure form occur only in few places in the world. Polish deposits 
belong to the richest ones. 


Keywords: clay, clay minerals, dielectric relaxation, structure, EMR materials 
Słowa kluczowe: glina, minerały ilaste, relaksacja dielektryczna, struktura, materiały 
EMR 
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HA 

HA-FA 
HA-NMFA 
HA-AA 
HA-DMSO 
HA-UREA 
HA-ACOK 
HA-FA-IMI 


HA-NMFA-EDTA 


HA-NMFA-AKA 


WYKAZ STOSOWANYCH SKROTOW 


promieniowanie elektromagnetyczne (ang. electro- 
magnetic radiation) 

względna zespolona przenikalność elektryczna 
rzeczywista składowa względnej zespolonej przeni- 
kalności elektrycznej 

ujemna składowa względnej zespolonej przenikalno- 
ści elektrycznej 

względna zespolona przenikalność magnetyczna 
rzeczywista składowa względnej zespolonej przeni- 
kalności magnetycznej 

ujemna składowa względnej zespolonej przenikalno- 
ści magnetycznej 

współczynnik (stała) tłumienia 

haloizyt 

interkalat haloizytu z formamidem 

interkalat haloizytu z N-metyloformamidem 
interkalat haloizytu z acetamidem 

interkalat haloizytu z dimetylosulfotlenkiem 
interkalat haloizytu z mocznikiem 

interkalat haloizytu z octanem potasu 

interkalat wtórny haloizytu z imidazolem (otrzymany 
z HA-FA) 

interkalat wtórny daloizytu z etylenodiaminą (otrzy- 
many z HA-NMFA) 

interkalat wtórny haloizytu z akryloamidem (otrzy- 
many z HA-NMFA) 
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WPROWADZENIE 


Pytanie o szkodliwość promieniowania elektromagnetycznego coraz częściej 
przewija się w dyskusjach prowadzonych w środowiskach naukowych, jest też 
interesujące dla opinii publicznej. Odpowiedź nie jest niestety łatwa i mimo podej- 
mowanych prób brak rozstrzygających dowodów. Zwolennicy katastroficznej wizji 
przyszłości argumentują, że smog elektromagnetyczny nigdy nie rósł tak szybko jak 
obecnie i prawdopodobnie nasze organizmy nie są przygotowane na takie warunki 
życia. Oponenci, będąc posądzani o sprzyjanie koncernom telekomunikacyjnym, 
twierdzą, że w środowisku naturalnym jest (lub było) dostatecznie wiele źródeł pro- 
mieniowania elektromagnetycznego, aby organizmy mogły się na nie uodpornić, 
lub, jeśli jest ono niebezpieczne, wykształcić sensory ostrzegające przed nim. Obok 
aspektów zdrowotnych podnoszone są również pytania o bezpieczeństwo transmisji 
danych lub propozycje zastosowań militarnych, jak choćby zabezpieczenie przed 
lokalizacją radarową. Przegląd, który pragniemy tutaj zaprezentować dotyczy próby 
zastosowania materiałów naturalnych lub ich chemicznych modyfikacji do tłumie- 
nia fal elektromagnetycznych. Materiały, o których mowa, to zwykła z pozoru glina, 
której znaczenie w przyrodzie i właściwości dopiero teraz odkrywamy. Minerały 
ilaste, będące podstawowym składnikiem gliny są niezwykłymi materiałami, które 
przyroda uformowała w naturalne nanostruktury. Dopiero nauka XX i XXI wieku 
zaczyna rozumieć ich rolę, sposób działania i możliwości. Minerały te, a zwłasz- 
cza kaolinit i haloizyt, bo nimi będziemy się tu zajmować, mimo że szeroko roz- 
powszechnione, to jednak w postaci czystej występują tylko w kilku miejscach na 
świecie. Polskie złoża należą do najbogatszych. 


1. HALOIZYT, OGÓLNA CHARAKTERYSTYKA I GENEZA 


Haloizyt jest dioktaedrycznym minerałem ilastym, glinokrzemianem warstwo- 
wym typu 1:1 należącym do grupy kaolinitu — serpentynu, do podgrupy kaolinitu 
[1]. Określenie „dioktaedryczny” oznacza, że w warstwie oktaedrycznej ośmiościany 
zajmują dwa miejsca na trzy możliwe przypadające na jednostkę strukturalną, trze- 
cie miejsce jest wolne. Typ 1:1 określa budowę pakietu haloizytu: na jedną warstwę 
oktaedryczną przypada jedna warstwa tetraedryczna, czyli pakiet jest zbudowany 
z tych dwu warstw ze sobą połączonych. Minerał ten został po raz pierwszy opi- 
sany przez Berthiera w 1826 roku [2], natomiast Hofmann ze współpracownikami 
[3] wykazał jego krystaliczny charakter. Haloizyt jest głównym składnikiem gleb 
pochodzenia wulkanicznego [4]. 

Skład chemiczny haloizytu zależy od miejsca jego występowania. Stąd też 
próbki tego minerału pochodzące z różnych miejsc będą różniły się zawartością pro- 
centową poszczególnych składników. Wzór sumaryczny minerału będzie zależał od 
tego, o jakiej odmianie haloizytu jest mowa. Minerał ten występuje bowiem w przy- 
rodzie jako mieszanina hydrohaloizytu i bezwodnego haloizytu. Wzór sumaryczny 
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hydrohaloizytu można przedstawić jako Al,Si,O.(OH),:2H,O, a haloizytu jako 
ALSi,O.(OH), [5]. Pod względem budowy haloizyt i kaolinit są podobne do siebie. 
Różnią się tym, iż w haloizycie powierzchnie są zakrzywione lub zwinięte wskutek 
niedopasowania warstw oktaedrycznych i tetraedrycznych. W skali molekularnej 
warstwa silikatowa i grupy OH warstwy gibbsytowej obu minerałów są identyczne 
[6]. Podstawową jednostką budulcową haloizytu jest pakiet. Każdy pakiet składa się 
z dwóch warstw: tetraedrycznej warstwy silikatowej [Si,O,]” i oktaedrycznej gib- 
bsytowej [Al,(OH),]**. Nawarstwiajace się pakiety dają pełną strukturę minerału. 
Warstwa silikatowa składa się z tetraedrów (w centrum atom krzemu, w narożach 
— atomy tlenu), które łączą się narożami z trzema sąsiednimi tetraedrami w ten spo- 
sób, że ich środki ciężkości tworzą zdeformowane sześcioboki. Czwartym atomem 
tlenu warstwa silikatowa łączy się z warstwą oktaedryczną — tak więc obie warstwy 
są ze sobą związane poprzez wspólny atom tlenu [1]. Warstwa gibbsytowa jest zbu- 
dowana z oktaedrów, w których środku znajduje się atom glinu, a w narożach cztery 
grupy hydroksylowe i dwa atomy tlenu. Pakiety są połączone wiązaniami wodoro- 
wymi między tetraedrycznymi tlenami warstwy silikatowej i zewnętrznymi grupami 
hydroksylowymi warstwy gibbsytowej (z przyległego pakietu). Jedna grupa wodo- 
rotlenowa (OH,,) leży między warstwami tego samego pakietu, natomiast pozostałe 
trzy (OH,,,) znajdują się pomiędzy dwoma sąsiednimi pakietami i tworzą wiązania 
wodorowe między warstwą silikatową i warstwą gibbsytową sąsiedniego pakietu 
[7, 8]. Jak już wspomniano, naturalny haloizyt może występować jako mieszanina 
hydrohaloizytu (odstęp międzypłaszczyznowy d = 10 + 0,2 A, zwany haloizytem 
10A) zawierającego słabo związane międzypakietowe cząsteczki wody i haloizytu 
(odstęp międzypłaszczyznowy d = 7 + 0,2 A, zwany haloizytem 7 A) bez między- 
pakietowej wody. Obecność wody międzypakietowej w haloizycie 10 A jest istotną 
cechą wyróżniającą ten minerał spośród innych przedstawicieli podgrupy kaolinitu 
(kaolinit, dickit, nakryt). Jej zawartość w hydrohaloizycie wynosi ok. 12,3% [9]. 
Hydrohaloizyt jest bardzo niestabilny, łatwo i nieodwracalnie traci wodę międzypa- 
kietową. Przy łagodnym ogrzaniu lub nawet w temperaturze pokojowej w powietrzu 
o małej wilgotności przechodzi w haloizyt 7 A. Nie można haloizytu 7 A przekształ- 
cić w jego uwodnioną formę poprzez np. mieszanie go z wodą. Istnieje jednakże 
metoda odzyskiwania hydrohaloizytu. Polega ona na interkalowaniu haloizytu 7 A 
octanem potasu, a następnie na wymianie octanu potasu na wodę podczas długo- 
trwałego mieszania interkalatu z wodą. W ten sposób otrzymać można haloizyt 10 A 
jako tzw. interkalat wtórny [10]. Cząsteczki wody międzypakietowej w hydrohalo- 
izycie tworzą wiązania wodorowe między sobą, a również z tlenami warstwy sili- 
katowej i z grupami hydroksylowymi warstwy gibbsytowej [11]. Costanzo i współ- 
pracownicy zidentyfikowali dwa typy wody w haloizycie: wodę zaadsorbowaną na 
powierzchni minerału oraz zasocjowaną wodę międzypakietową [12-14]. 

Haloizyt przyjmuje zwykle postać wydłużonych rurek, jednak tworzy również 
rurki krótkie, formy sferyczne i płaskie [2]. Morfologia tego minerału zależy od 
warunków krystalizacji i występowania geologicznego [15]. Średnica form kulistych 
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wynosi od 0,05 do ok. 0,50 um. Długość rurek zawiera się w przedziale od 0,02 do 
ponad 30 um, a ich grubość wynosi od ok. 0,05 do 0,2 um [16]. Morfologię rurkową 
posiada np. polski haloizyt z kopalni „Dunino” - ok. 30% struktury polskiego mine- 
rału to proste nanorurki, które znajdują się między pakietami płytek [17]. Pomię- 
dzy morfologicznymi odmianami haloizytu występują różnice w strukturze. Rurki 
haloizytowe buduje pięć „zrolowanych” pakietów [18], natomiast formy sferyczne 
przypominają strukturę cebuli [19]. 


2. WŁAŚCIWOŚCI FIZYKOCHEMICZNE 


Powierzchnia właściwa haloizytu zależy od pochodzenia próbki. Haloizyt 
z kopalni „Dunino” na Dolnym Śląsku ma powierzchnię właściwą równą 60.9 mg” 
[17], haloizyty z Nowej Zelandii i Australii - od 50 do 137 me" [19, 20]. Pojemność 
wymiany kationów haloizytu zawiera się w przedziale od 2 do 60 cmol(+)/kg i jest 
większa od pozostałych minerałów z podgrupy kaolinitu (1-15 cmol(+)/kg) [2]. 
Polski haloizyt wykazuje pojemność wymiany kationów na poziomie 55 cmol(+)/kg 
[17]. Pojemność wymiany kationów dla hydrohaloizytu jest większa niż dla haloizytu 
7 A [2]. Stwierdzono również, że wskutek obecności w rurkach kationów występu- 
jących w postaci rozpuszczalnych soli, rurkowe formy haloizytu osiągają większe 
wartości pojemności kationowej od pozostałych form [21]. Na pojemność jonową 
haloizytu mają wpływ niewysycone wiązania na krawędziach blaszek minerału. Im 
mniejszy jest rozmiar blaszek, tym większy jest wzrost pojemności wymiany katio- 
nów [22, 23]. Ładunek ten jest równoważony przez uwodnione kationy, sąsiadujące 
z tetraedryczną stroną przestrzeni międzypakietowej. Podstawienia w warstwie gib- 
bsytowej również powodują wzrost pojemności wymiany kationów (zarówno sub- 
stytucje jonów Al” jonami Fe” , jak i jonów Al” jonami Fe*’) [15]. 

Interesującą własnością minerałów ilastych jest ich zdolność do oddziaływa- 
nia z substancjami nieorganicznymi i organicznymi. Może to zachodzić na dro- 
dze adsorpcji na powierzchni ziaren minerału, wymiany kationów lub interkalacji 
[24]. Interkalacja haloizytu polega na wnikaniu cząsteczek gościa między pakiety 
minerału, co powoduje zerwanie wiązań wodorowych pomiędzy tlenami warstwy 
tetraedrycznej a wodorami warstwy oktaedrycznej i powstanie nowych wiązań mię- 
dzy gospodarzem (minerałem) a gościem (interkalantem organicznym lub nieorga- 
nicznym). Proces ten powoduje zwiększenie odległości między pakietami wzdłuż 
krystalograficznej osi c, co jest obserwowane na dyfraktogramach proszkowych. 
Interkalaty nie są stechiometrycznymi związkami i można je zaliczyć do bertolidów. 
Powstałe w wyniku omawianego procesu interkalaty mają inne właściwości fizyko- 
chemiczne niż wyjściowe materiały. 

Interkalacja solami jonowymi nosi nazwę intersalacji. W przypadku hydroha- 
loizytu proces polega na zastąpieniu wody międzypakietowej przez sól. Jony soli 
formują pojedynczą warstwę w przestrzeni między sąsiednimi pakietami. Niektóre 
sole, jak np. octan potasu, mogą też tworzyć warstwę podwójną. Interkalaty haloizytu 
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10 A zawierające sole jonowe otrzymuje się poprzez zwykłe mieszanie minerału ze 
stężonymi wodnymi roztworami soli, w których rolę kationu pełnią najczęściej 
jony jednowartościowe. Interkalat haloizytu 10 A i haloizytu 7 A z octanem potasu 
można również uzyskać w wyniku procesu mechanochemicznego, polegającego na 
ucieraniu minerału ze stałym interkalantem. Intersalacja zachodzi łatwo, jeśli sole 
spełniają następujące warunki [10]: a) składają się z dużych kationów o małej energii 
hydratacji (K*, NH}, Rb’, Cs"); b) zbudowane są z dużych anionów (np. CI, POJ, 
CH,COO, C,H,COO ). Kationy soli mogą być uwięzione w di-trygonalnych lukach 
warstwy silikatowej. Aniony ustawiają się dookoła schowanego kationu, neutralizu- 
jąc w ten sposób nadmiar ładunku dodatniego [2]. 

Związki organiczne posiadające zdolność do interkalacji haloizytu są polarne 
i zawierają grupy donorowe i akceptorowe (jedną z nich lub obie naraz). Dzięki 
temu mogą tworzyć wiązania wodorowe z warstwą gibbsytową i/lub silikatową. 
Interkalujące do przestrzeni międzypakietowej cząsteczki gościa mogą się z nim 
wiązać na różne sposoby, w zależności od rodzaju posiadanych grup funkcyjnych. 
Interkalanty podzielono na trzy grupy ze względu na sposób ich oddziaływania 
z minerałem [25]. Do pierwszej grupy należą molekuły, które tworzą silne wiązania 
wodorowe z warstwą krzemo-tlenową, np. hydrazyna, mocznik, formamid, aceta- 
mid. W drugiej grupie są cząsteczki z silnymi oddziaływaniami dipolowymi, które 
mogą oddziaływać z warstwą krzemotlenową, np. dimetylosulfotlenek. Grupę trze- 
cią stanowią sole metali alkalicznych z krótkołańcuchowymi kwasami karboksylo- 
wymi, przede wszystkim z kwasem octowym i propionowym. 


* a, béi + DMSO 


Rysunek 1. Sposoby interkalacji dimetylosulfotlenku między warstwy halloizytu 
Figure 1. Ways of intercalation of dimethylsofoxide into haloysite layers 
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Interkalacja molekuł do haloizytu 10 A jest bardziej zależna od własności 
interkalantów niż cech minerału [2], z kolei interkalacja haloizytu 7 A zależy przede 
wszystkim od własności minerału: jego pochodzenia i morfologii, zawartości żelaza, 
obecności i typu zanieczyszczeń [16, 26]. Duże rozmiary blaszek, sfer czy rurek 
haloizytu i mała zawartość żelaza sprzyjają interkalacji [2]. Wynika z tego, że różne 
próbki haloizytu będą wykazywały różną skłonność do interkalacji. Na Rysunku 1 
przedstawiono schematycznie proces interkalacji cząsteczek DMSO do hydrohalo- 
izytu i haloizytu. 

Wyróżniamy zatem interkalaty pierwotne i wtórne. Interkalaty pierwotne 
- powstają poprzez bezpośrednie interkalowanie cząsteczek pomiędzy pakiety 
haloizytu; haloizyt tworzy interkalaty pierwotne np. z mocznikiem, dimetylosul- 
fotlenkiem, octanem potasu, formamidem, N-metyloformamidem, acetamidem. 
Interkalaty wtórne — powstają w wyniku wymiany wcześniej zainterkalowanych 
molekuł przez inne molekuły. Iterkalację wtórną stosujemy jeśli nie udaje się otrzy- 
mać interkalatu w procesie bezpośrednim. Z interkalatu haloizytu z N-metylofor- 
mamidem można uzyskać interkalaty z akryloamidem i etylenodiaminą; z interka- 
latu z formamidem otrzymuje się interkalat z imidazolem. Jako interkalat wtórny 
uznać można również hydrohaloizyt uzyskany z interkalatu z octanem potasu po 
wymianie cząsteczek octanu potasu na wodę. 


3. WYBRANE ZAGADNIENIA TEORII ABSORBERÓW EMR 


Absorber EMR to materiał pochłaniający fale elektromagnetyczne. Przy opisie 
takich substancji podstawowe znaczenie mają dwie wielkości: przenikalność elek- 
tryczna e* = e - je” i przenikalność magnetyczna H" = W — jw. Są to liczby zespo- 
lone, gdzie j jest jednostką urojoną, € i w są składowymi rzeczywistymi związanymi 
z gromadzeniem energii, natomiast e” iu” to składowe urojone opisujące dyssypację 
energii elektrycznej i magnetycznej, wynikającej z przewodnictwa, rezonansu i zja- 
wisk relaksacyjnych. Głównymi wymaganiami, które winien spełniać dobry absor- 
ber EMRsą [27-32]: a) minimalizacja odbicia promieniowania od jego powierzchni; 
b) wysokie wartości e” i u”; c) absorpcja promieniowania EMR w szerokim zakresie 
częstotliwości; d) brak konieczności użycia zewnętrznego pola magnetycznego lub 
elektrycznego. Do opisu właściwości absorbujących substancji używa się różnych 
wielkości fizycznych. Dla materiałów dielektrycznych najlepiej stosować stałą tłu- 
mienia a, będącą rzeczywistą składową współczynnika propagacji y [27]: 


2 
a = Re(y) = Re (IŻ Jue) (1) 


gdzie f jest częstotliwością fali, c prędkością światła w próżni. Uwzględniając część 
rzeczywistą i urojoną przenikalności elektrycznej i magnetycznej dla materiałów 
niemagnetycznych (œ = 1, u” = 0) wzór przybiera prostą postać: 
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2 , 
a = AG +ye?+ e) (2) 

Im większa stała tłumienia o materiału, tym lepsze jego właściwości absorbu- 
jące. Jak wynika z powyższego wzoru, a jest największe przy dużych wartościach e” 
i małych e. Należałoby, zatem poszukiwać materiałów o takiej właśnie charaktery- 
styce. W rzeczywistych dielektrykach wysokie wartości e” są związane z wysokimi 
wartościami e (co obniża a). 

Absorbery EMR mają wiele potencjalnych zastosowań. Oprócz ochrony przed 
niekorzystnym działaniem niejonizującego promieniowania elektromagnetycz- 
nego na środowisko i organizmy żywe (elektromagnetyczna kompatybilność dla 
urządzeń elektronicznych, redukcja tła elektromagnetycznego), można je wyko- 
rzystać w sektorze wojskowym i cywilnym. Dzięki pochłanianiu promieniowania 
elektromagnetycznego absorbery EMR są niewidoczne dla radarów, mogą więc być 
z powodzeniem stosowane jako materiały maskujące. Materiały EMR znalazły też 
zastosowanie w konstruowaniu komór ciszy elektromagnetycznej i ochrony kabli 
przesyłających dane [27]. Badania nad absorberami EMR są ważną gałęzią współ- 
czesnej nauki o materiałach. Konstruowane absorbery są oparte na kompozytach 
polimerowych wypełnionych alotropowymi odmianami węgla, kompozytach fer- 
rytu strontu (SrFe,,O,,), sadzy i kauczuku nitrylowego [33], czy też zmodyfikowa- 
nym grafenie [35]. Właściwości dielektryczne i magnetyczne kompozytów poli- 
merowych wypełnionych alotropowymi odmianami węgla mogą być dosyć łatwo 
modyfikowane poprzez zmianę węglowego wypełniacza lub matrycy polimerowej. 
Właściwości matryc polimerowych zależą nie tylko od składu chemicznego, ale też 
od morfologii wypełniacza węglowego i jego składu. Dodatkowym atutem materia- 
łów polimerowych jest ich mały ciężar i łatwość formowania. Polimery z wypełnie- 
niem opartym na alotropowych odmianach węgla wykazują zdolność absorpcji pro- 
mieniowania elektromagnetycznego w stosunkowo dużym zakresie częstotliwości. 
Materiały te są jeszcze w fazie badań laboratoryjnych [27]. 


4. BADANE MATERIAŁY 


Haloizyt, którego badania opisano w niniejszej pracy pochodził z kopalni 
„Dunino” na Dolnym Śląsku. Właścicielem kopalni jest firma PTH INTERMARK. 
Kopalnia odkrywkowa „Dunino” to jedna z trzech kopalni haloizytu na świecie 
(pozostałe dwie znajdują się w Stanach Zjednoczonych i Nowej Zelandii). Polskie 
złoża haloizytu są szacowane na 10-12 mln ton. Minerał używany do badań został 
nieodpłatnie przekazany przez firmę PTH INTERMARK. Firma ta prowadzi prze- 
róbkę kopalin. Do badań relacjonowanych w tym przeglądzie używano zarówno 
haloizytu naturalnego, jak i minerału oczyszczonego. Badanie skały surowej niesie 
niebezpieczeństwo, że obserwowane właściwości będą spowodowane obecnością 
niekontrolowanych zanieczyszczeń. Haloizyt, zgodnie z opisem zamieszczonym 
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wcześniej, jest materiałem nietrwałym. Szczególnie dotyczy to jego uwodnionej 
odmiany — hydrohaloizytu. Z tego względu dołożono starań, aby badane próbki nie 
były w istotny sposób modyfikowane w trakcie przygotowania. W czasie oczyszcza- 
nia naturalna skała była rozcierana, wielokrotnie przemywana wodą i kilkakrotnie 
bezwodnym etanolem. Po odwirowaniu materiał suszono w temperaturze pokojo- 
wej. 

Podczas otrzymywania interkalatów używano haloizytu oczyszczanego w wyżej 
opisany sposób. Interkalaty haloizytu otrzymywano głownie metodą interkalacji 
z cieczy, w niektórych zaś przypadkach stosowano również metodę mechanoche- 
miczną. 

W celu potwierdzenia skuteczności interkalacji stosuje się rutynowo rentge- 
nografię proszkową (XRD). Minerał przyjmujący do przestrzeni międzypakietowej 
cząsteczki gościa zwiększa stałą sieciową d,,, i jest to dowodem na powodzenie pro- 
cesu. Rysunek 2 przedstawia zestawienie wybranych widm proszkowych. Badania te 
wykonano na Wydziale Chemii Uniwersytetu Wrocławskiego przy wykorzystaniu 
dyfraktometru proszkowego D8 ADVANSE i promieniowania o długości 1,5418 A. 


! units 
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intensit 
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Rysunek 2. Widmo proszkowe XDR: 1 - naturalny haloizyt, 2 — haloizyt suszony, 3 - interkalat z acetamidem, 
4 — interkalat z octanem potasu, 5 — interkalat z dimetylosulfotlenkiem, 6 — interkalat z formami- 
dem, 7 — interkalat z mocznikiem, 8 — interkalat z N-metyloformamidem 

Figure 2. XRD spectrum of: 1- natural haloysite, 2 — dryed haloysite, 3 - intercalate with acetamide, 4 - 
intercalate with potassium acetate, 5 — intercalate with dimethylsulfoxide, 6 — intercalate with 
formamide, 7 — intercalate with urea, 8 - intercalate with N-methylformamide 
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Dla ułatwienia porównania otrzymanych rezultatów otrzymane zależności ska- 
lowano do sygnału od d. odjęto również tło. 

Otrzymane interkalaty były mieszaniną wyjściowego materiału (haloizytu) 
i interkalowanych pochodnych. Zawartość interkalatów była oszacowana na podsta- 
wie porównania intensywności sygnałów pochodzących od niezmienionego mate- 
riału (d,,, ok. 7 A) i interkalatu. Oszacowane w ten sposób wydajności zilustrowane 
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Rysunek 3. Wydajność interkalacji: HA-FA haloizyt z formamidem, HA-NMFA haloizyt z N-metyloforma- 
midem, HA-AA haloizyt z acetamidem, HA-DMSO haloizyt z dimetylosulfotlenkiem, HA-UREA 
haloizyt z mocznikiem, HA-FA-IMI interkalat wtórny haloizyt-formamid-imidazol, HA-NMFA- 
ETDA interkalat wtórny haloizyt-N-metyloformamid-etylenodiamina, HA-NMFA-AKA inter- 
kalat wtórny haloizyt-N-metyloformamid-akryloamid, HA-ACOK haloizyt z octanem potasu 

Figure 3. Efficiency of intercalation: HA-FA halloysite with formamide, HA-NMFA halloysite with 
N-methylformamide, HA-AA halloysite with acetamide, HA-DMSO halloysite with dimethylosul- 
foxide, HA-UREA halloysite with urea, HA-FA-IMI secondary intercalate halloysite-formamide- 
imidazole, HA-NMFA-ETDA secondary intercalate halloysite-N-methylformamide-ethyleno- 
diamine, HA-NMFA-AKA secondary intercalate halloysite-N-methylformamide-acryloamide, 
HA-ACOK halloysite with potassium acetate 


Widma proszkowe XRD pozwalały nie tylko na kontrolę wydajności procesu 
interkalacji, ale również na wyznaczenie przyrostu odległości międzypakieto- 
wej w wyniku wprowadzenia cząsteczek gościa. Zestawienie rozsunięć liczonych 
jako różnica odległości d,,, dla interkalatu i dla haloizytu przedstawione jest na 
Rysunku 4. Najmniejsze rozsunięcie pakietów zaobserwowano dla interkalatu z for- 
mamidem, największe dla interkalatu z imidazolem. 
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Rysunek 4. Rozsuniecie pakietów haloizytu spowodowane interkalacją. Oznaczenia takie same jak na 
Rysunku 3 
Figure 4. Displacement of halloysite layers due to intercalation. Captions are the same as in Figure 3 


Porównanie wymiarów gościa i rozsunięcia pakietów jest wskazówką dotyczącą 
sposobu ulokowania interkalanta w warstwie glinokrzemianu, w szczególności 
sposobu jego wiązania z warstwami. Pozwala to również na wnioskowanie co do 
swobody ruchów cząsteczki gościa pod wpływem zewnętrznego pola elektrycznego. 
Rysunek 5 przedstawia porównanie rozsunięcia pakietów z osią długą cząsteczek 
(obliczoną za pomocą programu Hyperchem). 
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Rysunek 5. Porównanie rozsunięcia warstw haloizytu z osią długą interkalowanych cząsteczek 
Figure 5. Comparison of displacement of halloysite layers and long axis of the intercalated molecules 
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Jak widzimy, rozsuniecie warstw dobrze koreluje z wymiarami cząsteczek. Od 
korelacji odbiega jedynie octan potasu, dla którego rozsunięcie warstw jest większe 
niż wynikałoby to z jego rozmiarów. Przypuszczalnie jest to związane z jonowym 
charakterem octanu potasu i dalekim zasięgiem oddziaływań kulombowskich. Przy- 
puszczamy, że warstwy są utrzymywane dzięki przyciąganiu między jonem potasu 
ulokowanym we wnęce warstwy gibbsytowej i jonem octanowym związanym z war- 
stwą glino-tleno-wodorotlenową. 


5. HALOIZYT I JEGO INTERKALATY JAKO ABSORBERY 
PROMIENIOWANIA ELEKTROMAGNETYCZNEGO 


Aby ocenić zdolność materiału do tłumienia fal elektromagnetycznych 
konieczne jest poznanie przenikalności elektrycznej i magnetycznej, i ich zależności 
od częstotliwości. Haloizyt i jego interkalaty mogą wykazywać pewne właściwości 
magnetyczne, założyliśmy jednak, że przenikalność magnetyczna badanych mate- 
riałów jest praktycznie równa przenikalności próżni. Badania dielektryczne prowa- 
dzono w zakresie częstotliwości od 0,1 kHz do 1 MHz i dla temperatur w zakresie od 
ok. -150 do 120°C mierząc rzeczywistą (£) i urojoną (e) składową przenikalności 
elektrycznej [36, 37]. Pomiary prowadzono zarówno dla próbek proszkowych, sto- 
sując technikę zastosowaną wcześniej dla kaolinitu [38], jak i dla próbek uformo- 
wanych w postaci pastylek. Dla próbek proszkowych doświadczalne wartości prze- 
nikalności elektrycznej korygowano ze względu na upakowanie. Stwierdzono, że 
dla haloizytu wartości przenikalności maleją w szeregu: haloizyt naturalny, haloizyt 
przemywany wodą i suszony na powietrzu, haloizyt przemywany etanolem i suszony 
na powietrzu. Spadek ten jest prawdopodobnie związany ze zmniejszeniem ilości 
wody międzyziarnowej, wody zaadsorbowanej na powierzchni ziaren, ale przede 
wszystkim wody międzypakietowej, która — jak się wydaje — ma decydujące zna- 
czenie dla obrazu dielektrycznego haloizytu. Proces relaksacyjny zaobserwowano 
w próbkach haloizytu, jak również interkalatu haloizytu z formamidem i N-mety- 
loformamidem. W pozostałych interkalatach w postaci proszków nie stwierdzono 
relaksacji w badanym przedziale częstotliwości. 

Dla pastylek otrzymano nieco inny obraz zjawisk dielektrycznych. Dla prosz- 
kowych próbek haloizytu naturalnego przenikalność elektryczna mierzona dla np. 
200 Hz w temperaturze pokojowej wynosiła ok. 1300 jednostek, tymczasem w mate- 
riale przygotowanym w postaci pastylki w podobnych warunkach przenikalność 
osiąga zaledwie 100 jednostek. Przykładowy rezultat badań przeprowadzonych dla 
haloizytu suszonego przed pomiarem przedstawiają Rysunki 6 i 7. 
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Rysunek 6. ` Rzeczywista składowa przenikalności elektrycznej suchego haloizytu 
Figure 6. Real component of the electric permittivity of dry halloysite 
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Rysunek 7. _ Ujemna urojona składowa przenikalności elektrycznej suchego haloizytu 
Figure 7 Negative imaginary component of the electric permittivity of dry halloysite 


Relaksację niskoczęstotliwościową tłumaczymy polaryzacją przestrzenną zwią- 
zaną z migracją ładunków powierzchniowych obecnych na ziarnach. Przypusz- 
czamy, że dla materiału sypkiego ruch ładunków ograniczony jest do rozmiarów 
ziarna, natomiast proces prasowania próbek prowadzi do pojawienia się kontaktów 
między ziarnami, pozwalając na migrację ładunków na znacznie większe odle- 
głości. Duża droga migracji musi powodować przesunięcie dyspersji do niższych 
częstotliwości leżących prawdopodobnie poza zakresem stosowanym w badaniach. 
Przypuszczamy, że jest to główna przyczyna dużych różnic pomiędzy rezultatami 
otrzymanymi dla materiałów proszkowych i pastylek. Silny wzrost przenikalności, 
podobny do tego obserwowanego dla próbek proszkowych, powinien być obecny dla 
niższych częstotliwości. Postulowane przesunięcie obszaru występowania dużych 
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polaryzacji do niskich częstotliwości otwiera możliwość badania procesów mole- 
kularnych maskowanych poprzednio polaryzacją przestrzenną i powierzchniową. 
Dla wszystkich materiałów badanych w postaci pastylek zaobserwowano relaksa- 
cję dielektryczną w przedziale częstotliwości 0,1-1000 kHz. Rysunek 8 przedstawia 
porównanie częstotliwości relaksacyjnych, a Rysunek 9 inkrementów dielektrycz- 
nych obserwowanych procesów. 
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Rysunek 8. Częstotliwość relaksacyjna haloizytu (HA przemyty), haloizytu suszonego (HA suszony) i bada- 
nych interkalatów. Oznaczenia takie jak na Rysunku 3 

Figure 8. Relaxation frequency of halloysite (HA przemyty), dry halloysite (HA suszony) and intercalates. 
Captions are the same as in Figure 3 
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Rysunek 9.  Inkrementy dielektryczne dla haloizytu suszonego (HA suszony) i badanych interkalatów. 
Oznaczenia takie jak na Rysunku 3 

Figure 9. Dielectric increments of dry halloysite (HA suszony) and intercalates. Captions are the same as in 
Figure 3 


f, [kHz] 


ch = 
© 
© 


CO 
© 


d 
© 


a 
© 


dielectric increment 


HALOIZYT - NATURALNY ABSORBER PROMIENIOWANIA ELEKTROMAGNETYCZNEGO 303 


Widzimy, że wprowadzenie cząsteczek gościa znacznie zwiększa średni czas 
relaksacji (zmniejsza częstotliwość relaksacyjną). Od schematu odbiega octan 
potasu, w tym jednak przypadku można to tłumaczyć obecnością międzypakietowej 
wody, która zwykle towarzyszy octanowi potasu w procesie interkalacji. Korelacja 
między wzrostem średniego czasu relaksacji i obecnością cząsteczek gościa jest ewi- 
dentna, wydaje się jednak, że istotną rolę w relaksacji odgrywa też woda, i to nie 
tylko w przypadku octanu potasu. Wnioskować można o tym na podstawie dziw- 
nych zmian czasu relaksacji dla wyższych temperatur. 

Dla większości badanych materiałów, zwłaszcza badanych w postaci pastylek, 
ale też jako proszki, zaobserwowano, że w wyższych temperaturach czas relaksacji 
rośnie wraz ze wzrostem temperatury. Ilustruje to Rysunek 10. 
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Rysunek 10. Wykres Eyringa dla haloizytu: logarytm iloczynu temperatury i czasu relaksacji w funkcji odwrot- 
ności temperatury 

Figure 10. Eyring plot for halloysite: logarithm of a product of temperature and relaxation time in function 
of inverse of temperature 


Jest to sprzeczne ze wzorem Eyringa, sugeruje bowiem ujemną energię akty- 
wacji. Za obserwowany efekt „cofania się” maksimum absorpcji przy wzroście tem- 
peratury może odpowiadać utrata wody znajdującej się w przestrzeniach między- 
warstwowych. Powoduje to zmiany strukturalne w obszarze warstwy, a to powinno 
skutkować zwiększeniem energii aktywacji dla ruchów cząsteczek aktywnych 
w relaksacji dielektrycznej. Pogląd ten potwierdzają widma termogawimetryczne 
pokazujące spadek masy próbek w przedziale temperatur, w których dochodzi 
do anomalnych zmian czasów relaksacji. Wspomnianych zmian czasu relaksacji 
nie zaobserwowano dla interkalatów haloizytu z akryloamidem, etylenodiaminą 
i N-metyloformamidem. Przypuszcza się jednak, że są one obecne w wyższych tem- 
peraturach niż stosowane w eksperymencie. Dla wyjaśnienia anomalnych zmian 
czasu relaksacji zaproponowano model, w którym energia aktywacji jest funkcją 
temperatury. Jeśli wraz ze wzrostem temperatury przebiega proces starzenia próbki 
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lub zmiany jej składu i związane z tym utrudnienie ruchów cząsteczek odpowie- 
dzialnych za relaksację, wówczas energia aktywacji rośnie, a to powoduje obser- 
wowany wzrost czasu relaksacji ze wzrostem temperatury. Podsumowując badania 
dielektryczne należy podkreślić, że niezwykłe właściwości dielektryczne posiada już 
sam haloizyt — naturalny minerał. Chemiczne próby modyfikacji minerału nie pro- 
wadzą do jakościowych zmian, okazuje się jednak, że ciekawe właściwości relaksa- 
cyjne są uwydatnione i mogą być obserwowane w większym przedziale temperatur. 
Celem modyfikacji było wprowadzenie do przestrzeni międzypakietowej haloizytu 
cząsteczek zdolnych do polaryzacji dipolowej, a przez to zwiększających rzeczywi- 
stą i urojoną składową przenikalności. Czy się to udało? W kilku przypadkach tak. 
Znaczaco wyższe wartości przenikalności otrzymano dla interkalatu z formamidem 
i imidazolem, zaś wyraźna relaksacja dielektryczna była obserwowana w proszko- 
wej próbce formamidu. Interesujący i ważny z punktu widzenia zastosowań był 
interkalat z octanem potasu. Dla niego otrzymane wartości były szczególnie duże. 
Nie wiąże się tego w sposób jednoznaczny z polaryzacją elektryczną, choć i ten 
mechanizm jest możliwy do rozważenia. Przypuszcza się jednak, że zasadniczą rolę 
odgrywa higroskopijność octanu potasu, co powoduje znaczny wzrost przewod- 
nictwa powierzchniowego ziaren, przewodnictwa stałoprądowego oraz duże efekty 
polaryzacyjne. Obok celów poznawczych, prezentowana praca miała również cel 
praktyczny, zmierzający do wskazania efektywnych i tanich absorberów promienio- 
wania. 
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Rysunek 11. Współczynnik tłumienia (neper/m) w funkcji częstotliwości: 1 — haloizyt, 2 - HA-NMEFA, 3 - HA- 
ACOK, 4 - suchy haloizyt, 5 - HA-FA-IMI, 6- HA-FA, 7- HA-UREA, 8 - HA-AA, 9 - HA-NMFA- 
AKA, 10 - HA-NMFA-EDTA. Oznaczenia takie jak na Rysunku 3 

Figure 11. Attenuation coefficient (neper/m) in function of frequency: 1 - halloysite, 2 - HA-NMFA, 3 - HA- 
ACOK, 4 — dry halloysite, 5 - HA-FA-IMI, 6- HA-FA, 7- HA-UREA, 8 - HA-AA, 9 - HA-NMFA- 
AKA, 10 - HA-NMFA-EDTA. Captions the same as in Figure 3 
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Tłumienie elektrycznej składowej promieniowania zależy zarówno od rzeczy- 
wistej, jak i urojonej składowej przenikalności. Korzystając z rezultatów przepro- 
wadzonych badań obliczono współczynnik tłumienia dla haloizytu i wszystkich 
badanych interkalatów, biorąc pod uwagę doświadczalne wartości przenikalności 
dla temperatury ok. 25°C. Rysunek 11 przedstawia rezultaty obliczeń. 


UWAGI KOŃCOWE 


Otrzymane rezultaty uczą nas pokory wobec Natury. Największy współczynnik 
tłumienia posiada naturalny haloizyt zawierający słabo związaną wodę wewnątrz 
przestrzeni międzywarstwowych. Spośród przebadanych modyfikacji haloizytu naj- 
lepsze pod względem tłumienności są interkalaty zawierające N-metyloformamid. 
Porównanie przedstawione na Rysunku 11 dotyczy haloizytu przemywanego wodą 
i stosunkowo krótko suszonego na powietrzu. Dłuższe suszenie, w bardziej agre- 
sywnych warunkach, drastycznie zmniejsza tłumienność minerału. W odróżnieniu 
od haloizytu, właściwości tłumienne interkalatów są bardziej stabilne. Ogrzewane 
nawet do umiarkowanie wysokich temperatur (40-50°C) nie zmniejszają swojej tłu- 
mienności w tak dużym stopniu jak minerał naturalny. Biorąc to pod uwagę, otrzy- 
mane interkalaty wydają się być obiecującymi materiałami do konstrukcji tanich 
absorberów promieniowania elektromagnetycznego. 
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PODZIĘKOWANIA 


„Człowiek jest wielki nie przez to co posiada, 
lecz przez to, kim jest, 

nie przez to, co ma, 

lecz przez to, czym dzieli się z innymi” 

JAN PAWEŁ II 


Komitet Redakcyjny „Wiadomości Chemicznych” składa najserdeczniejsze 
podziękowania Jego Magnificencji Rektorowi Uniwersytetu Wrocławskiego — Panu 
Profesorowi Adamowi Jezierskiemu oraz Dziekan Wydziału Chemii Uniwersytetu 
Wrocławskiego — Pani Profesor Annie Trzeciak za wsparcie finansowe udzielone 
na wydanie zeszytu poświęconego Profesorowi Aleksandrowi Kollowi z okazji 
80. rocznicy Jego urodzin. 


Komitet Redakcyjny „Wiadomości Chemicznych” składa również serdeczne 
podziękowania Panu Profesorowi Kazimierzowi Orzechowskiemu za ogromne 
zaangażowanie w powstanie zeszytu poświęconego Profesorowi Aleksandrowi Kol- 
lowi z okazji 80. rocznicy Jego urodzin. 


prof. dr hab. Zdzisław Latajka 
Redaktor Naczelny 
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MECENASI WIADOMOSCI CHEMICZNYCH 


* Hydrolab to polska firma z ponad 20-letnim doświadcze- 
hs niem na rynku. Jestesmy dobrze zorganizowanym, nowocze- 
EES śnie zarządzanym przedsiebiorstwem. 


Produkujemy laboratoryjne demineralizatory, zaprojek- 
towane zgodnie z wytycznymi polskich i europejskich norm, wykorzystując naj- 
nowsze technologie oczyszczania wody. 

Jesteśmy w stanie zaplanować całą gospodarkę wodną w laboratorium, z pełną 
dokumentacją kwalifikacyjną. 

Hydrolab aktualnie posiada w swojej ofercie ponad sto modeli urządzeń do 
oczyszczania wody. 


Selwa jest firmą założoną przez dobrze znanych Państwu 
SLLWA ludzi, działających w branży sprzętu laboratoryjnego i apara- 
E, tury naukowo-badawczej produkowanej przez światowych 

liderów. 

Głównym celem firmy jest zapewnienie Państwu profesjonalnego doradztwa 
w zakresie doboru aparatury, optymalizacji konfiguracji i wykorzystania dokładnie 
zgodnie z Państwa potrzebami. 

Nasze ponad 15-sto letnie doświadczenie w tej branży pozwala dopasować 
naszą ofertę dokładnie do Państwa potrzeb. 

Świadczymy również jako jedna z niewielu firm na rynku, kompletne rozwiąza- 
nia łącznie z wdrożeniem metod pod konkretne Państwa wymagania. Sprzęt przez 
nas sprzedany będzie, po skończeniu wdrożenia, gotowy do pracy bez konieczności 
wykonywania dodatkowych działań z Państwa strony i poświęcania Waszego cen- 
nego czasu. 

Z naszymi partnerami zawarliśmy umowy na zasadach wyłączności i upraw- 
niające nas do sprzedaży, serwisowania a przede wszystkim wsparcia technicznego 
i aplikacyjnego oferowanego sprzętu na terenie Polski. 

Oferowany przez naszą firmę aparatura posiada certyfikaty zgodności CE, 
spełnia warunki GLP oraz GMP i jest wytwarzana przez producentów, którzy mają 
wdrożone systemy zarządzania jakością ISO 9000. 
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ZMARŁ PROF. DR HAB. INŻ. JACEK NAMIESNIK, 
REKTOR POLITECHNIKI GDAŃSKIEJ 


F ké Z Politechniką Gdańską związany był przez całą 
I swoją karierę zawodową. W 1967 roku rozpoczął studia 

na Wydziale Chemicznym. Od 1995 roku kierował Kate- 

ka A drą Chemii Analitycznej, jednego z najprężniejszych 


i najlepiej rozpoznawalnych na arenie międzynarodo- 
wej polskich ośrodków naukowych w zakresie chemii analitycznej i środowiskowej. 
Prof. Namieśnik był jednym z najbardziej aktywnych polskich chemików, zarówno 
pod względem naukowym, dydaktycznym, kształcenia młodej kadry, jak i organiza- 
cji życia akademickiego. Wypromował 66 doktorów (14 jako współpromotor). 

Był specjalistą z zakresu chemii analitycznej i chemii środowiska. Jego dorobek 
cieszy się dużym zainteresowaniem w środowisku naukowym. To autor lub współ- 
autor ponad 800 prac opublikowanych w czasopismach z listy JCR oraz około 200 
prac, które ukazały się w innych czasopismach. Prof. Namieśnik ma na swoim kon- 
cie 10 przyznanych patentów, w tym jeden europejski. 

Na macierzystym Wydziale pełnił funkcję prodziekana ds. kształcenia 
oraz dwukrotnie funkcję dziekana. 

Od 2016 roku Rektor Politechniki Gdańskiej oraz wiceprzewodniczący Konfe- 
rencji Rektorów Polskich Uczelni Technicznych. 

Za swoją działalność wielokrotnie nagradzany i odznaczany. Dwukrotnie uho- 
norowany tytułem doktora honoris causa, Nagrodą Naukową Miasta Gdańska im. 
Jana Heweliusza, licznymi Nagrodami Prezesa Rady Ministrów oraz Ministra Nauki 
i Szkolnictwa Wyższego. Odznaczony Krzyżem Komandorskim Orderu Odrodze- 
nia Polski i innymi najwyższymi odznaczeniami Panstwowymi. 

Prof. Jacek Namieśnik był również wiceprzewodniczącym Konferencji Rekto- 
rów Polskich Uczelni Technicznych (KRPUT). Należał do Komitetu Chemii Anali- 
tycznej PAN, a przez dwie kadencje pełnił funkcję przewodniczącego tegoż komi- 
tetu. W latach 2007-2016 udzielał się jako członek Centralnej Komisji do spraw 
Stopni i Tytułów. 

Urodził się 10 grudnia 1949 r. w Mogilnie. Prywatnie mąż Krystyny, ojciec 
dwóch córek: Katarzyny i Justyny, dziadek dwojga wnucząt. Interesował się historią 
Francji, chętnie jeździł na rowerze. 
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INFORMACJE REDAKCJI „WIADOMOŚCI CHEMICZNYCH” 


Redakcja czasopisma „Wiadomości Chemicznych” zwraca się do wszystkich 
Autorów z prośbą o dołączanie do zgłaszanych prac, numerów ORCID (dla wszyst- 
kich autorów danej pracy). Numery te będą zamieszczane przy Państwa nazwiskach 
w tekstach, które będą publikowane. 


CENY PRENUMERATY NA ROK 2019 


Redakcja miesięcznika PTChem „Wiadomości Chemiczne” zawiadamia, że 
wysokość prenumeraty rocznej „Wiadomości Chemicznych” za 2019 r. będzie wyno- 
siła 252 zł dla instytucji i niezrzeszonych prenumeratorów indywidualnych. Dla 
członków PTChem 20 zł. Należność za prenumeratę prosimy przekazywać na konto: 


PTChem Warszawa, ul. Freta 16 
Bank BGZ 
54 2030 0045 1110 0000 0261 6290 


Prenumerata „Wiadomości Chemicznych” dla członków PTChem, połączona 
z opłatą składek członkowskich, jest znacznie niższa i przedstawia się następująco: 
- prenumerata „Wiadomości Chemicznych” na rok 2019 wraz ze składką 
członkowską, wynosi 70 zł (składka — 50 zł, prenumerata — 20 zł); 
- emeryci, doktoranci oraz studenci płacą 35 zł (składka — 15 zł, prenume- 
rata — 20 zł); a nauczyciele szkół średnich i podstawowych płacą 40 zł 
(składka — 20 zł, prenumerata — 20 zł). 


Członkowie PTChem, którzy zechcą zaprenumerować „Wiadomości Chemiczne” 
na podanych tu warunkach, proszeni są o wnoszenie opłat na konto: 


PTChem Warszawa, ul. Freta 16 
Bank BGŻ 
54 2030 0045 1110 0000 0261 6290 
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XIII Kopernikańskie Seminarium Doktoranckie 
16 - 18 czerwca 2019, Bachotek k. Brodnicy 


1. Tytuł konferencji 

XIII Kopernikańskie Seminarium Doktoranckie (XIII KSD; www.ksd.umk.pl) 

2. Termin i miejsce 

16-18.06.2019r., miejsce: Ośrodek Konferencyjno-Wypoczynkowy Uniwersytetu Mikołaja 
Kopernika, Bachotek, k. Brodnicy 

3. Dane organizatora 

Wydział Chemii Uniwersytetu Mikołaja Kopernika w Toruniu 

Przewodniczący Komitetu Organizacyjnego IX KSD, Dziekan Wydziału Chemii UMK 

— prof. dr hab. Edward Szłyk 

Komitet organizacyjny 

. prof. dr hab. H. Kaczmarek - Kierownik Studiów Doktoranckich WCh UMK 

. dr hab. I. Łakomska, prof. UMK 

. dr hab. U. Kiełkowska 

. dr hab. J. Scianowski 

. dr hab. J. Nowaczyk 

. dr hab. A. Nowaczyk 

. dr A. Filipiak-Szok 

. dr A. Kmieciak 

. mgr S. Grabska 

. mgr S. Kowalska 

. mgr N. Stachowiak 

4. Zasięg konferencji 

XIII Kopernikańskie Seminarium Doktoranckie (XIII KSD) ma zasięg ogólnopolski. Każdego 
roku przyjeżdżają do nas doktoranci z Uczelni oraz Jednostek Badawczych z całej Polski 
(podczas X KSD w 2016 roku gościliśmy doktorantów m.in. z UMK, CM UMK w 
Bydgoszczy, UKW, UTP w Bydgoszczy, UW, PW, WAT w Warszawie, AGH, UJ w 
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Krakowie, UŁ, PŁ w Łodzi, UG, GUM w Gdańsku, UŚ w Katowicach, UW, UM we 
Wrocławiu, UP, UM w Poznaniu, UP w Lublinie, UWM w Olsztynie, ZUT, PUM w 
Szczecinie, UJK w Kielcach, PCz w Częstochowie, UZ w Zielonej Górze oraz Instytuty PAN z 
Warszawy, z Sopotu, z Lublina). 


5. Przewidywana liczba uczestników 


Przewidywana liczba Uczestników, zaproszonych Gości, członków Jury, Kadry Naukowej, 
Komitetu Organizacyjnego w liczbie 200 osób). 


6. Charakter konferencji (cykliczna, jednorazowa) 


konferencja charakter cykliczny. W tym roku organizowana jest już po raz trzynasty, a swój 


początek ma 2008 roku. 


7. Opis merytoryczny 

XIII Kopernikańskie Seminarium Doktoranckie (XIII KSD) jest skierowane do doktorantów, 
młodych doktorów, przedstawicieli świata nauki reprezentujących wszystkie uczelnie wyższe, 
jednostki PAN z całej Polski, a także Europy. XIII KSD ma charakter interdyscyplinarny i 
kierowane jest do młodych adeptów nauki realizujących swoje badania w zakresie chemii, fizyki i 
nauk technicznych, farmacji, medycyny, biologii i nauk przyrodniczych. Podczas Seminarium 
uczestnicy prezentują wyniki w formie komunikatu ustnego lub posteru na trzech równolegle 
trwających sekcjach: nauk chemicznych, biologiczno-medycznych i techniczno-fizycznych. 
Wszystkie wystąpienia są oceniane przez specjalistów z danej dziedziny, którymi są pracownicy 
naukowi wydziałów naszego Uniwersytetu oraz innych szkól wyższych. Najlepsze wystąpienia 
ustne oraz postery są nagradzane przez komisje oceniające dane dziedziny wiedzy. Dzięki temu 
spotkaniu naukowcy nie tylko mogą zaprezentować wyniki badań własnych i zapoznać się z 
tematyką badawczą kolegów, ale także jest to wyjątkowa okazja do wymiany pomysłów, 
wspólnych dyskusji i nawiązania współpracy międzyuczelnianej. 
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XIII pozońka cake Seminarium Doktoranckie 
16 - 18 czerwca 2019, Bachotek k. Brodnicy 


KSD odbywa się pod patronatem honorowym J.M. Rektora Uniwersytetu Mikołaja Kopernika 
- prof. dr hab. A. Tretyna, Dziekana Wydziału Chemii UMK - prof. dr hab. E. Szłyka, 
Prezydenta Miasta Torunia - Michała Zaleskiego, Marszałka Województwa Kujawsko- 
Pomorskiego — P. Całbeckiego, PTChem Oddział Toruński, SITPChem. 


8. Wnioskowana forma wsparcia: 


- wsparcie finansowe; 
- przekazanie materiałów promujących firmę (m.in. notatniki, długopisy, teczki, torby, itp.); 
- ufundowanie nagród przyznanych za najlepsze wystąpienia ustne i posterowe. 


W zamian za wsparcie finansowe zaoferować możemy: 


- umieszczenie logo firmy/insynuacji Sponsora i Patrona na stronie internetowej Konferencji 
(www.ksd.umk.pl), w broszurach i plakatach informacyjnych, w materiałach konferencyjnych; 
- możliwość uczestnictwa przedstawicieli firmy w uroczystym rozpoczęciu bądź zakończeniu 
Konferencji, lub podczas trwania Konferencji (stoisko firmowe); 

- umieszczenie dostarczonego baneru, plakatu w Sali Obrad na Wydziale Chemii; 


- umieszczenie nazwy Patrona i Sponsora w pisemnym sprawozdaniu z konferencji 
Jednocześnie zobowiązujemy się do przedstawienia materiałów konferencyjnych, broszur, 
plakatów, itp. oraz sprawozdań, w których ukaże się relacja z konferencji. 


Łączę wyrazy szacunku 


Prof. dr hab. Edward Szłyk 
Dziekan Wydziału Chemii UMK 
Przewodniczący Komitetu Organizacyjnego XIII KSD 


BIURO KONFERENCJI: 
UNIWERSYTET MIKOŁAJA KOPERNIKA W TORUNIU 


Wydział Chemii ` UNIWERSYTET 
ul. Gagarina 7, 87-100 Torun MIKOŁAJA KOPERNIKA 
e-mail: ksdQumk.pl W TORUNIU 

© 


www.ksd.umk.pl 
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Biochemia 


Jeremy M. Berg * Lubert Stryer Autor: 


John L. ymoczko * Gregory J. Gatto Jeremy M. Berg, Lubert Stryer, 


BIOC H È M IA John L. Tymoczko, Gregory J. Gatto 


Wyd. V / październik 2018 / s. 1118 
Typ oprawy: miękka 
Cena: 249 zł 


Wydawca: Wydawnictwo Naukowe 
PWN 


Najpopularniejszy podręcznik biochemii na świecie! Najnowsza wersja podręcznika 
zawiera ogrom informacji dotyczących biochemii, genetyki i fizjologii organizmów, 
ze szczególnym uwzględnieniem procesów zachodzących w organizmie człowieka. 
Autorzy jeszcze bardziej skupili się na klarownym i sprawnym przedstawianiu zagad- 
nień biochemicznych oraz omówili ważne dla biochemii nowe odkrycia. W całym 
podręczniku jaśniej zdefiniowano wiele podstawowych pojęć wzmacniając ich opisy 
licznymi zastosowaniami w nowych badaniach naukowych. Każdy z 36. rozdziałów 
rozpoczyna się krótkim spisem treści danego rozdziału, a kończy podsumowaniem 
omawianego materiału. Poszczególne rozdziały zawierają zadania umożliwiające 
samoocenę zrozumienia zagadnień, testy wielokrotnego wyboru, testy, w których 
udziela się krótkich odpowiedzi, testy dopasowania, zadania problemowe oraz roz- 
wiązania wszystkich typów zadań. Przystępny tekst i porządek przedstawianych 
zagadnień uzupełniają rysunki obrazujące poszczególne zagadnienia. 
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Poligon rachunkowy dla chemikow. 
Renata Jastrząb Zbiór zadań z podstaw chemii 
Romualda Bregier-Jarzebowska 


Małgorzata T. Kaczmarek Autor: Małgorzata Kacz- 
Ce marek, Renata Jastrzab, 


POLIGON oe 


ra ( h U n kowy Wyd. I / sierpień 2018 / s. 218 
dla chemików Typ oprawy: miękka 


Cena: 69 zł 


ZBIÓR ZADAŃ Z PODSTAW CHEMII 


Wydawca: Wydawnictwo Naukowe PWN 


Ta książka to niezbędnik, który pozwoli przetrwać studentom pierwszych lat stu- 
diów chemicznych oraz kierunków pokrewnych. Jego celem jest ugruntowanie 
nabytej już wiedzy, jak również wprowadzenie zagadnień, które nie są realizowane 
w szkole średniej, a w wielu przypadkach są niezbędne w pracy chemika. Materiał 
zebrany w zbiorze obejmuje zagadnienia z: 

- obliczeń bazujących na podstawowych prawach chemicznych, 

- stężeń oraz ich przeliczania, 

- mieszania roztworów, 

- rozpuszczalności, 

- pisania równań reakcji redukcji i utleniania oraz 

- zadań opartych o równania reakcji. 

Wszystkie zadania poprzedzone są krótkim wstępem teoretycznym wraz z szeroko 
opisanymi przykładami. 
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Chemia organiczna. 
Marcin Kazmir Testy egzaminacyjne z rozwigzaniami 
Katarzyna Koroniak:Szejn Autor: Marcin Kaźmierczak, 
Tomasz Cytlak, 


© h = m i A Katarzyna Koroniak-Szejn, Hen- 


ryk Koroniak 
Wyd. I / sierpień 2019 / s. 200 


Testy egzaminacyjne T r: : mi kk 
z rozwiązaniami HN op awy E a 
Cena: 59 zł 


Wydawca: Wydawnictwo Naukowe PWN 


| ==" | Bo 


Publikacja zawiera testy egzaminacyjne wraz z metodologia ich rozwiazywania. 
Prowadzi czytelnika krok po kroku, tak, by nauczyt sie korzystania z uzyskanej 
wiedzy przy rozwiązywaniu konkretnych problemów. Książka skierowana do stu- 
dentów kierunków chemicznych wszystkich uczelni wyższych. Może być pomocą 
dydaktyczną dla wykładowców chemii organicznej. 
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Chemia obliczeniowa 
Autor: Jeremy Harvey 


Wyd. I/ lipiec 2019 / s. 220 
Typ oprawy: miękka 
Jeremy Harvey Cena: 69 zł 


Wydawca: Wydawnictwo Naukowe PWN 


CHEMIA 


obliczeniowa 


PWN 


Publikacja wprowadza studentów w różnorodne komputerowe metody oblicze- 
niowe. Zawiera liczne przykłady i aplikacje pokazujące znaczenie obliczeń kompu- 
terowych dla badań i przemysłu chemicznego, pytania na końcu każdego rozdziału, 
rozwiązania do ćwiczeń oraz różnorodne zastosowania chemii obliczeniowej zade- 
monstrowane poprzez liczne przykłady. 


Redakcja „Wiadomości Chemicznych” informuje, ze są u nas do nabycia 
następujące pozycje „Biblioteki Wiadomości Chemicznych”: 


Nomenklatura chemii nieorganicznej. Zalecenia 1990, red. Z. Stasicka, cena 25 zł 

Podstawowa terminologia stereochemii oraz Słownik podstawowych terminów 
w nauce o polimerach. Zalecenia 1996, red. O. Achmatowicz, B. Szechner 
i P. Kubisa, cena 12 zł 

Nomenklatura węglowodanów. Zalecenia 1996, tłum. i red. T. Sokołowska 
i A. Wiśniewski, cena 18 zł 

I.Z. Siemion, Bronisław Radziszewski i lwowska szkoła chemii organicznej, cena 


18 zł 

K. Maruszewski, Fizykochemia molekuł zamkniętych w zeolitach i zol-żelach, 
cena 18 zł 

Praca zbiorowa, Uporządkowane materiały mezoporowate, red. B. Burczyk, 
cena 18 zł 


Skorygowana nomenklatura rodników, jonów, jonorodników i podobnych indywi- 
duów chemicznych. Zalecenia 1993, red. T. Sokołowska i A. Wiśniewski, 


cena 15 zł 

I.Z. Siemion, Lutum sapientiae, czyli Notatek chaotycznych część pierwsza, 
cena 18 zł 

M. Zabłocka-Malicka, Ruchliwość jonów w podwójnych układach stopionych soli, 
cena 8 zł. 


Praca zbiorowa, Nanomateriały, red. D. Hreniak, W. Łojkowski, W. Stręk, 
M. Suszyńska, cena 25 zł. 

Praca zbiorowa, Ogniwa paliwowe — nowe kierunki rozwoju, red. H. Drulis, 
J. Hanuza, D. Hreniak, M. Miller, G. Paściak, W. Stręk, cena 20 zł 

Glosariusz nazw klas związków organicznych i reaktywnych produktów pośrednich 
oparty na strukturze (Zalecenia IUPAC 1994), red. i tlum. F. Kaźmierczak 
i J. Gawroński, cena 16 zł. 

Od substacji prostych do życia. Świat RNA — początki życia na Ziemi, Zdzisław 
Chilmonczyk - NAKŁAD WYCZERPANY. 

Profesor Bogusława Jeżowska-Trzebiatowska. 1908-1991 w setną rocznicę uro- 
dzin, cena 12,00 zł. 

Chemia koordynacyjna w Polsce. Część I. Nakład wyczerpany, dostępna wersja 
elektroniczna. 

Chemia koordynacyjna w Polsce. Część Il. Nakład wyczerpany, dostępna wersja 
elektroniczna 

Chemosensory optyczne oraz materiały rozpoznawcze dla jonów metali w roztwo- 
rach, Krzysztof Kledzik, cena 22,00 zł. 

Obliczenia teoretyczne stałej ekranowania magnetycznego i stałych strzężeń spi- 
nowo-spinowych. Teobald Kupka cena 20,00 zł. 


Książki wysyłamy na koszt zamawiającego. Zamówienia prosimy kiero- 
wać pod adresem: Redakcja „Wiadomości Chemicznych”, ul. F. Joliot-Curie 14, 
50-383 Wrocław lub na e-mail: czasopisma@ptchem.pl. Opłaty należy wnosić na kon- 
to: PTChem Warszawa, ul. Freta 16, Bank BGZ 54 2030 0045 1110 0000 0261 6290. 


REGULAMIN I INFORMACJE DLA AUTOROW PUBLIKUJACYCH 
W CZASOPIŚMIE „WIADOMOŚCI CHEMICZNE” 


1. Informacje ogólne 

„Wiadomości Chemiczne” są recenzowanym czasopismem naukowym Polskiego Towarzystwa Chemicznego, 
które publikuje przede wszystkim artykuły przeglądowe. Ponadto publikowane są tutaj inne wartościowe mate- 
riały o charakterze edukacyjno-informacyjnym takie jak: artykuły oparte na pracach doktorskich lub habilitacyj- 
nych, które zostały wyróżnione przez Rady Wydziałów, przed którymi toczyły się odpowiednie procesy; mate- 
riały informacyjne na temat uczonych oraz jednostek naukowych/firm chemicznych lub pokrewnych chemii; 
materiały o aktualnych osiągnięciach w szeroko pojętych naukach chemicznych. 

Dodatkową ofertę Wydawnictwa stanowi seria, „Biblioteka Wiadomości Chemicznych gdzie publikowane są 
dłuższe artykuły przeglądowe lub monografie poświęcone ważnym i aktualnym problemom współczesnej che- 
mii. Autorzy, którzy chcieliby takie prace napisać, powinni wcześniej skontaktować się z Redakcją, a następnie 
przesłać wstępnie przygotowaną publikację (redagowaną na wzór artykułów w czasopiśmie „Wiadomości Che- 
micznych”) lub informację na temat przygotowywanej pracy - tytuł przygotowywanej publikacji, przybliżoną 
liczbę stron, tabel, rysunków. W chwili obecnej Redakcja nie posiada środków na finansowanie prac w serii 
„Biblioteka Wiadomości Chemicznych”. W zależności od sytuacji finansowej Wydawnictwa, Redakcja zastrzega 
sobie prawo negocjacji kosztów druku z autorami lub Instytucjami zlecającymi druk. 

„Wiadomości Chemiczne” wydawane są zarówno w wersji drukowanej jak i elektronicznej. Wersja elektro- 
niczna udostępniana jest bezpłatnie w Internecie. 

Czasopismo jest indeksowane/abstraktowane w kilku bazach danych: Chemical Abstracts, Polska Bibliografia 
Naukowa, BazTech, Polska Bibliografia Lekarska, Index Copernicus, Baza ARIANTA. 


2. Informacje dla autorów na temat wymagań i zasad publikowania prac 

e Prace nie były wcześniej publikowane, ani nie są złożone w redakcji innego czasopisma. 

e Autorzy prac stosują się do wymagań praw autorskich tzn. w przypadku zamieszczania rysunków, tabel 
itp., pochodzących z opracowań opublikowanych w innych czasopismach lub publikacjach zwartych, 
posiadają pisemną zgodę na ich przedruk. 

e Opublikowana raz praca bez zgody Redakcji, nie może być wydawana gdzie indziej. 

e Autorzy przysyłający prace po raz pierwszy powinni podać swój numer telefonu oraz adresy poczty tra- 
dycyjnej i elektronicznej. Jest to niezbędny warunek sprawnego przebiegu opracowania redakcyjnego 
tekstu. 

e Autorzy zobowiązani są do wykonania korekty tekstu. W pracach przyjętych do druku Redakcja ma 
prawo dokonywania niezbędnej korekty. 

» Jeżeli autorzy nie zastrzegą inaczej w momencie zgłoszenia pracy, wydawca nabywa ogólnych praw autor- 
skich do wydrukowanych prac (w tym prawo wydawania na nośnikach elektronicznych oraz w Interne- 
cie). Tytułem powyższego wykorzystania utworów autorom nie są wypłacane honoraria. 

e Wszystkie nadsyłane prace są poddawane wstępnej ocenie, która określa czy odpowiadają randze i pro- 
filowi „Wiadomości Chemicznych” oraz czy zostały przygotowane zgodnie z formalnymi wymogami 
MNiSW oraz Redakcji. 

e Po uzyskaniu pozytywnej wstępnej oceny wszystkie prace są recenzowane przez co najmniej dwóch nie- 
zależnych recenzentów, zgodnie ze wskazówkami zawartymi w broszurze informacyjnej Ministerstwa 
Nauki i Szkolnictwa Wyższego, http://www.nauka.gov.pl/fileadmin/user_upload/ministerstwo/Publika- 
cje/20110216_MNISW_broszura_210x210.pdf. 

e O przyjęciu pracy do druku decyduje Komitet Redakcyjny. 

e Prace, które Komitet Redakcyjny na podstawie uzyskanych recenzji stwierdził, że nie należy przyjąć do 
druku w czasopiśmie, po uwzględnieniu sugestii recenzentów mogą być powtórnie przesłane do czaso- 
pisma. W takim przypadku praca traktowana jest jako nowy tekst i ponownie przechodzi pełną proce- 
durę recenzowania. 

e Ponadto Komitet Redakcyjny informuje, że tzw. „ghostwiting” (ktoś wniósł znaczący wkład w powstanie 
publikacji, a nie został przedstawiony jako współautor lub też nie został wymieniony w podziękowaniu 
zamieszczonym w publikacji) lub „guest authorship” (udział autora jest znikomy lub też w ogóle nie miał 
miejsca, a mimo to jest współautorem publikacji) są przejawem nierzetelności naukowej. Wszelkie prze- 
jawy nierzetelności naukowej, łamania i naruszania zasad etyki obowiązującej w nauce będą ujawniane, 
włącznie z powiadomieniem jednostek zatrudniających autorów. 

e Autorzy mają prawo do zaproponowania co najmniej trzech niezależnych recenzentów, jednak ostatecz- 
nego wyboru anonimowych recenzentów dokonuje Redakcja. 


3. Koszty 

Autorzy czasami mogą ponosić częściowe koszty wydania swoich artykułów. Tak jest w przypadku tzw. stron 
nadliczbowych tj. powyżej 25 stron. Za każdą rozpoczętą nadliczbową stronę jest naliczana opłata w wysokości 
około 50 zł. Najczęściej kwota ta pokrywana jest z funduszy pozyskiwanych przez Autorów lub przez Wydziały 
które wspomagają wydawanie „Wiadomości Chemicznych”. Niezależnie od rodzaju pracy opłata pobierana jest 
również za strony drukowane w kolorze (zgodnie z faktycznym kosztem druku). 

Redakcja zastrzega sobie możliwość zmiany wysokości opłat, w zależności od wielkości dofinansowania 
z MNiSW oraz wypracowanych środków własnych. Faktura wystawiana jest po ukazaniu się pracy. 

W przypadku prac w serii „Biblioteka Wiadomości Chemicznych”, Redakcja nie posiada środków na finanso- 
wanie i zastrzega sobie prawo negocjacji kosztów druku z autorami lub Instytucjami zlecającymi druk. 


4. Informacje szczegółowe dotyczące przygotowania maszynopisu do druku 

4.1. Wymagania merytoryczne 

Tekst należy napisać zwięźle, prostym stylem, według zasad pisowni polskiej, z zachowaniem poprawnego i 
obowiązującego nazewnictwa fachowego. Nie należy zamieszczać nadmiaru szczegółów odsyłając Czytelnika do 
piśmiennictwa oryginalnego, które to powinno uwzględniać najnowsze informacje, dotyczące napisanej pracy. 
Literaturę należy cytować ze źródeł oryginalnych. 

4.2. Wymagania techniczne składu tekstu 

e W przypadku prac współfinansowanych przez autorów, liczba stron oraz forma kolorystyczna manu- 
skryptu nie jest ograniczona (wymagane jest wcześniejsze uzgodnienie z Redakcją). 

e Maszynopisy prac autorów którzy nie chcą ponosić dodatkowych kosztów, nie powinny przekraczać 25 
stron całej pracy (po wydruku w czasopiśmie) oraz drukowane będą w wersji czarno białej. 

e Główny tekst nadsyłanych prac powinien być napisany w edytorze Word, czcionką Times New Roman, 
12p z zachowaniem interlinii 1,5 oraz z 5 cm marginesem z prawej strony. Przy podziale tekstu należy sto- 
sować numerację cyfrową wielorzędową. Numerujemy tylko tytuły rozdziałów, nie numerujemy działów: 
Abstract, Wykaz stosowanych skrótów, Wprowadzenie, Uwagi końcowe, Podziękowanie, Piśmiennictwo 
cytowane. Jednolity sposób numeracji konsekwentnie stosuje się wewnątrz tekstu (w całym tekście tj. 
zarówno przy numerowaniu rozdziałów, przy przytaczaniu piśmiennictwa cytowanego oraz odwoływaniu 
się do tabel rysunków itp., nie należy stosować odsyłaczy hipertekstowych). 

e Tekst powinien być napisany poprawnym językiem, wszystkie skróty muszą być wyjaśnione, oznaczenia i 
jednostki miar należy podawać według układu SI, pozycje cytowanej literatury należy oznaczać numerami 
umieszczonymi w nawiasach kwadratowych, w kolejności cytowania wg wzorów [1, 5, 7] (dla prac 1, 5i7) 
lub [1-5, 7] (dla prac od 1 do 5 oraz pracy 7). 

e Jeśli w artykułach znajdują się przedruki rysunków, czy innych elementów prac cudzych, w opisach (pol- 
skich i angielskich) należy zamieścić stosowną informację. 

e Zaleca się umieszczać w tekście pracy rysunki, tabele oraz podpisy (jeśli są przygotowane w edytorze 
Word), jednak w przypadku plików o bardzo dużych rozmiarach należy zaznaczyć miejsca na ich umiesz- 
czenie (zob. Pliki jakie należy przekazać do Redakcji). 

e Pierwsza strona pracy powinna zawierać kolejno: 

- tytuł pracy w języku polskim (Times New Roman, 14 p, pogrubiony, WERSALIKI), i angielskim 
(Times New Roman, 14 p, WERSALIKI), 

— pełne imiona i nazwiska autorów wraz z numerami ORCID, 

— pełne imię i nazwisko autora (autorów) pracy (Times New Roman, 15p, pogrubione), 

— pełne nazwy ośrodków przypisane do autorów pracy (wraz z adresem ośrodka i adresem e-mail autora 
korespondującego (Times New Roman, 10,5, kursywa), 

— spis treści pracy z zastosowaniem następującego podziału: 

Abstract 

Wykaz stosowanych symboli i oznaczeń 

Wprowadzenie 

1. Tytul rozdziału 

1.1. Tytuł podrozdziału itp. 

Uwagi końcowe 

Podziękowanie 

Piśmiennictwo cytowane 

e Kolejne strony pracy powinny zawierać: 

- notki o autorach pracy wraz z tytułami naukowymi (można dołączyć osobno pliki z fotografiami 
autorów (zob. Pliki jakie należy przekazać do Redakcji), 


- obszerne streszczenie pracy w języku angielskim (od 1800 do 2700 znaków ze spacjami) z uwzględ- 
nieniem cytowanego piśmiennictwa oraz odsyłaczami do tabel, rysunków zamieszczonych w tekście 
(Rys. 1, Tab. 1-2, Schemat 1) oraz słowa kluczowe - nie więcej niż 6, uzyskane najlepiej z bazy haseł 
przedmiotowych podawane w języku angielskim i polskim, 

- wykaz stosowanych skrótów — w przypadku niewielkiej liczby skrótów lub akronimów nie jest 
konieczne zamieszczanie tej pozycji, wówczas, skróty wyjaśniamy w tekście przy pierwszym użyciu. 
Angielskie skróty należy podać i wyjaśnić wg poniżej podanego wzoru lub w oparciu o inne prace 
zamieszczone w „Wiadomościach Chemicznych”. Przykład: dla skrótu SSRI - selektywne inhibitory 
zwrotnego wychwytu serotoniny (ang. Selective Serotonin Reuptake Inhibitor), 

- dalszy tekst pracy zgodny z podawanym wcześniej spisem treści. 

e Tabele, rysunki, fotografie 

Tabele i rysunki powinny być zamieszczone w przesyłanym tekście oraz dodatkowo (po zatwierdzeniu 

pracy do druku, na etapie przygotowywania szczotki) dołączone w postaci osobnych plików zapisanych 

w formacie pdf , doc, jpg, tiff. 

Tabele i rysunki powinny być przejrzyste, zawierać informacje niezbędne do zrozumienia treści, bez 

konieczności poszukiwania objaśnień w tekście pracy, należy je numerować cyframi arabskimi oraz podać 

tytuł (polski/angielski, nad tabelą, pod rysunkiem, Times New Roman, 10 p). 

Wszystkie fotografie — należy przesłać w postaci plików zapisanych w formacie tif, jpg lub podobnym, 

każdą zapisać w oddzielnym pliku o rozdzielczości co najmniej 300 dpi. 

e Piśmiennictwo cytowane 

Piśmiennictwo należy zestawić numerycznie według kolejności cytowania w tekście, należy cytować 

wyłącznie pozycje istotne dla treści pracy w sposób precyzyjny. 

W przypadku artykułów z czasopism tradycyjnych, opis powinien zawierać kolejno następujące ele- 

menty: inicjały imion i nazwisko autora (autorów), skrót tytułu czasopisma zgodny z przyjętymi normami, 

rok wydania, numer wolumenu zaznaczony pogrubioną czcionką, numer pierwszej strony cytowanej 
pracy, np. 

[1] J. Kowalski, Wiad.Chem., 2007, 61, 473. 

[2] W. Kowalski, A. Nowak, Przem. Spoż. 2010, 51, 3. 

W przypadku książek najprostszy opis powinien zawierać: inicjały imion i nazwisko autora (autorów), 

tytuł książki, nazwę wydawcy, miejsce wydania, rok wydania, np. 

[1] J. Malinowski, Tytuł książki, PWN, Warszawa, 2004. 

[2] W. Kowalski, Tytuł książki, Volumed, Wrocław, 1999 

W przypadku zasobów Internetowych najprostszy opis powinien zawierać: inicjały imion i nazwisko 

autora (autorów), tytuł (artykułu) dokumentu online, [dostęp], wydawca, [data dostępu]. Warunki 

dostępu, np. 

[7] J. Kowalski, Tytuł artykułu. [online], wydawca, [dostęp: 2010-05-20]. Dostępny w Internecie: 


4.3. Materiały jakie należy przygotować w celu przesłania pracy do Redakcji 

Przed podjęciem decyzji o zakwalifikowaniu pracy do druku w celu oceny merytorycznej należy przesłać 
jeden plik kompletnej pracy zredagowany zgodnie z wymaganiami Redakcji. 

Po uzyskaniu pozytywnej recenzji i po ustosunkowaniu się do uwag Recenzenta oraz Redakcji należy prze- 
słać ostateczną wersję pracy w następującej postaci: 

e 1 plik tekstu zredagowany zgodnie z wymaganiami Redakcji; 

e 1 plik zawierający krótkie notki biograficzne o autorach nadesłanej pracy (każda notka do 150 wyra- 
zów powinna zawierać: tytuł naukowy, miejsce pracy oraz inne ważne informacje o autorze); 

e pliki zawierające zdjęcia portretowe autorów, w nazwie powinny wskazywać autora, którego zdjęcie 
dotyczy (dobrowolne, przesłanie plików jest jednoznaczne ze zgoda na jego opublikowanie); 

» 1 plik zawierający: stronę tytułową, streszczenie (abstrakt), słowa kluczowe, podpisy pod rysunki, 
tabele, schematy (wszystko w obu wersjach językowych); jeśli zachodzi potrzeba to również oddzielne 
pliki z rysunkami, schematami, tabelami (zob. Tabele, rysunki, fotografie). 

Prace nie odpowiadające wyżej wymienionym wymaganiom nie będą przyjmowane do druku. Redakcja 
zastrzega sobie prawo dokonywania poprawek stylistycznych i skrótów. Autorzy są zobowiązani do wykonania 
korekty artykułu i jego zwrotu do Redakcji w ciągu kilku dni od otrzymania. 

Na etapie przygotowania szczotki, w przypadku przesyłania prac z kolorowymi stronami prosimy o zaznacze- 
nie stron, które w formie druku mają być kolorowe Brak tej czynność będzie skutkował czarno-białym wydru- 
kiem wersji papierowej. W przypadku zmian w wersji drukowanej kolorowych rysunków na czarno-białe, 
prosimy o przesłanie dostosowanych do tego celu rysunków. 


Prace prosimy przesyłać pocztą elektroniczną na adres: czasopisma@ptchem.pl, zaś dokumenty wymagające 
podpisów autorów (np. list intencyjny, oświadczenia autorów, kopie zgody na przedruk potwierdzone za zgod- 
ność z oryginałem) pocztą tradycyjną na adres Redakcji. 
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